00:06 | So uh let me bring up your from last time. OK. So |
|
00:23 | first one is from um uh Li and she has um two questions, |
|
00:33 | concerning uh three questions, one con of them concerning presentation eight. |
|
00:40 | So she says slide 28 we consider core fluid pressure is uniform within the |
|
00:49 | , this implies seismic band waves, sonic or ultrasonic. And she |
|
00:55 | why is it not Sonic? So you know, that's a very good |
|
00:59 | . So let me ask you uh , what do you think is the |
|
01:03 | to that question? Yeah. She she doesn't know her voice is very |
|
01:09 | voice. And so I'm sure that folks uh uh online cannot hear |
|
01:14 | Uh So let me turn to Carlos. Uh why, why did |
|
01:19 | not include Sonic? Carlos? Are with me? Yeah. Yeah. |
|
01:28 | . Professor. I just um I thinking but I am OK. So |
|
01:31 | want you to think, I want to think out loud so we can |
|
01:34 | you thinking. Yeah. Yeah. . No, no, I I |
|
01:37 | know the answer. Professor. OK. So we want to consider |
|
01:41 | case where the um where the uh pressure is uniform within each max, |
|
01:49 | contains many grains. But it it's not as so big as a |
|
01:53 | , but it's uh uh uh uh the size of a hand sample has |
|
01:57 | grains. And we want to have uh consider the case where the pore |
|
02:03 | is uniform on that scale. So gonna mean if we're doing this with |
|
02:09 | uh waves of course waves are always . But uh if we consider low |
|
02:15 | uh frequency, then we can uh can, we can um trigger that |
|
02:21 | low frequency seismic waves satisfy the assumption local core pressure uniformity on every piece |
|
02:32 | the formation, you know, the of a hand sample. So uh |
|
02:37 | you know, uh it, it's good question because nobody really knows what |
|
02:41 | means. A me we always assume that includes uh uh uh uh seismic |
|
02:49 | . We always assume that low frequency uh it, it includes uh uh |
|
02:54 | seismic band. And so we can these results at the uh at the |
|
03:00 | for our seismic data. But think it, you know, doesn't that |
|
03:10 | implies something about permeability. If the is very high, like in a |
|
03:17 | , I, I guess it's probably that the seismic band is low enough |
|
03:22 | . But what happens if it's a , if, if it's a |
|
03:26 | you know, we have intrinsically low . And so maybe even in the |
|
03:33 | band for shales, um that's not enough frequency, maybe because of the |
|
03:40 | permeability uh during the passage of a seismic band wave in a |
|
03:47 | maybe that's not enough time for the to equalize everywhere. This is |
|
03:52 | a question, a problem which is uh well understood, hardly ever addressed |
|
03:58 | everybody asked that question. And so question is, is not quite |
|
04:03 | but it's, you know, it's . Now, let's turn to her |
|
04:07 | about Sonic. So Sonic means um of about 1000 cycles per second. |
|
04:14 | so, uh you can do some uh is that low enough or not |
|
04:21 | enough? And so when you make estimates, you're gonna assume something about |
|
04:25 | rotting. And uh it's uh uh think it's usually true that if uh |
|
04:32 | have a permanent rotate like in a sonic frequencies are too high to really |
|
04:39 | the, the uh um that the to really establish that uh uh ferocity |
|
04:50 | that the uh pore pressure is uniform every hand sample. So that's a |
|
04:54 | of a theoretical answer. Better answer come from our own rock physics lab |
|
05:01 | in um here at the University of . Uh uh uh maybe do you |
|
05:09 | uh long tang, you know? . So he's uh he's, he's |
|
05:15 | graduate student. Uh uh is he postdoc? He's a graduate student? |
|
05:21 | . And he's uh an advanced graduate and he's been doing a lot of |
|
05:25 | of uh experiments in our laboratory And uh he's using a piece of |
|
05:31 | which was developed by a, a uh student now graduated. And this |
|
05:38 | is uniquely capable of doing experiments through wide band of frequencies as low as |
|
05:46 | like 10 Hertz, I think. as high as um uh as high |
|
05:52 | 1000 Hertz, I think. And I'll, I'll hi his experiments. |
|
05:58 | , I'll verify, yeah, theoretical that I just gave for sandstones. |
|
06:06 | I'm not sure what he has learned shale, but I suspect that for |
|
06:11 | there might be a problem even so uh even a seismic frequencies just because |
|
06:18 | preme is so low uh as the is passing through that there's not time |
|
06:24 | for the fluid to squirt around inside rock to equalize out the core pressure |
|
06:29 | the, on the local scale. But uh that's uh you know, |
|
06:34 | are complicated and the best answer to like that is always experimental. And |
|
06:40 | I, I just uh gave you , the results of the experimentation. |
|
06:45 | let me uh now turn back to question um slide 35 and 36. |
|
06:53 | could you tell the gas mine equation the wrong assumption from that plot? |
|
06:59 | . So we've got to uh we've got to look at this, |
|
07:01 | got to look at that. I am going to um bring up |
|
07:10 | from, OK, here is um eight. Then I'm gonna go to |
|
07:38 | 35. OK. OK. So am now going to show um |
|
07:51 | uh I'm going to show this to goodly has. And so uh first |
|
07:58 | all, put it at presentation mode then I'll um, I hope go |
|
08:08 | the Z window and shoot her that and I'm gonna share a screen |
|
08:20 | OK. So you should be able see this now, I think. |
|
08:27 | . So uh uh this slide is taken from an early work by |
|
08:32 | 1985. Wow, that was a time ago. That was for 40 |
|
08:37 | ago, maybe before some of you were born. And so uh when |
|
08:43 | did this work, I didn't had idea how it was gonna come |
|
08:47 | Uh And so you can see we've here a lot of data, uh |
|
08:54 | probably uh 50 or 100 data Uh And they're all uh uh and |
|
09:01 | just uh uh the B let me my cursor. So it, it |
|
09:08 | the bolt modulus K as a function uh pressure for lots of rocks. |
|
09:14 | sort of the list of, of rocks. And uh I should |
|
09:20 | out that I did not do make of these measures. My colleagues at |
|
09:25 | did that and uh most of that already done. Uh uh before I |
|
09:31 | joined the company and they had uh the, the data, you |
|
09:36 | on uh pieces of paper in in the, uh uh in the |
|
09:45 | . We didn't really have computers uh in the way of computers in those |
|
09:51 | . Computers are very primitive, modern and we didn't have much in the |
|
09:56 | of databases. So, uh as matter of fact that when I first |
|
10:01 | the company, uh my boss OK, now we want you to |
|
10:05 | build uh to build the database out all this data. So we can |
|
10:12 | it out of the drawers and put in the computer so we can actually |
|
10:15 | it. That, that was a plan. That was my first uh |
|
10:19 | work at NFL. And out of came this project. So you see |
|
10:23 | their observations around this axis and all predictions from um uh from gas mon |
|
10:30 | on this axis. And if the were accurate, all the um observations |
|
10:36 | lie on this line. But you , they don't, they lie to |
|
10:40 | , to the right of this they lie at higher values of um |
|
10:45 | is observed to be higher values of mean faster velocities than Gasman predicted. |
|
10:53 | that's the uh uh the answer. how we can tell from this slide |
|
10:58 | Gas Mont was uh uh was wrong uh I shouldn't say it this |
|
11:05 | Uh uh I shouldn't say this proves is wrong. It proves that we |
|
11:11 | not apply Gas Mon's theory to ultrasonic . And the reason of course is |
|
11:18 | the ultrasonic data has some frequencies too . So that gas models assumption is |
|
11:23 | met. So this uh graph does prove the gas model is wrong. |
|
11:28 | just proves that we were, that were wrong to um apply his low |
|
11:34 | theory to high frequency data. let's look at the next slide. |
|
11:44 | . So uh this uh uh is single sample of Bria sandstone well understood |
|
11:51 | as a function of pressure from uh friend Arthur Chang here, the function |
|
11:56 | pressure here. And this is what measured and this is what the theory |
|
12:02 | . It should be. What's the based on uh uh based on uh |
|
12:07 | we measure um uh you measure the of uh the dry rock. And |
|
12:14 | we uh we measure the properties of fluid and we measure the amount of |
|
12:19 | porosity and we uh uh make sure properties of the core of are |
|
12:27 | And so, uh then on the of that, that's all you need |
|
12:30 | do for gas ma make this um prediction. And the prediction is way |
|
12:38 | . You see the prediction is um a function. The the error is |
|
12:44 | function of uh pressure. The error small here and it gets bigger and |
|
12:48 | at lower pressure. And so what means is that uh uh uh the |
|
12:55 | is at least in part due to closing of crack. So, if |
|
13:00 | look in gas man's theory, there's mention of cracks. His theory is |
|
13:04 | to apply to crack. And of , I guess man would have known |
|
13:16 | there would be cracks in some but he, he must have assumed |
|
13:21 | the effects of those cracks are um by the measurement of the dry |
|
13:29 | And uh but here you see, not true. So this is |
|
13:34 | this is a proof that uh that . Uh huh. Well, let's |
|
13:42 | . Now, I guess it's not . Oh, th this again proves |
|
13:48 | um uh we shouldn't be applying gasoline to ultrasonic. It could be that |
|
13:54 | crack, this effect of the cracks you see here only happens at ultrasonic |
|
14:01 | frequencies. And it, it could that if you did this uh uh |
|
14:07 | experiment on the use of sandstone rock proper low frequencies, the gas mont |
|
14:13 | be uh validate. Well, uh my knowledge that has not been done |
|
14:21 | we have recent theory that says that mine is wrong even at the lowest |
|
14:26 | . I think I talked with you that. So that answers uh Lena's |
|
14:33 | question. So now let's look at third question. Rich and um my |
|
14:43 | not. Oh And she asked her 35 what does the incidental reflector in |
|
14:52 | mean when the real part is So uh uh let me go to |
|
14:58 | slide 35 and you taste not like 35 it's not too far. |
|
15:42 | 34 35. OK. So, , let me put this into presentation |
|
15:56 | and then let me um share that the class particular stop sharing and then |
|
16:41 | gonna start sharing again. Yeah, time I always have this problem. |
|
17:08 | come over here and help me. , I'm not seeing that. |
|
17:15 | Thank you. You're reminded. Now, um uh so uh a |
|
17:26 | after the, the previous discussion leading to this slide, oops, let |
|
17:33 | get the pointer. After the previous , we answer to the slide. |
|
17:38 | , we uh discovered that the uh car efficient in the case of the |
|
17:44 | or have attenuation in the uh has an imaginary part to it. And |
|
17:50 | it shows explicitly and here's the real that you are accustomed to looking |
|
17:54 | And, and uh all the rest your professional life is this uh jump |
|
17:59 | uh impedance divide by twice. The impedance only thing is here. It's |
|
18:03 | it says uh um the real part specified that we're talking about the real |
|
18:10 | of that. Well, I'm sure whenever you looked at that normal |
|
18:16 | reflection coefficient before ever in your you've always just assumed it was |
|
18:21 | But now here we're saying, uh we're looking only at the real |
|
18:26 | if there's an imaginary part to it's not inside here because this is |
|
18:30 | the real part. And so now now what the sliders occur, we're |
|
18:35 | at when the real part goes to . And so what Lily is asking |
|
18:39 | what does it mean to have real ? Equal zero? Well, you |
|
18:43 | , that, uh that is just , a matter of imagining lithology contrast |
|
18:50 | that boundary uh where uh the, know, the fluids in the, |
|
18:56 | fluids are about the same and the are about the same and um uh |
|
19:03 | microstructure is about the same and everything about the same. So it's a |
|
19:07 | weak reflection, that's all, all means. And so uh uh uh |
|
19:11 | , I think it's easy to imagine um we have strong reflectors and weak |
|
19:17 | and this is just uh uh to the issues in your mind, it |
|
19:22 | imagine a case where we have a weak reflector and normal incidence. And |
|
19:28 | uh and when I say we, mean, the real parts of uh |
|
19:32 | , of OK velocity which go into are all are very similar on both |
|
19:39 | of the, of the reflector, that doesn't say anything about the |
|
19:47 | And so um uh the aeration is left when we assume this uh uh |
|
19:53 | is pretty small, we're still left the imaginary parts and uh that implies |
|
19:58 | the reflective wavelength uh is face shifted 90 degrees. Uh So what, |
|
20:05 | sort of Ortho could we have Uh uh uh You, you could |
|
20:10 | uh uh where it's uh uh sandstone e on either side of the, |
|
20:17 | the reflector and uh um, one is gas and one side is uh |
|
20:25 | brine. So you, you might thinking to yourself, well, how |
|
20:30 | that be? Uh surely the uh sides are sandstone then uh um uh |
|
20:38 | , the, the um gas is cross the boundary or the liquid is |
|
20:44 | cross the boundary. It's, it's homogenize. But think about this, |
|
20:49 | suppose you have a, a AAA with a gas cap in it and |
|
20:57 | the top of the uh uh reserve is gas and uh then um in |
|
21:02 | lower part is uh brine or maybe . So, uh uh so |
|
21:08 | in that case, the uh litho is plausibly the same above and |
|
21:14 | but that gas cap, the thickness the gas cap was simply established by |
|
21:19 | much gas got into the formation, ? Who says that the formation should |
|
21:23 | be filled up with gap? So it's only half filled up a |
|
21:27 | In that case, it's gonna be flat interface running right across the middle |
|
21:31 | or somewhere in the middle of the and it's gonna be flat because gravity |
|
21:36 | gonna make it flat. And uh of the uh uh the liquid below |
|
21:42 | interface is gonna be uh denser than gas. So that's an easy example |
|
21:47 | uh to imagine where uh uh uh real parts are uh small. But |
|
21:55 | so the imaginary part is not We could say that, oh I |
|
22:00 | , this is um this is for shale. But you could imagine also |
|
22:04 | a gas for the uh uh the on one side has a Q or |
|
22:09 | that's an ordinary Q and the Q uh on the other side uh is |
|
22:14 | low because of the gap. Then put those numbers into here and you |
|
22:19 | a pretty big um uh reflection of . It's a, it's a 4.5% |
|
22:26 | so 4.5% is a lot less than . But, you know, normally |
|
22:30 | deal with reflection coefficients with which are lot less than one. So normally |
|
22:36 | a strong reflection would be a number 10% or, or, or |
|
22:41 | Uh Those numbers are uh typical and here, we have 4.5% in the |
|
22:48 | uh order of magnitude. But it's this eye here, which means that |
|
22:52 | reflected wavelength comes back up is gonna phase shifted. And you can see |
|
22:59 | in your data, if you see uh uh reflections from uh interfaces above |
|
23:05 | horizon have uh uh are all uh pretty much the same shape of the |
|
23:10 | . And then you suddenly see at , at this uh target horizon, |
|
23:15 | see one that looks very different, you say ahuh that could be due |
|
23:20 | this kind of effect, we don't assume this is zero. Just assume |
|
23:25 | uh less than this. If this is comparable to this, then we're |
|
23:29 | get a, a phase shift. example, let me just if |
|
23:32 | if, if this imaginary part is the same magnitude as the real |
|
23:37 | then you get a 45 degree phase . And so that makes a different |
|
23:42 | coming out than coming in coming And furthermore, different from other nearby |
|
23:49 | which don't have this large jumping tube . So, um oh to |
|
23:55 | that's always been a very interesting possibility a new way to find gas. |
|
24:02 | , and the reason it's so interesting we as, as a result of |
|
24:09 | , we can call it AAA cue . We didn't lose any high frequencies |
|
24:17 | the reflected wave never got down into lower medium with high um uh high |
|
24:25 | , I mean, the lower meeting low Q. So it did, |
|
24:28 | didn't have any extraordinary attenuation, but simply change, change its shape because |
|
24:36 | this effect. OK. So let now go to uh OK. So |
|
24:50 | is Carlos. So uh so he says uh is QP and Q |
|
24:57 | incorporated in ci processing to compensate for attenuation. Uh So, uh I |
|
25:04 | to tell you that the uh the to that is not very satisfactory. |
|
25:09 | uh And the reason is because it's so easy to estimate what values to |
|
25:13 | , to put in. So I would have to say that that |
|
25:17 | current processing uses a lot of informal to uh deal with uh uh uh |
|
25:26 | attenuation. Uh We're not doing and nobody argues that we're doing it |
|
25:33 | . But what people argue, maybe uh uh it's good enough. |
|
25:39 | um obviously, the most important thing do is to uh uh uh we |
|
25:49 | gain correction. So you're accustomed to at uh your workstation and seeing maybe |
|
25:56 | second or maybe several seconds of reflection there. And they're all uh um |
|
26:02 | or less the same. But you that those long reflection times correspond to |
|
26:08 | uh long distances of travel and those when they were coming into the |
|
26:14 | they must have been a lot weaker the uh waves uh uh at the |
|
26:20 | of your work station screen. And we correct for that by what we |
|
26:24 | it, automatic gain correction. And we just uh uh uh we just |
|
26:29 | boost up those long arrival times without the uh the frequency content at |
|
26:36 | We just boost the amplitude, boost amplitude. Of course, we boost |
|
26:41 | noise as well. As we boost signal. So we have to be |
|
26:44 | little bit careful how we do But uh uh uh that's normally done |
|
26:50 | in the process. So you don't who did that. Uh What, |
|
26:54 | they were thinking was before it ever to you. There's been some automatic |
|
26:58 | control using uh you know, details that um of that calculation which you |
|
27:05 | not know about, but that's why can't ever uh uh believe and take |
|
27:12 | the amplitudes at different uh well, at, at different times in your |
|
27:21 | , which are uh separated a If they're only separated by half a |
|
27:26 | , maybe it's not a big But if they're separated by two or |
|
27:29 | seconds, then you know that the reflection got boosted up a lot and |
|
27:34 | doesn't show on your screen, you , that's so you can see it |
|
27:37 | that boosting hadn't been done, you not see it well with your |
|
27:42 | So that's how we do it. uh there's an example of correcting for |
|
27:47 | without correcting for frequency. And if know that that's true because uh uh |
|
27:56 | can see with your eyeball that the arriving reflections have lower frequency content. |
|
28:03 | whatever has been done with frequencies wasn't . Exactly. And of course, |
|
28:12 | reason is it's, it's hard to and uh and uh it's hard to |
|
28:18 | the cues. And then once you've the queue, there's a lot of |
|
28:22 | and so on. So you uh we don't tax free. So that's |
|
28:29 | remembering when you're doing a bo The question is from uh Rosa, she |
|
28:44 | in the attenuation lecture S slide 13 14, why is the sign of |
|
28:49 | imaginary part of m different for the elastic properties now? So let us |
|
28:55 | at that. I have that slide . OK. So let's just go |
|
29:05 | here. I forgot uh Rosada. , what uh number were you looking |
|
29:21 | this? The slide 20. Um me check in slides 13 and |
|
29:35 | That's what I'm gonna um go back . Uh Yes, I think it |
|
29:45 | here. This is the one Um mhm Now, what is your |
|
29:53 | ? Uh I uh I think it's the next slide. It says that |
|
29:58 | new and for ps that it's similar these are positive but for VP it's |
|
30:10 | . No sponsor. Oh So positive . And then it turned negative here |
|
30:18 | this is for the inverse, this for the slowness. Oh The |
|
30:27 | I wouldn't do that to you. these are not the Yeah. So |
|
30:34 | I had in the case of the , but in the case of |
|
30:37 | vs, it's the inverse of the here. And also in both cases |
|
30:42 | the inverse of Q and I thought become the became the inverse because we |
|
30:48 | the velocity as Well, you uh thank you for that. Uh |
|
30:55 | What I did to you, there um uh a typo uh you know |
|
31:01 | I, you know, obviously did I copied and pasted. And so |
|
31:05 | should be the real part of vs the positive power, positive one are |
|
31:11 | like this. So, uh uh you for that. You uh uh |
|
31:17 | , you found a typographical error Ok, thanks. Ok. |
|
31:24 | So I, I will uh correct and, and you should correct that |
|
31:28 | your notes. I will correct it re and repost that uh file |
|
31:37 | Ok. So that's good. So uh we finished. Oh, |
|
31:51 | and I also had another question. . Get a free question because you |
|
31:56 | so good on that one. You a free question, huh? |
|
32:02 | I send you the, the question well in the email. OK. |
|
32:07 | let, let me find that and it. Oh, yeah. Another |
|
32:23 | . The strength of the squirt mechanism attenuation depends on the fluid and the |
|
32:28 | b. What about the attenuation caused crustal heat? Hm, attenuation caused |
|
32:35 | the crustal heat? What are the values on top of the crust compared |
|
32:42 | kites and sediments or gas? I'm sure what you mean by attenuation caused |
|
32:48 | the crustal heat. Tell me about . Uh Well, I've seen in |
|
32:54 | images where the, where the crust very shallow, the sediments on top |
|
33:04 | attenuated compared to the ones that are . And there's a more, |
|
33:10 | more lithology underneath that they, where the cross is more buried. |
|
33:17 | . So, uh yeah. uh of course, rocks get hotter |
|
33:22 | we go down. Uh And that's you're talking about. The crustal heat |
|
33:28 | the way. So, here's a uh uh oh, here's a puzzle |
|
33:33 | you before I answer this question that ask. Let me pose to you |
|
33:37 | question. So we have the the earth is hot. Everybody knows |
|
33:41 | you got your volcanoes and everybody knows interior is hot and you know that |
|
33:47 | heat flows from hot to cold. the heat is flowing out of the |
|
33:51 | over the logic time um uh up the top and then uh radiated away |
|
33:57 | outer space. OK. So that's obvious. Now, how many of |
|
34:03 | have uh been out in the field gone into a cave whenever you go |
|
34:10 | a cave, the cave is right? I think everybody here has |
|
34:17 | inside of a cave, maybe a cave, maybe a small cave |
|
34:20 | whatever cave you go into, it's strong. That means that we have |
|
34:31 | full layer between the hot interior and outer space and, and the upper |
|
34:40 | . So why doesn't the um heat from the interior? Why doesn't it |
|
34:45 | up the um the cave? So cave is um uh almost as |
|
34:53 | uh, you know, almost as as, as the interior. How |
|
34:57 | we have a temperature minimum in the subsurface of the earth? Have you |
|
35:03 | thought about that? Uh, I think it's related to the heat |
|
35:09 | the sun, right? Because, , very good. So the sun |
|
35:14 | up the surface from above. and, uh, uh, |
|
35:19 | so then, uh, if you're the desert, uh, uh |
|
35:23 | and uh, the sun is beating , it's hot. But if you're |
|
35:27 | to a cave, you go in cave and it's, uh, cool |
|
35:32 | the heat from the sun has um, heated up the cave. |
|
35:37 | , uh, the heat from uh, the interior of the earth |
|
35:40 | heated up the cave either. It out that this is a trick question |
|
35:45 | , uh, it, it happens the, um, uh, the |
|
35:49 | inside the cave is equal to, , uh, the, the, |
|
35:55 | just a, a little bit warmer a little bit warmer than the |
|
36:01 | And, uh, uh, that, that contradicts your, |
|
36:06 | what I just said, you're standing the desert and getting heated from |
|
36:10 | but at night time, the sun not there and in the wintertime it's |
|
36:15 | . And so the, the, temperature in the cave is equal to |
|
36:19 | , uh, the, the yearly day and night of the, of |
|
36:25 | temperature at the top. And in , it's just a little bit warmer |
|
36:29 | that because the heat is continually flowing of the interior through the cave to |
|
36:35 | surface and then to outer space. there's no violation of the uh second |
|
36:42 | of thermos. OK. Now, to your question, what about the |
|
36:46 | caused by the crustal heat also, makes you think that uh the crustal |
|
36:52 | was causing that attenuation that you saw the slides that you saw, you |
|
37:00 | that it was uh uh some kind attenuation deeper than the shall or |
|
37:05 | But uh did uh uh did the who is uh presenting this material? |
|
37:11 | he say anything about what caused the ? No, no, I, |
|
37:17 | just saw that and I suspected that was because of the heat. |
|
37:22 | I think I've read in some, some places but not deep. I |
|
37:27 | , comparing the same section in portions the basement is shallower, very shallow |
|
37:35 | compared to where it's deeper, it's of the, the sediments just above |
|
37:43 | on top of that shallow cross or shallow basement are more attenuated. You |
|
37:51 | , I'm not familiar with this uh . Um Could you uh uh is |
|
37:59 | written down somewhere? Is this uh in, in a paper somewhere? |
|
38:05 | . Uh Yeah. Uh So it's for me to uh so, uh |
|
38:09 | know, um if the attenuation was , more generation means lower Q just |
|
38:18 | the basalt than elsewhere in the uh . Uh surely that would have attracted |
|
38:24 | attention of the expert and he would said something about that. You don't |
|
38:29 | what he said. No, nobody said anything about it. Yeah. |
|
38:35 | , if you come across that slide , I'd be happy to, |
|
38:38 | if you send it to me and me think about that. I |
|
38:42 | I can imagine that there might be , not for that, but it |
|
38:46 | be complicated because, you know, subsurface is complicated and I, I'd |
|
38:52 | to think about uh first, I think in my mind if you showed |
|
38:57 | that slide, I would think I would say, is that really |
|
39:00 | ? How did, how did he about this? What makes him think |
|
39:04 | the attenuation is high above the He better give a good explanation. |
|
39:09 | ought to tell me how he came that conclusion. And then once he |
|
39:14 | me that's uh uh really the then I'm prepared to think of. |
|
39:19 | how could this happen? And so just think here uh situation where the |
|
39:26 | sediments are normal except that they're Well, so maybe the fluids are |
|
39:35 | , you know, maybe um oh , I would have to be motivated |
|
39:42 | some data to, to think about uh what to think about explaining this |
|
39:51 | in a, in a complicated scenario that. So I'm gonna beg off |
|
39:55 | tell you uh send me uh uh slide that you're remembering. OK. |
|
40:02 | your next question is what are the values on top of the crust compared |
|
40:07 | Q values in sediments of gas? , that's a really good question. |
|
40:11 | , at the top of the let's uh think about the, the |
|
40:14 | 100 m. Sorry, it's mostly it's dirt. And uh we |
|
40:21 | uh figure is always almost always less consolidated than in the deeper segment |
|
40:29 | So it's gonna have a larger value and it's gonna have less cementation. |
|
40:35 | so uh uh uh furthermore, at uppermost part of the earth, it |
|
40:41 | be partially saturated with air, Why not? And so all of |
|
40:47 | things are going to increase the attenuation to the deeper cost. So we |
|
40:53 | be happy that uh uh we get seismic energy at all through the upper |
|
41:00 | because uh uh the upper, the worlds layers are really highly attenuated. |
|
41:05 | And uh luckily, uh our uh don't have to execute very many wavelengths |
|
41:12 | this highly attenuated, highly attenuated material um they would never make it to |
|
41:20 | receiver. Now, the sediment uh uh the rest of your question asked |
|
41:30 | two values in sediments or gas. , all those are also high, |
|
41:34 | high generation two numbers like five or less are you can expect in a |
|
41:44 | saturated sediment even if, even if well consolidated. Oh And I'll tell |
|
41:50 | this Yeah, this is a good to bring this up. Why do |
|
41:55 | explore with uh uh for oil and using c why don't we use some |
|
42:03 | thing? Why don't we use, example, a radar? So, |
|
42:09 | anybody have any ideas, why we explore for oil and gas? The |
|
42:15 | ? So Utah, you have an . Say it again. That's |
|
42:24 | So what Utah says is the seismic can travel deeper than the radar |
|
42:28 | And why is that? Because the wave has attenuated very high, |
|
42:34 | it uh it's got a Q value than one the radar. OK. |
|
42:41 | let me then uh pursue you on uh uh Utah radar uh has high |
|
42:48 | to it. And so why don't overcome this objection by doing low frequency |
|
42:55 | exploration? So have very low And I can tell you that the |
|
43:01 | who do this for a living, look at frequencies of about uh one |
|
43:06 | or less lot, a lot lower than radar. And so such low |
|
43:13 | do penetrate a lot deeper than the . Uh They can go a couple |
|
43:20 | kilometers down. And so uh uh yeah, the uh the attenuation is |
|
43:33 | . Uh uh The, the ainu by the time the signal goes down |
|
43:38 | uh uh a kilometer in back, has uh decreased in amplitude by over |
|
43:43 | million uh over a factor over you can barely detect it, but |
|
43:49 | have good detectors, they can detect . And so, um I can |
|
43:56 | you that that oh so electromagnetic wave inside the earth in a way is |
|
44:05 | similar to psychic waves with two One, the attenuation is very high |
|
44:12 | electromagnetic wave and it was very low of Q. And so they overcome |
|
44:17 | with uh uh low frequency. And and you know, from our uh |
|
44:23 | last time that since, since the is high, the dispersion is also |
|
44:31 | for electromagnetic waves. So different frequencies uh uh electromagnetic waves uh travel uh |
|
44:40 | different velocities depending uh on the And so that means that the uh |
|
44:47 | analysis of electromagnetic wave is very different the analysis for. Uh so for |
|
44:52 | sound wave, even though they're pretty the same, the only differ in |
|
44:57 | electromagnetic waves have higher generation and highest but if you account for that, |
|
45:04 | you can do a pretty good job , of e exploring with electromagnetic |
|
45:09 | And we had, I had a here uh a few years ago at |
|
45:13 | University of Houston who basically taught us to do that. And uh uh |
|
45:19 | ideas have been slow to uh penetrate the community of uh uh electromagnetic |
|
45:28 | But I think they're right. And think that uh uh eventually his ideas |
|
45:33 | become mainstream and uh and uh his will be famous. It hasn't happened |
|
45:42 | . So that's it. So, , the same way. But it's |
|
46:07 | uh um yeah, in, in same way. Yeah. Yup. |
|
46:16 | . So thanks for these questions, . And now let us go to |
|
46:30 | let us go to the anti oxy that dancing coldest we lost. |
|
47:04 | Uh So, uh uh I will uh uh sir, uh Utah was |
|
47:08 | happy with my, was, was satisfied with my previous explanation. And |
|
47:15 | goes further and he says in um waves, uh the the waves are |
|
47:21 | the particles that makes friction and the makes heat. And so that makes |
|
47:25 | heat. Well, so in the place, that's not a good uh |
|
47:29 | uh a good description of the way waves make heat. Uh When we |
|
47:36 | friction, it's uh uh uh sort meaning moving stuff like on a, |
|
47:42 | a table top that that's friction. uh So uh we don't have that |
|
47:52 | of friction in uh sound waves at amplitude, right? We, |
|
47:59 | we're not scraping atoms past other atoms all as uh in a sound wave |
|
48:06 | uh were uh sort of maybe uh uh putting a, an infinite testable |
|
48:12 | on the, on the grains and fluids in the rock. And uh |
|
48:18 | that is naturally making. Uh uh basically, I think the, the |
|
48:22 | to do the way to say it because of the second law. The |
|
48:26 | law of thermodynamics says that whenever you a change, uh uh it's gonna |
|
48:32 | increase the entropy uh meaning uh uh of uh of energy to heat. |
|
48:40 | so, when you think about it way, basically, it's the same |
|
48:43 | electro magnetic waves in electromagnetic waves, moving electrons around in uh inside the |
|
48:50 | the cloud of electrons surrounding each nucleus the rock. And uh uh |
|
48:59 | the second law of thermodynamics is going apply that we, we are going |
|
49:03 | uh convert some of that energy to the entropy. And oh no, |
|
49:13 | we assume that let's back up early in this course, we assumed that |
|
49:19 | sound propagation was adiabatic. That, , that none of the oh energy |
|
49:29 | lost the heat got outside of the . So we're here, we're not |
|
49:35 | about using heat outside the sample. uh the heat is still inside the |
|
49:41 | , but it's changing from uh mechanical to random uh vibration of the atoms |
|
49:48 | heat uh in, in every little tiny piece of the rock. And |
|
49:54 | think you can say the same thing electromagnetic uh uh uh exhortation uh uh |
|
50:02 | is, is uh deforming the cloud electrons around every nucleus in the |
|
50:10 | And that also is gonna lead to local increase of entropy according to the |
|
50:16 | law. OK. What's up? Now forgot where we left off. |
|
50:28 | we begin the talk about anisotropy last we met. So we didn't. |
|
50:35 | so we did a, a few . Uh uh So let me see |
|
51:01 | . OK. So, uh did see this slide before? OK. |
|
51:07 | , uh uh let me uh go and we saw these slides, I |
|
51:19 | it's maybe a good idea to uh repeat this. OK. So, |
|
51:27 | start here. And um so the slide says that uh Helmholtz, let |
|
51:37 | get it. Um uh Let me them. OK. So Helm's uh |
|
51:43 | uh theorem is still valid. It means that uh for these anisotropic |
|
51:48 | you can separate uh into a curl part and uh uh divergence free |
|
51:54 | But each of those is a mixture PNS. So you can't separate with |
|
51:59 | uh PNS uh uh in that But let's, let's uh um uh |
|
52:08 | the Helmholtz thing and uh uh look the scalar potential which uh is uh |
|
52:14 | curl free and it obeys this expression . So, uh uh uh uh |
|
52:25 | see what it's got here, it's lots of derivatives. See we got |
|
52:28 | derivatives with respect to the space and to space. And over here, |
|
52:33 | the left side, we have one with respect to uh space and two |
|
52:40 | respect to time. That's not the equation, isn't it? And, |
|
52:44 | we did use this expression and go uh um in your earlier notes, |
|
52:49 | see where we use this expression to uh a separate P and wave. |
|
52:57 | so we're gonna have to uh uh more clever. Yeah. So first |
|
53:03 | we do is we look at uh the stiffness a tensor. And remember |
|
53:10 | is not a tensor. This is matrix which represents the uh tensor, |
|
53:15 | has four indices. This has only indices. They count from 1 to |
|
53:20 | . And uh so here they And uh uh a as it turned |
|
53:25 | , the simplest case that we are with in geophysics ha has this |
|
53:31 | we call that polar an isotropic. it's for uh un fractured shale and |
|
53:38 | bedded sequences where all the dias directions the same, but they're different from |
|
53:43 | vertical direction. And so when you a stiffness element, a stiffness matrix |
|
53:49 | this, and we're only showing the triangle and lots of zeros and |
|
53:57 | Uh Yes. So uh this uh is sort of pointed in the wrong |
|
54:03 | . Uh uh This, th this , it should be pointing here and |
|
54:09 | point should be pointing here. I know exactly how uh it got messed |
|
54:15 | , but uh the 11 direction and 22 directions uh uh are the same |
|
54:21 | from the three directions. So we uh uh have this pattern. This |
|
54:30 | the same as this different from And down here in the shears, |
|
54:34 | have the same thing. This one different from this, the same as |
|
54:38 | and different from this, that makes different parameters. Here's 1/5 1. |
|
54:43 | then this one up here is calculated ones already developed five independent criminals. |
|
54:50 | one is gonna govern the vertebral P velocity. And you know that because |
|
54:57 | uh the, the first three uh we got two subscripts, |
|
55:01 | 1st subscript means that uh it's a unit area and the unit area is |
|
55:08 | in the three direction means it's uh the, the unit area is, |
|
55:13 | horizontal with its normal vector pointing in three directions. And the way was |
|
55:19 | going in three directions uh in the section that governs the uh P, |
|
55:25 | P wave vertical Vasi this governs the V these two go horizontal here. |
|
55:34 | got two different sheer moduli uh And one a a the main uh uh |
|
55:39 | also. And uh what we uh from this on the last time was |
|
55:46 | we can't uh separate the uh the modes and the sheer modes from each |
|
55:54 | using these two curl free parts and divers free parts uh uh Helm doesn't |
|
56:02 | us, heals is still correct, he doesn't help us. So what |
|
56:08 | do is we go back to this full tensor equation of motion. |
|
56:13 | we have this placement here. This almost like the wave equation, doesn't |
|
56:18 | , it's got uh two derivatives with to space two days, we respect |
|
56:24 | time of the displacement here and here our, our stiffness, a stiffness |
|
56:31 | with four indices. And uh we how to deal with this. |
|
56:37 | uh uh we've done this so many . We guess the answer and then |
|
56:40 | verify the guest. So we're gonna that the displacement uh vector is given |
|
56:48 | a con uh uh uh an amplitude in front. It's a vector. |
|
56:53 | you got a subscript I oh And see this is not very clever of |
|
57:00 | . I see this is not clever me because here I've used ISF as |
|
57:05 | hope effective component here. I used as a over to minus one. |
|
57:13 | this means that I there is another where I need to go back and |
|
57:20 | the uh materials so that you don't confused when you read about this |
|
57:26 | But for now you uh you'll recognize it, when we have a subscript |
|
57:30 | a vector component. And so the is a vector pointed in the uh |
|
57:36 | the I direction. And then uh have AAA re SOTO factor or phase |
|
57:45 | of either the IO mega plus or K dot X. So since we |
|
57:50 | K dot X, this can be way of going in any direction. |
|
57:54 | can we say about the wave vector ? Well, if it has a |
|
58:02 | length is given by this and if has a length uh uh related to |
|
58:07 | frequency, this gas is gonna solve equation. So that's what we, |
|
58:12 | put this into the equation of And uh now you see no more |
|
58:17 | because we've got make, we've executed derivatives and we got an omega |
|
58:25 | we got a minus omega squared on left side and we got a minus |
|
58:29 | on the right side and I'm not the minuses here. I don't need |
|
58:34 | anymore. OK? And you can that there's um a subscript I here |
|
58:41 | three paired subscripts JM and net. we, we're gonna sum with |
|
58:48 | And so that means there's gonna be term on the right side of this |
|
58:54 | when you spell out all the sums equations like that. So now we're |
|
59:01 | divide by the square of the length the wave vector both sides. And |
|
59:09 | uh this omega square to a case that is gonna be equal to the |
|
59:16 | velocity square, we know that and are gonna be what we call direction |
|
59:21 | that it gives the and the case of the wave factor in a ratio |
|
59:30 | the amplitude uh the magnitude of So that's AAA what we call a |
|
59:36 | cosine. And so when we put change of notation um uh into the |
|
59:44 | expression, we get this expression which in matrix notation looks like |
|
59:53 | And this is called an eigenvalue This is a famous equation in mathematical |
|
60:00 | . Turns turns up in lots of in mathematical physics. And what, |
|
60:05 | is this um operator L here? , here is the definition of uh |
|
60:11 | the operator here, you know, got that just from here. And |
|
60:16 | mm so uh the eye uh and uh the L, the L operator |
|
60:24 | two indices and that comes from or over uh uh J and N here |
|
60:33 | over J and N. And uh , more notation here, this identity |
|
60:40 | is just like the Chronicler delta it's a zeros off Dagle and one's |
|
60:47 | the Nagle. So this is a way to write that expression. And |
|
60:54 | we know how to solve this Physicists have learned, you know, |
|
61:00 | century ago how to solve an eigenvalue this equation. And uh uh immediately |
|
61:07 | and what are the unknowns here? The unknowns are the uh uh the |
|
61:13 | the three components of displacement. And immediately you can see uh a solution |
|
61:21 | this solution to this is U equal because it, if uh the displacement |
|
61:27 | zero automatically um salt, but we want a nontrivial solution. And |
|
61:36 | uh it's been known for, you , 100 years that uh you only |
|
61:41 | a nontrivial solution when the determinant of matrix is zero, not the L |
|
61:46 | , but the L minus V squared . So um do people know |
|
61:54 | what we mean by a determinant, . Do you, do you know |
|
62:01 | it? What you mean by a ? Yeah. So, um, |
|
62:07 | and I'm guessing that Carlos doesn't know . He, he's keeping a discreet |
|
62:13 | over there because Carlos is a geologist uh uh this is really, |
|
62:19 | this is really a, a AAA idea. And so, since it's |
|
62:25 | mathematical idea, let us uh go the support slide. So what I'm |
|
62:30 | do is stop sharing this and I'm to bring up. Um Thank |
|
62:41 | What I wanna do it. I think what I wanna do because |
|
62:53 | wanna go to uh the glossary. Well, I'm, I'm in the |
|
63:05 | file. You, you can't see but OK. Hm. Sorry about |
|
63:36 | looking for where I talk about the . OK. So I just, |
|
63:58 | uh I thought I just turn that . I know who's trying to call |
|
64:14 | there, but I don't want to to him here. So I'm gonna |
|
64:17 | close this and I'm gonna go to 101. Yes. OK. So |
|
65:22 | gonna show you this slide right Cool. Let's start off with this |
|
65:54 | in the spring. OK. Um you see this, this is um |
|
66:16 | M 60 from math one or And so uh we're, we're dropping |
|
66:25 | the middle of the discussion about tensors and here's an example of a |
|
66:31 | A three by three tensor. And it says is you can rotate, |
|
66:35 | can always find a principal coordinate system uh uh the uh stress tensor looks |
|
66:42 | this simple form like this. Remember all of these, these uh indices |
|
66:47 | , they refer to directions in the system. The coordinate system is something |
|
66:51 | decide on. And so maybe the doesn't know or care what you're |
|
66:56 | It has its own cord system, matter what it is, you can |
|
67:00 | find a cord system where the, tensor reduces to this diagonal form. |
|
67:07 | that's called the principal current system. these are called the principal values of |
|
67:13 | , the 10. Now, next , I said principal values, the |
|
67:21 | word is eigenvalues. And so we use that word eigenvalues even today, |
|
67:28 | those of us who are not And then it says the new directions |
|
67:33 | after we rotate from this uh into , the, the directions of the |
|
67:38 | system which leads to this diagonal form called the eigenvectors. OK. So |
|
67:48 | , what can we say um uh this uh in what we wanna do |
|
67:56 | say, what independent of the court ? What can you say? |
|
68:01 | uh uh look at this is the one says that if you sum the |
|
68:06 | eigenvalues, that's the same, that's uh uh that's invariant to the uh |
|
68:12 | court system. Let me back up . So the sum of these three |
|
68:17 | is exactly the same as the sum these three terms. So you don't |
|
68:21 | to find uh uh that magic court and then sum them up, you |
|
68:26 | sum these things up uh from the data because we know that uh um |
|
68:35 | , the, it's called the That's some of the, of the |
|
68:39 | terms. That's the same independent of cord system. And then here is |
|
68:45 | um quantity which is independent. But is still, I used, I'm |
|
68:50 | go directly to this third one. called the de the, the |
|
68:55 | And so what is this determined? the uh mhm trying to the three |
|
69:08 | back up here. So if you these three together, that uh is |
|
69:14 | number all the determinant. And uh uh so that's not the same as |
|
69:22 | these three together because these are not uh these are not the eigenvalues. |
|
69:28 | what is the de determine? Uh can you calculate the determinant from this |
|
69:36 | finding this first? OK. Without the principle value, the deterrent is |
|
69:42 | by this. So in two it's uh uh it's simple. It's |
|
69:47 | uh a 11 times a 22 minus product was mixed. These are |
|
69:53 | the off diagonal. If we had two di uh uh uh two by |
|
70:01 | matrix, then the determinant is, defined in this way in three |
|
70:07 | Uh It's more complicated. So let's uh make three column vectors, we |
|
70:14 | them a one vector, a two and a three vector. And then |
|
70:19 | determinant is, is defined by this product A one dot A two cross |
|
70:26 | three. So uh you can see uh number one, it's complicated and |
|
70:34 | um two, it's uh uh uh , the determine is gonna depend for |
|
70:40 | 3d uh tensor. It's gonna uh uh depend on uh it's gonna have |
|
70:46 | dimensions of uh of the cube, cube of, of the individual |
|
70:53 | And uh this is a CIC OK. So that's all I wanna |
|
71:00 | about the determinant. You can go to math 101 and uh uh uh |
|
71:06 | that again at your leisure. So now I'm gonna stop sharing and I'm |
|
71:12 | to oops. So c me up to um I'm sorry, what, |
|
71:53 | . This is where I'm gonna come , this is where we left |
|
71:57 | And so share that again with It's great. Very true. |
|
72:41 | OK. So this is where we off when we first introduced the term |
|
72:46 | determinant. So, uh when, we made the excursion back to math |
|
72:54 | , we found that the determinant is property of any matrix uh um uh |
|
73:00 | tensor. And we can take this the tensor inside this, this double |
|
73:08 | oh double bars here. And here means we're gonna make the determinant of |
|
73:15 | tensor and I gave a formula for , but it was complicated. And |
|
73:21 | what it says. In fact, here, it says it's a |
|
73:28 | it's a cubic equation in the elements this um uh of this 10, |
|
73:37 | right here. And since it's a equation, there's gonna be three solutions |
|
73:43 | uh uh uh the solutions, uh is the unknown or the unknown is |
|
73:48 | squared. It's not V itself, the square of V. So we |
|
73:52 | three different uh three different solutions for squared. And looking ahead, you |
|
73:57 | think now, OK, this is be V squared for PV squared for |
|
74:02 | and then another V squared for S be fast and slow values for S |
|
74:07 | , that's what we're gonna find, there yet. So those three solutions |
|
74:11 | called eigenvalues and the eigenvectors are gonna the polarization of those three ways. |
|
74:20 | that's what it says here. In case, we got velocity squared and |
|
74:25 | corresponding polarization directors. So uh let uh pause here for a moment and |
|
74:32 | well, thank you very much and a little quiz. So it says |
|
74:39 | um uh I'm gonna start with So pay attention here. It says |
|
74:45 | an isotropic elasticity is more complicated than elasticity because A B uh C or |
|
74:53 | and notice the last one is none the above. And so, um |
|
74:58 | let me ask you about um uh is a, a legitimate answer to |
|
75:05 | question. You think? Me? do you think DD for David? |
|
75:29 | ? Well, OK. So uh let's do this one at a |
|
75:33 | Uh uh But you, you think not a, yeah, we, |
|
75:36 | don't care about uh uh uh hetero or, or not uh on the |
|
75:42 | of a hand sample, we can uh anisotropic rocks, of course. |
|
75:48 | also you could have anisotropic crystals we it's on a slider though. So |
|
75:53 | not a matter of heterogeneity. Um And it's uh uh but uh do |
|
76:00 | think it's b um the C uh obviously not? Right. Uh I'm |
|
76:06 | go uh with, I'm gonna go B, I'm gonna uh uh agree |
|
76:10 | you on that. OK. Next , it comes to you Carlos. |
|
76:16 | need to use different mathematics for the case because ABC or none of the |
|
76:23 | . So, uh uh how about Carlos? Uh uh Do you agree |
|
76:27 | a Carlo? Yeah. Yeah, , I, yeah, I, |
|
76:44 | am not sure is valued for isotropic so that I wouldn't agree. I |
|
76:49 | that. Well, uh actually we about this a lot. Uh And |
|
76:54 | Helmholz was a mathematician and he didn't or care about um elasticity. Uh |
|
77:02 | isotropy is, but he didn't know of that. He just did his |
|
77:07 | and it's valid for uh any uh factor which varies as you know, |
|
77:15 | as function of position like displacement. uh And uh it doesn't uh uh |
|
77:22 | , he doesn't know or care whether rocks are isotropic or what. |
|
77:27 | uh this answer is not correct. Yeah. OK. So, um |
|
77:36 | , um turning to you, uh you got BC or D, |
|
77:47 | . 1st, 1st, let me you about C I is C |
|
77:52 | No, no, that's, that's not correct that the anti equation of |
|
77:57 | we showed you uh before uh it's , it is linear. It has |
|
78:01 | uh the unknown UO only appears to first B. So, so C |
|
78:07 | wrong um Now, but uh how B it says the scalar potential is |
|
78:13 | a solution to the equation of Yeah. Uh uh I'm gonna say |
|
78:36 | this is, I would say this a big of an ambiguous question. |
|
78:40 | , what do you mean by the of motion of? We did show |
|
78:44 | equation there as it had the uh potential in there. But was that |
|
78:48 | equation of motion? You know, think we meant we described that in |
|
78:53 | um uh in different ways. In past, I'm thinking this term is |
|
78:57 | bit, this one is a bit ambiguous. So I'm gonna put on |
|
79:02 | one And I'm going to uh uh you bra you don't have to answer |
|
79:08 | one because it's uh ambiguous. So get the next one which is, |
|
79:16 | this true or false. The so-called equation is just a special case of |
|
79:22 | simultaneous linear equations in three unknowns which homogeneous since all of the terms contain |
|
79:29 | of the unknowns. So there, , yeah, that's true. There's |
|
79:35 | in that equation like we had in uh uh uh I in the inhomogeneous |
|
79:42 | equation with the source term, there's source term there. Uh uh All |
|
79:46 | terms contain one of the unknown. that's true. Yeah, very |
|
79:50 | So now let's look at the So this is the uh the uh |
|
79:58 | stiffness matrix corresponding to the stiffness tensor the simplest case which we call polar |
|
80:05 | . The old fashioned name for this transverse isotropy. Uh I guess it |
|
80:10 | vertical, transverse isotropy because you can that it's this vertical direction which is |
|
80:17 | and the one direction and two direction equivalent. So what it says here |
|
80:22 | this uh uh uh cubic equation for determinant is too complicated and nobody here |
|
80:29 | how to solve a cubic equation. in instead of uh uh you |
|
80:33 | uh slogging our way through the uh uh let's try to be more |
|
80:38 | . And so um uh we know in isotropic seismic sh waves play a |
|
80:46 | role and, and they, they , are simple, they're decoupled from |
|
80:50 | P and SV waves. So, this as a clue, let's look |
|
80:56 | sh waves in this anisotropic case. . So uh we're gonna have, |
|
81:02 | gonna look for a solution which propagates the 13 plane, but it has |
|
81:07 | a cross line component U two So when we have these six, |
|
81:13 | this is gonna be now the, square velocity for the second Eigen value |
|
81:20 | here. And it's gonna be multiplying only the U two um component |
|
81:27 | of displacement because that's what it says . We're gonna be looking for solutions |
|
81:32 | only have a U two component. so, uh we're gonna call this |
|
81:38 | eigenvalue vsh square. And in terms the uh uh the uh right side |
|
81:46 | the equation, um all we have do is put into this position and |
|
81:52 | uh distance matrix two. And then have to perform all these songs. |
|
82:00 | ? And also we're gonna ha uh gonna insist that the wave is traveling |
|
82:04 | the K two direction. Uh That uh uh uh it's in the 13 |
|
82:10 | . So there's no uh wave vector the two direction. So that's gonna |
|
82:15 | things. And so, uh uh did we give this uh uh that |
|
82:20 | essence, you know, by OK. Now, putting uh we're |
|
82:26 | turn out to be a uh uh we do all these sums, it's |
|
82:30 | simplify a lot. So, uh , when we assume that K |
|
82:39 | we're gonna sim simplify the previous song to only these four terms. And |
|
82:57 | , here we go because of the and the stiffness sensor, some of |
|
83:02 | get uh um uh eliminated. So left with only a 2121 term and |
|
83:07 | 2323 term. So, for uh here we have a 23, |
|
83:12 | can't have M equals two. uh uh uh We can't have a |
|
83:16 | uh a 21 number here because that element is zero. And so, |
|
83:23 | of the zeros uh uh that's gonna it like itself. And so uh |
|
83:32 | changing to the two index notation, is gonna be a 66 and this |
|
83:37 | gonna be a 44. And will recognize these are um direction cosine? |
|
83:43 | so, uh uh the great thing it is because of these simplifications, |
|
83:49 | is only um on the right hand , we have only uh the twos |
|
83:56 | uh um uh you, you one and you three disappeared, we have |
|
84:01 | you two. And so we can by you two out. And so |
|
84:06 | get this solution directly divide out the two. And so immediately we find |
|
84:13 | the, the square of the uh , of the vsh velocity times the |
|
84:21 | is given by C 66 times sine plus C +44 times cosine square. |
|
84:28 | uh what are those angles, those are from the uh from the |
|
84:34 | OK. So, Um Let me at, at, at this figure |
|
84:40 | carefully, remember X two is out the plane in the plane is X |
|
84:45 | X three. So K is in plane and so this is the angle |
|
84:48 | right here. So uh uh and this is our Eigen Valley, we |
|
84:56 | very clever to um uh simplify the problem by recognizing that maybe sh is |
|
85:05 | play a special role. And sure it did. And what is the |
|
85:10 | for find this eigenvalue? Well, the displacement vector uh in the uh |
|
85:15 | the two directions. And we're gonna uh take it it, it's a |
|
85:21 | is one. OK. So that us the other two components coupled |
|
85:30 | Mm oh I see lots of um of um uh mistakes in the uh |
|
85:42 | of these slides. Obviously, this should all be on one line. |
|
85:47 | But uh when you look at these equations, it has the form of |
|
85:51 | two by two eigenvalue equation which is easier to solve than the three by |
|
85:57 | equation because the, the, the the determinant of this two by two |
|
86:03 | is uh so much easier. So it is the determinant of that two |
|
86:08 | two matrix is uh uh the product the, of the diagonal terms minus |
|
86:16 | product of the off diagonal terms. , you know, the, the |
|
86:22 | diagonal term is gonna have an L and A V squared times one that's |
|
86:30 | times the other diagonal term, which L 33 minus V squared times |
|
86:36 | That's this one. And then we to, to get the determinant, |
|
86:40 | have to track off the product of uh off diagonal terms. And so |
|
86:45 | is a diagonal uh excuse me, is a quadratic equation in the unknown |
|
86:50 | square which is um um uh a easier to solve than a cubic |
|
86:57 | And what it says here is uh second row of L has tensores |
|
87:04 | So the, the three showed up even though it's a two by two |
|
87:08 | because uh uh we eliminated um the right. OK. Now, the |
|
87:15 | to this are uh easy to work . Everybody knows how to solve a |
|
87:20 | equation. And so here are uh squared and VS squared uh uh |
|
87:25 | in terms of, you know, from uh the solution and it looks |
|
87:32 | easy, doesn't, it, it looks like it's some simple trigonometry |
|
87:36 | uh something called D or uh what that down here for? Sh look |
|
87:41 | this. Um It's the same simple here. And the only difference between |
|
87:49 | two is the minus instead A plus D is. Well, here's |
|
87:54 | wow. D is complicated. He complicated. It's got square roots, |
|
88:02 | got squares and fourth powers and uh a mess. And so these are |
|
88:10 | solutions for VP and VA SV. we still have a solution for |
|
88:16 | So uh that's simple. Now, want you to notice that each one |
|
88:21 | these expression and of course, we're be mostly interested in this one. |
|
88:28 | one of these contains four distance Let's count them 123. And the |
|
88:35 | one is inside D where is Here? It is out here. |
|
88:42 | , this should be a um a bit disturbing to you here is the |
|
88:48 | for VP. And you know that for vertical traveling VP, it's gonna |
|
88:54 | on C 33. you know how gonna work. But look right in |
|
88:58 | is AC 44, that's supposed to a sheer wave modules. But it's |
|
89:06 | here in the VP solution also down by the way. And also by |
|
89:12 | way inside here, let's take a here. So here it is inside |
|
89:19 | . OK. So let's draw now picture of these three solutions. That's |
|
89:23 | they look like. So here is uh the, the symmetry axis and |
|
89:28 | symmetry plane and those three waves all down at the same angle theta. |
|
89:35 | um the P wave here's the P traveling in the KP direction. |
|
89:40 | And here is the um um And can you see this he wave |
|
89:50 | not exactly in the longitudinal direction is exactly the displacement is not exactly pointed |
|
89:59 | the direction of propagation as you it is for isotropic meeting. So |
|
90:05 | can assure you that this little angle is not a, a graphical |
|
90:09 | That's um uh that's real and the of that little angle there uh depends |
|
90:17 | um uh the amount of anisotropy. course, uh I can tell you |
|
90:22 | nobody has ever uh figured out a to make any use out of |
|
90:27 | So we're gonna call the, we be calling these quasi P waves, |
|
90:31 | we're gonna call them P waves OK. So the next one is |
|
90:37 | SV wave, it's going down in same direction, it's following the KSV |
|
90:43 | , same direction here. Uh but different velocity. So uh uh it's |
|
90:49 | vector and this one is polarized in plane perpendicular to the direction, but |
|
90:57 | closely. Does that look to you it's not quite perpendicular? Yeah, |
|
91:02 | true in the same way. This not quite large Jual, this one |
|
91:07 | not quite perpendicular, but again, has ever figured out how to use |
|
91:13 | at all. And then the third is the sh wave. And you |
|
91:18 | see right here there is another one those uh font errors that uh mysteriously |
|
91:24 | up. That's supposed to be a The same size as this dot Just |
|
91:29 | dot Indicating that it's polarized out of plane. So I'm gonna have to |
|
91:34 | back and uh and repair this I don't know how it got me |
|
91:39 | up. Actually, I do know , uh, last, uh, |
|
91:43 | , a few months ago my hard , um, fail. It |
|
91:50 | no. And so I had a struggle to recover all my data |
|
91:57 | um, uh, download new powerpoint and everything. It's a big |
|
92:04 | Uh, no, the reason I you this is to, uh, |
|
92:08 | , warn you, you all should up your computer re religiously, |
|
92:13 | like every week on a schedule Friday before you go to bed, you |
|
92:18 | up your computer to an external uh, actually should, uh, |
|
92:23 | , back it up to an external so that if your hard drive fails |
|
92:29 | a Saturday morning or on Thursday uh, you won't have too much |
|
92:34 | to do to, uh, recover of your previous work. Especially those |
|
92:38 | you who have a lot of, , of, uh, uh, |
|
92:43 | on your computer which is gonna support supposed graduation from University of Houston. |
|
92:49 | , uh, I, if your drive were to fail you t you |
|
92:54 | be shit out of luck. So particularly, you back up your stuff |
|
93:00 | , with an external drive, you buy it for 50 for 50 bucks |
|
93:04 | , at the store, uh, , and back it up so that |
|
93:08 | your hard drive on your computer fails mark you're not lost. Otherwise I |
|
93:14 | you would be lost. You, lose a year of your life, |
|
93:18 | , uh, uh, you wouldn't able to graduate in time, |
|
93:22 | uh, you know, and, , maybe lose your scholarship apart. |
|
93:27 | would be terrible. You probably wanna your throat. Uh, uh, |
|
93:31 | do that drive and invest in in a, an external hard drive |
|
93:36 | back up your stuff every week or day. Probably not every hour, |
|
93:43 | a week is probably good enough. I learned that the, uh, |
|
93:47 | hard way, uh, you we, we get accustomed to technology |
|
93:52 | and when it works it's wonderful and it doesn't work, it's terrible and |
|
93:57 | broken heart cars are terrible. So was able to replace my hard drive |
|
94:02 | only about 50 bucks. It was cheap part. I lost all of |
|
94:08 | information of it. Terrible. So, uh, and the reason |
|
94:14 | I was struggling because I had been about backing up. So, whatever |
|
94:19 | I had were way out of date a lot of work to recover. |
|
94:24 | , don't let it happen here. . So, um, uh, |
|
94:30 | , your, uh, your turn says anisotropic rocks have ABC D or |
|
94:37 | of the above. Uh, or yeah. D all of the |
|
94:43 | . Ok. Uh Now, what's mean? It says at least five |
|
94:48 | constants. You got it. Why it say at least we showed five |
|
94:52 | , in, in the example we that was five. It has at |
|
94:58 | F what does that mean? here's the answer we considered only the |
|
95:06 | case of anisotropy. Only the simplest fractured shales well over anisotropy that has |
|
95:13 | . But uh the real earth is turn out to be more complicated than |
|
95:18 | . So it's gonna have more elastic system content. OK. Now, |
|
95:22 | you do, do have I, uh uh if it's a fractured |
|
95:27 | say with a single set of single set of vertical fractures in an |
|
95:32 | un fractured chair that makes it orthorhombic nine constants. Wow. And uh |
|
95:40 | now tell me this, look at no point, it says two |
|
95:44 | sheer body ways, not one, it doesn't say at least two sheer |
|
95:48 | . So if we got nine different constant, how many different waves are |
|
95:52 | gonna be having? Still only Because it's three dimensions, right? |
|
95:58 | the three comes from three dimensions. uh uh So for these orthorhombic |
|
96:06 | we're gonna have three wave types and one is gonna be complicated. |
|
96:13 | So um let us turn to the of weak polar anisotropy. So um |
|
96:32 | three people in history have understood the implications of the equations that I just |
|
96:38 | you five minutes ago. So the guy was a guy named Maurice |
|
96:44 | he was a professor in Polo and died 100 years ago. And, |
|
96:51 | he was the first uh person to call himself professor of geophysics. And |
|
96:58 | , um his research specialty was anisotropy which he understood polar anisotropy and he |
|
97:06 | everything. Uh and then he And, um the next guy who |
|
97:12 | it up was Klaus Helbig, who still alive, still with us. |
|
97:16 | is about 90 now, maybe over . So, um uh full of |
|
97:21 | , lives in Germany, mostly retired , but not completely. I hope |
|
97:25 | see him this summer. And uh he took it up and then |
|
97:31 | uh the uh the cause is smarter most of us. And so he |
|
97:36 | things in a way that nobody else understand only clout. And then the |
|
97:41 | guy who understood it was my Amaco , Joe Dellinger. You might know |
|
97:46 | name Joe is younger than me. uh uh joined uh Amaco when I |
|
97:52 | there. Still, he still works BP. He was the smartest guy |
|
97:56 | Amaco and now is the smartest guy VP. Uh And uh uh uh |
|
98:03 | those three guys are the only ones own, who understood those equations. |
|
98:10 | you might know that uh uh that am well known for anisotropy, but |
|
98:14 | not on this list. This is these three guys because they understood, |
|
98:20 | a lot more than I did. I found was that uh what I |
|
98:28 | was the following because these solutions which just showed you are so complicated. |
|
98:35 | gonna have to make uh some approximations for man, most of the last |
|
98:42 | , the popular approximation was to uh what so-called elliptical an isotopy. We're |
|
98:48 | talk, define that later. But not a good assumption because most rocks |
|
98:53 | not like that. A better approximation that the anisotropy is weak. Think |
|
99:00 | it. We have um we have found an awful lot of oil and |
|
99:06 | by assuming the anti start to be zero, assuming that isotropy is good |
|
99:13 | . And it was good enough to enormous amounts of oil and gas. |
|
99:16 | by the way, it was good to find uh all the features of |
|
99:20 | deep interior of the earth that we about, but maybe we can do |
|
99:27 | today. And so if you look at the previous expression, you can |
|
99:33 | the following combinations, two of them have the dependence of velocity, you |
|
99:39 | , was the the square root of modulus divided by density and uh three |
|
99:46 | dimensional parameters. You can see how three are non dimensional. I don't |
|
99:51 | , I, you know this one this one, obviously this one also |
|
99:56 | dimensional. And furthermore, you can that in the case where the anisotropy |
|
100:02 | zero, this one's gonna be a because if the anti be a |
|
100:08 | this element is the same as this . So this is a zero, |
|
100:12 | way this is a zero. And they think about it a little bit |
|
100:15 | see that one's also a zero in case of isotropy. So now let's |
|
100:23 | the case of weak anisotropy or we that these terms are not zero, |
|
100:30 | they're small compared to one. And we're gonna put that assumption into |
|
100:35 | um uh we're gonna rewrite the previous in terms of these things, you |
|
100:41 | there's five and all. So uh is as you look at it as |
|
100:45 | a rep parameterization of the exact So now we're gonna assume that these |
|
100:50 | small compared to one. And we're do a first order tailor expansion returning |
|
100:56 | terms which are linear in these three . And so when we do that |
|
101:03 | happens suddenly those three equations um are simple that anybody can understand. You |
|
101:12 | , for example, for the, P wave velocity, it's got |
|
101:16 | a reference velocity and too simple, isotropic, um too simple trigonometry, |
|
101:23 | isotropic um contributions. And these are be small, see we assumed uh |
|
101:32 | that uh delta is small and epsilon small. So uh uh these are |
|
101:36 | make uh only small changes here. how did uh uh how do we |
|
101:41 | them or the tailor of approximation? , before we go on, let |
|
101:47 | point out to you that this combination you see right in here, which |
|
101:52 | repeated right in here, that's gonna that a new name SMA. It's |
|
101:57 | uh not a new parameter. It's AAA new name. So now with |
|
102:03 | new notation, let us look at policy, these are the three equations |
|
102:08 | we just um uh looked at and apply to waves traveling at some uh |
|
102:15 | angle in the subsurface. Now, vertical propagation oops for vertical propagation, |
|
102:27 | we have sate equals zero. So term goes away and this term goes |
|
102:33 | and we're left only with the VP . So it wasn't that clever for |
|
102:37 | to call this reference velocity VP zero indicates uh uh propagating at the, |
|
102:43 | zero angle vertical. OK. look the neck at VSV. So |
|
102:50 | signed eight equals zero. So we're with VS zero just like it says |
|
102:55 | , FVSH against Syed equals zero for propagation. Again, we left with |
|
103:01 | it's the same sheer velocity for both modes for propagation. OK. So |
|
103:11 | us um look at horizontal propagation. there we have uh a sine theta |
|
103:20 | one, but course the equals So this term goes away again, |
|
103:26 | term becomes a one. And here have epsilon plus one times VP |
|
103:30 | That's the horizontal V velocity. Let's look at SV uh at, |
|
103:37 | , at the SV mode. Uh got um cos state equals zero. |
|
103:43 | this term is zero again. So left with only the VS zero for |
|
103:49 | it's the same as vertical. But course, in between, it's different |
|
103:52 | between is given by this. And on to the sh mode, we |
|
103:58 | scient eight equals one. So we gamma plus one times vs zero equals |
|
104:05 | . So uh you can see that waves traveling horizontally, uh shear waves |
|
104:11 | horizontally have different velocities whereas shear waves vertically have the same. Now, |
|
104:19 | of our data is P wave. let's concentrate on that. Here is |
|
104:24 | wavefront emanating from here and I'm showing two D cross section. And uh |
|
104:32 | so um it says that the wavefront a homogeneous layer is not a |
|
104:37 | So here's the wavefront and um uh is an isotropic circle to guide your |
|
104:45 | . And both of these have the vertical velocity. So this is a |
|
104:49 | taken after uh a certain number of . So the wave can go down |
|
104:54 | far, doesn't necessarily mean 2000 It just means it went down that |
|
105:00 | . You see the circle uh uh you have the 2000 on there. |
|
105:03 | circle comes out here at also horizontal , the wavefront is ahead. |
|
105:11 | So, uh uh so it's the VP zero here, but uh the |
|
105:17 | is not traveling along this line. is a, the isotopic circle. |
|
105:22 | is um oh this is the So the wavefront has gotten out ahead |
|
105:32 | horizontally. So it's not a circle it's not an ellipse. So just |
|
105:38 | show you that I'm gonna put on a perfect ellipse in green and this |
|
105:43 | an ellipse drawn um with the uh given to me by Mr Bill |
|
105:50 | And you can see it as um , it's tied to invertical and it's |
|
105:55 | at the horizontal, but in between doesn't tie. So we call that |
|
106:00 | epsilon ellipse. Why it's because it the uh el electricity is epsilon. |
|
106:09 | an ellipse has uh has only one of uh a degree of elliptic and |
|
106:15 | epsilon in this case. Now, so what we establish is the wavefront |
|
106:25 | not an ellipse. I'm gonna show now another ellipse which I call them |
|
106:31 | ellipse. That's this one in You see it's also matched here and |
|
106:36 | the original of the initial variation is to the wavefront. So you see |
|
106:45 | , the wavefront stays closer to the ellipse than it does to the green |
|
106:52 | . And then eventually at large enough it crosses over from the red ellipse |
|
106:57 | the green ellipse and ends up on green ellipse. So what is this |
|
107:03 | ellipse? It's an ellipse with elliptic . That's what it is. So |
|
107:09 | have the right to ask, uh . So um uh uh what does |
|
107:14 | mean here? Uh The, it's obviously not the same as the |
|
107:19 | velocity because that's what, what velocity this correspond to the real velocity is |
|
107:25 | here? Right? The, the wavefront has come and gone. |
|
107:29 | uh it's out here. And what this point here? Well, for |
|
107:34 | , it's simply the, uh the the notional idea that's where this d |
|
107:41 | comes out at the uh at the , the delta lips is defined. |
|
107:47 | it matches the wave front here, matches the initial uh variation, |
|
107:52 | an isotopic variation from the circle. the wavefront stays close to the delta |
|
108:00 | and then crosses over and ends up the absolut curve. So that's um |
|
108:05 | way the way P wave funds So now, before we, uh |
|
108:14 | we pass on, let me just uh ref refer this in two |
|
108:19 | One is gonna, I'm going to the and uh two dimension identities, |
|
108:27 | know, the cosine squared equals one sine squared. So just put that |
|
108:31 | here and we're left with this expression where, where the delta, |
|
108:35 | the minus delta sine fourth is now here. So we're gonna have, |
|
108:39 | have a name for this. Uh And so this arrow is supposed to |
|
108:44 | pointing here. That's the so-called an elliptic. A, an ecliptic |
|
108:50 | departures from elliptic. So if uh , if epsilon equal to delta, |
|
108:58 | would be a zero and we'd have these two terms and that would be |
|
109:02 | perfect ellipse. So the, a estimate uh uh uh measures the departures |
|
109:12 | electricity. OK. Now, another we can do is if we restrict |
|
109:19 | to small angles. The then we write this in this way only and |
|
109:25 | this term whether or not the coefficient a zero, if we restrict ourselves |
|
109:30 | small value. The because for small of data sine squared is small but |
|
109:36 | of the fourth is even small. we can uh simplify this for a |
|
109:41 | angle, you know, near vertical . So now let's ask ourselves, |
|
109:46 | does this show up in reflection move well. So here is the canonical |
|
109:52 | reflection problem. We got a homogeneous layer. You've seen pictures like this |
|
109:58 | , but those are isotropic pictures. so uh uh this is an an |
|
110:03 | layer. So here is the hyperbolic out equation and uh uh it has |
|
110:09 | Taylor series coefficient which is one over NMO velocity square. See the Taylor |
|
110:16 | efficient uh is, is the, derivative of the unknown here with respect |
|
110:22 | the small parameter. In this a small parameter X squared evaluated at |
|
110:29 | uh where the small parameter zero. that's the uh uh the Taylor series |
|
110:36 | . And so when you um uh this uh that you do this derivative |
|
110:43 | on the previous equations, you learn the move out velocity looks like |
|
110:48 | It is the vertical velocity times one delta. You recall that was the |
|
110:55 | mysterious uh uh uh velocity that we out before. That's this. |
|
111:06 | Now we're gonna go back for Now, here's a, here's a |
|
111:12 | good question for you. This equation evaluated near the origin for very small |
|
111:20 | , right? That's what it says here. So since this picture is |
|
111:25 | quite right, this uh uh this the wave vector should be coming down |
|
111:31 | almost vertically, right? Because we right here in the limit of small |
|
111:37 | . So the question is why isn't short spread move out velocity equal to |
|
111:42 | vertical velocity? You know, the the waves are going down like |
|
111:55 | I in the limit of small they should be traveling with uh |
|
112:01 | the, the a the vertical average . Well, this is a homogeneous |
|
112:05 | . It should be traveling with the vertical velocity, not the average |
|
112:17 | . This is a similar question to we had uh before uh uh uh |
|
112:22 | talk about uh from isotropic layers. said why is it the move out |
|
112:28 | A V zero vertical average instead of MS average? And we found out |
|
112:33 | that time, the reason is because what we measure is not this, |
|
112:38 | measure the horizontal move out that leads us the, the R MS |
|
112:44 | And so uh uh for the isotropic the anisotropic case, you're uh |
|
112:50 | homogeneous, no layers, homogeneous A isotropic case. Again, we're measuring |
|
112:56 | the horizontal move out right here. the uh uh as the delta in |
|
113:03 | and it's not just the vertical even though it's it's defined in terms |
|
113:08 | very short offsets, the figures uh in a bit uh defined in terms |
|
113:14 | very short offsets. But it it the delta in it right there is |
|
113:20 | delta in. So let's look at implications. This short spread move out |
|
113:29 | is hyperbolic even though the layer is isotropic. So when you see hyperbolic |
|
113:34 | out on your workstation, that does mean that the media is isotropic because |
|
113:42 | isotropic, meaning even if it's homogeneous leads to hyperbolic move out for short |
|
113:50 | set. So uh where is the I said, well, it's hidden |
|
113:57 | the move out velocity because what you on your workstation screen is this |
|
114:03 | You don't measure these two things you measure this so that uh the |
|
114:09 | ectopy is hidden inside what you And uh now it gets nasty. |
|
114:16 | you want to use the velocity, pick from your workstation to convert time |
|
114:22 | depth, you're gonna get the wrong because the depth is equal to the |
|
114:28 | time which you know times the vertical , it's not equal to the vertical |
|
114:35 | times the move out velocity. This is gonna be different than this one |
|
114:39 | of delta. So if you've ever a time to death this time, |
|
114:45 | only two explanations, you're screwed up or uh oh let's assume you didn't |
|
114:52 | that. The other assumption you might have made is you might have assumed |
|
114:56 | the media are isotropic. And we from the, from this argument that |
|
115:03 | media has uh has an isotropic move velocity even when we consider only very |
|
115:15 | offsets. Now here is a real right here. This one, it |
|
115:22 | the magni the an isotopy is magnified the move out because of the argument |
|
115:27 | . So follow me on this, the vertical velocity as a function of |
|
115:31 | angle uh offset. And let's consider only small angles. So we're |
|
115:37 | use the small angle approximation. And consider a case where uh the value |
|
115:43 | delta is 10% and consider a case all the angles are less than 30 |
|
115:50 | . OK. So uh the sine 30 degrees is one half square is |
|
115:56 | . So that means that all of rays in our gather are traveling with |
|
116:02 | um uh which differ from the um velocity by less than 2.5% which is |
|
116:10 | times 1/4. All the velocities are within 2.5% of the vertical velocity |
|
116:18 | So the move out velocity differs from vertical velocity by that full 10%. |
|
116:25 | see, because there's no signs for here science square time makes this a |
|
116:31 | uh uh uh uh small variations. even so the move out velocity differs |
|
116:39 | the vertical velocity by the full 10% them. So, you know, |
|
116:47 | out is our primary observable. The is a secondary observable primary observable. |
|
116:54 | one we use to make our images , that's the arrival times and the |
|
116:59 | of arrival time with offset which is out and in that move out the |
|
117:05 | shock we use in there right And it's magnifying in a magnified |
|
117:13 | So we derive this assumption, this using the assumption that the HM layer |
|
117:22 | uniform. Of course, that's not very good assumption. So let's uh |
|
117:25 | at a case where there's uh an uh uh layer one D me laterally |
|
117:35 | , but they're an anisotropic and then the move out velocity then is given |
|
117:41 | the R MS average of the vertical with an anti and I should copy |
|
117:49 | which I have a subscript on Uh It says RM SS uh a |
|
117:55 | on this average here. So that's in the lecture notes. Uh Not |
|
118:02 | election notes that you have. but uh in the, in the |
|
118:06 | notes which I published uh um for seg this is defined in detail. |
|
118:16 | I don't wanna say anything more about . But uh what I want to |
|
118:19 | is if you use, uh if say I want to know what are |
|
118:24 | interval losses. So if you find the interval VLO is either by migration |
|
118:32 | analysis or by diggs differentiation of of of velocities. Like we. So |
|
118:39 | before you find that you get for interval layer, you get an interval |
|
118:46 | which is contaminated by its own So for example, this subscript out |
|
118:52 | referred to the whole thing as the velocity for this third layer times one |
|
118:57 | delta for the third layer and the with all the other. So you |
|
119:01 | get away from the anisotropy by going velocity. You can't escape by computing |
|
119:11 | general velocity. And you can't consider escape it by considering only short offsets |
|
119:18 | in there. So you can measure by comparing VNMO with VP zero, |
|
119:24 | you obtain from a vertical VSP. to do that, you have to |
|
119:29 | drill a well, you have to 100 $100 million to drill a |
|
119:35 | the different, you always find differences which uh vary by the way uh |
|
119:40 | for delta as a function of de And uh because delta is right with |
|
119:50 | . So I think I I is possible for us to make this correction |
|
119:55 | growing? Well, that's all that's real bummer. Spent $100 million to |
|
120:00 | find the value of delta. Let's at longer offsets and observe the non |
|
120:05 | move out and use that to estimate Z. So I think I showed |
|
120:10 | this um uh picture before. Uh here we have a real data and |
|
120:16 | it's a, a common midpoint E gather flatly geometry. And a velocity |
|
120:23 | has been chosen to flatten the gathers short offsets, but at long |
|
120:28 | it doesn't uh work. So the thing to do is uh let's do |
|
120:36 | a higher order tailor expansion. This , have another term in tailor expansion |
|
120:41 | like we considered for the many layer case. And we'll make a physically |
|
120:47 | correction factor put in here uh so we're gonna choose this value a in |
|
120:53 | clever way. So that um at offsets uh this square and offset cancels |
|
121:02 | two of these. So again, end up uh T squared varying as |
|
121:06 | squared with the right velocity instead of wrong philosophy. Well, the way |
|
121:10 | do that is to define a in way uh uh uh but that requires |
|
121:16 | you need to know the horizontal So um when uh uh I, |
|
121:28 | was a co-author on the paper that uh pointed this out and I was |
|
121:33 | happy at the time, very pleased myself at the time, but I |
|
121:37 | realized that nobody else was particularly impressed that work. And the reason is |
|
121:42 | it required for processors to um determine each vertical arrival time. A, |
|
121:50 | short spirit move out of velocity, quantity, a four and a quantity |
|
121:56 | in order to flatten the gather. normally most um data does not allow |
|
122:02 | enough flexibility to determine three parameters at uh time. So a few years |
|
122:14 | , um uh somebody else made a invention and, and he found out |
|
122:19 | for a single polar ANAs tropic homogeneous polar, polar ANAs tropic |
|
122:25 | the previous formula simplifies to this. it looked the same as um I |
|
122:31 | showed you except that look here, got one, we renamed that uh |
|
122:36 | four to be a minus two A and look the same to a and |
|
122:42 | here. So for this formula, only have to uh uh determine at |
|
122:48 | uh uh vertical travel time it has travel uh to determine a short |
|
122:55 | move out velocity and an a, parameter and every uh there's two instead |
|
123:00 | three. So that suddenly makes it . And so what is this |
|
123:05 | in terms of quantities we already Well, it's epsilon minus delta over |
|
123:11 | plus two T one plus two delta um uh I should uh tell you |
|
123:22 | this guy uh Falcon uh was a at Amaco and then he went to |
|
123:27 | Colorado School of mines still there. his first student was this guy and |
|
123:33 | uh uh probably his best student was guy. And so al Khalifa came |
|
123:39 | with this and uh he, he that uh uh only have to make |
|
123:44 | single approximation, which is a good and you convert the previous column which |
|
123:49 | impossible to this one, which is . And uh uh uh compared to |
|
123:55 | with what we did before, uh is uh uh in lecture four, |
|
123:59 | came up with a similar equation, a similar notation. Uh But all |
|
124:05 | stuff that uh you don't need any this stuff. Now, the animal |
|
124:11 | don't need any of this stuff. so, uh uh uh now, |
|
124:17 | yeah, well, I th th shows the comparison between is should topic |
|
124:22 | we did before and I should topic we doing. Now here, we |
|
124:25 | two different parameters here to determine what's third one here, here, there's |
|
124:29 | two to determine that's really good OK. Now, that uh simplification |
|
124:41 | strictly valid only for um a sing single layer. But people use it |
|
124:46 | the time by putting in here an uh the parameter. And we determine |
|
124:52 | empirically by flattening the together for our of uh uh uh as a function |
|
125:01 | time. So when you do you, you uh transform the previous |
|
125:06 | uh problem. This, it's still perfect at, at f, at |
|
125:11 | offsets, we still have problems. before this major reflector was departing from |
|
125:18 | right about here. And now it flat for a lot further. |
|
125:26 | So I think this is a good time to stop for a break. |
|
125:33 | uh resume at 330 Kirsten type. uh uh I will see you |
|
125:41 | Uh Yes. Uh So, um Carlos, are you there? |
|
125:52 | Professor Rahi and uh Bruce, are there? It's not back yet. |
|
126:01 | Well, OK. So we're gonna anyway. Uh So this is where |
|
126:07 | left off. And uh what I to show you next is the effect |
|
126:12 | anisotropy on images. So, let show you what we have here. |
|
126:23 | We have some uh 2.5 D You know what I mean by 2.5 |
|
126:27 | modeling in the computer means that the the um the wave propagation in the |
|
126:33 | is three dimensions, but the model only two dimensions. So this model |
|
126:38 | into and out of the a figure variation. So that's called 2.5 D |
|
126:44 | . And so here is a, representation of the model in four different |
|
126:50 | . And you can see uh the shades of gray, we have the |
|
126:54 | the p the vertical P wave velocity in the upper left corner, you |
|
126:59 | see sedimentary layers here and you can see uh a salt body. And |
|
127:05 | so the uh uh the, the point of this uh study was |
|
127:13 | understand the effects of uh anisotropy in imaging. So now let's see what |
|
127:22 | we have uh we have here uh the upper right corner, we have |
|
127:27 | uh the representation of the model parameters delta. So you see insult of |
|
127:34 | zero. And then we have uh uh uh uh measures of uh uh |
|
127:39 | various values of delta in the very layers. Over here, we have |
|
127:46 | the uh the same parameters uh the the, the model representation for the |
|
127:53 | epsilon also in layers also zero in salt body. And here we have |
|
127:59 | a and you know that ADA is calculated from delta in epsilon. So |
|
128:04 | want to show you now uh uh these uh people on at all uh |
|
128:10 | for forward modeling, I think it Kirkoff Migra for forward modeling through this |
|
128:16 | 2 2.5 year model. And then Kirkoff migration of the results trying to |
|
128:25 | the uh uh an image of the model. Now, I think the |
|
128:29 | would have been the same if they'd a more modern form of migration. |
|
128:34 | That is to say it would be same an isotropic effects. Excuse |
|
128:42 | Can you repeat again, what is meaning of 1.5 a model? I |
|
128:49 | understand that part. What is uh mo why these are 2.5 models? |
|
128:56 | . Uh So the definition of 2.5 modeling 2.5 D is when you have |
|
129:01 | two D model, but you do wave propagation in there. So you |
|
129:06 | a source right here and the energy out in three dimensions. OK. |
|
129:10 | model is only two dimensions, so can call that 2.5 D. |
|
129:15 | Thank you. OK. That's uh that's kind of cute. But uh |
|
129:19 | I think that's a, a well uh notation. OK. Now I'm |
|
129:24 | show you the results of uh of in this um for this model with |
|
129:31 | different styles of migration case, I 44 different cells. So the first |
|
129:36 | here is uh uh an isotopic Uh uh So, and where they |
|
129:43 | they know what are the, the values of the parameters because they made |
|
129:48 | model that they know what they put . And so it's a pretty good |
|
129:51 | uh image and it's only for part the model, I'm gonna back |
|
129:56 | And uh uh it's this part of model right here, this part of |
|
130:01 | model right here. That's what we're on here. And I want to |
|
130:06 | out to you uh where this um uh this intersection between the fault and |
|
130:13 | betting is happening here. So, normally you cannot do this kind of |
|
130:19 | because you don't know what are the parameters can only do that uh uh |
|
130:26 | a uh kind of modeling environment. then they model in a different |
|
130:31 | then they model with an an isotropic and they fed the uh um uh |
|
130:39 | , they fitted the migration parameters uh with the uh move out velocity uh |
|
130:45 | you know, the flatten the gathers and find the best fit velocity for |
|
130:51 | the Gers. And then from flattening far offsets, they determined ADA and |
|
130:56 | they couldn't, they did not know is um uh the vertical velocity because |
|
131:02 | is uh fitted from the surface seismic . So they had to assume |
|
131:07 | So they assume that delta equals And so if you look at this |
|
131:12 | , uh it looks like a pretty image except that, and you see |
|
131:17 | intersection between the uh the uh the bed and the fault is happening at |
|
131:24 | lower level. This is a perfect line here. So because they uh |
|
131:30 | the model delta was not equal to , but here they were forced to |
|
131:34 | it's zero. So they got the bits. That's what I talked about |
|
131:39 | with you before you're gonna get the depths if you have the wrong |
|
131:45 | So the next time they uh next they did, they did isotropic migration |
|
131:50 | where they found the best fit but they assumed it's all isotropic |
|
131:56 | And you can see that this image uh distinctly inferior to this one. |
|
132:02 | It's got, got, you uh uh artifacts here and there. |
|
132:07 | uh uh uh so this is what state of the art at that time |
|
132:13 | the year 2000, I would say this is close to the state of |
|
132:18 | art today. Although they might make different assumption, you know, |
|
132:22 | if you're gonna assume delta, you have to assume delta equals zero, |
|
132:26 | could assume delta equals 5%. Why ? Maybe you have a nearby, |
|
132:32 | , uh that where they have a VSP in that. Well, and |
|
132:35 | actually measure uh uh delta. And , uh maybe from there, uh |
|
132:40 | you measure delta equals 5%. Uh why not just assume uh the same |
|
132:46 | uh nearby, probably wrong, but better than assuming delta equals zero. |
|
132:55 | in the fourth instance, what they was they drilled a well right here |
|
133:00 | they uh measured uh vertical velocities in well. And so it's a terrible |
|
133:06 | . You looked at this image and see all kinds of, of uh |
|
133:11 | uh uh uh artifacts in there because vertical uh uh well, it vertical |
|
133:17 | of uh velocity uh is uh is uh correctly known now because of the |
|
133:25 | uh from the VSB. But it's wrong velocity for imaging for imaging. |
|
133:31 | wanna use the NMO velocity. And that's why this image is wrong. |
|
133:37 | bad, but you see it's at proper depth, you see it's the |
|
133:42 | depth, whereas this one is deeper this one is deeper, this one |
|
133:46 | drilled in the right depth. It's bad image, but at least it's |
|
133:49 | right depth isn't that interesting. um uh these four panels give an |
|
133:57 | nice uh uh uh demonstration of the of uh including or uh ignoring anisotropy |
|
134:06 | imaging. And so, um uh days uh uh normally, what we |
|
134:12 | the best we can do is the of this is where we have the |
|
134:17 | fit parameters which we fit from the . But we know that's not enough |
|
134:21 | it doesn't give us doubt. So anticipate time to death. Misty, |
|
134:28 | you see here. So let's uh have a quick uh survey, uh |
|
134:35 | quick uh quiz here. Um For see, I think uh Carlos, |
|
134:41 | think it's your turn uh P velocity weakly or anti formations depends upon how |
|
134:50 | elastic parameters. All these choices. do you say, read it |
|
134:57 | We're talking about P wave velocity and an isotropic uh formations or assume it's |
|
135:05 | an isotopic. It depends upon how elastic parameters. No, Professor, |
|
135:13 | not sure. I be guessing. me hm. Well, I would |
|
135:20 | I would say three. But then , that's, that's exactly right. |
|
135:23 | . You saw the formula, it's zero, delta and epsilon. So |
|
135:28 | is exactly right. OK. So Brice, are you in, are |
|
135:31 | here now? Yes, yes, am. OK. So this one's |
|
135:36 | . Um uh It says here as statement, anisotropy effects move out velocity |
|
135:43 | uh ABC or D notice the D only about So, for, for |
|
135:49 | , I'm gonna do only a uh this is true no matter how short |
|
135:53 | the maximum of thats and the determination uh velocity. Is that true? |
|
136:00 | true. That's right. OK. uh uh le le uh how about |
|
136:05 | b also a trip? Um And Carlos uh uh c we got two |
|
136:16 | true. So uh uh either this , yeah, this better be |
|
136:20 | And in that case, we got of the above. OK. So |
|
136:25 | back to you uh Merce um says affects non hyperbolic move out. Now |
|
136:32 | this one is about a hyperbolic move . Uh So this is non hyper |
|
136:38 | of these statements describe this anisotropic Uh You see there's no, none |
|
136:43 | the above. So um uh it here, the non hyperbolic term depends |
|
136:49 | the near offset anisotropic parameter delta. two from sorry, I was |
|
137:05 | I think it is false. It be CC yeah. A better answer |
|
137:12 | C it de it depends on delta not delta alone, depends on |
|
137:17 | but not epsilon alone. It really upon a and has nothing to do |
|
137:21 | SV. So that's very good. . So my professor, I have |
|
137:27 | question if you can. Do you back to the example of that the |
|
137:32 | images uh migrated with different? So in four, the, the |
|
137:43 | in the isotropic migration and the isotropic , the velocities are the same. |
|
137:48 | the, but the anisotropy parameters are ones that are missing. That's |
|
137:55 | Here, we assumed the, anti zero, but here we found |
|
138:00 | best fit ada from the, non acro bolic move out. But |
|
138:06 | means that also the, the VN changes. Right. Well, the |
|
138:13 | is gonna be, uh, it's gonna depend upon, uh, |
|
138:16 | zero and delta in a certain combination you can measure this one from the |
|
138:22 | offsets uh independently of your measurement of one from the further offset. See |
|
138:28 | that works. OK. Yes. . No, no, I unders |
|
138:32 | understand now, thank you. now this is a very important topic |
|
138:38 | . Uh poor anisotropic A VL. so, uh uh I should tell |
|
138:43 | that I uh I should remind you I am the inventor inside chemical of |
|
138:49 | L back in 1981 I think. uh uh uh shortly after uh uh |
|
138:57 | and the way we discovered that is , we learned from uh a mistake |
|
139:02 | by uh uh uh uh uh a partner that they revealed to us |
|
139:08 | it could be done that a vo be done. And then uh |
|
139:12 | the uh well, hm the assignment figure out how to do it came |
|
139:18 | me, I was new in in-house . And so I quickly uh figured |
|
139:23 | how to do it and then we realized that that could be a great |
|
139:27 | way to reduce risk in drilling. uh we scrambled a tiger team of |
|
139:34 | experts uh went off site, looked a lot of data and confirmed that |
|
139:39 | vo you know, as you currently it, a vo uh we already |
|
139:44 | it basically in the same terms back 1981. And it can be a |
|
139:49 | way to reduce risk and prove that looking at a bunch of historical |
|
139:54 | So as soon as the importance of vo became obvious to everybody, they |
|
139:59 | the project out of my hands because was new hire and they were gonna |
|
140:04 | this over to one of those people have more confidence in. And so |
|
140:09 | what they did. Uh But I um um I continue to think about |
|
140:16 | unauthorized and I thought here's what I to myself. We're, we're doing |
|
140:20 | analysis. Um amplitude offset, amplitude with offset, it means amplitude variation |
|
140:30 | incident bank. But we're assuming that we're doing this, we are assuming |
|
140:35 | the rocks themselves are isotropic, shouldn't be considering a vo in the presence |
|
140:43 | velocity variation with OIE an iso? so I, I uh asked myself |
|
140:51 | question and I did a little, followed a little bit and decided the |
|
140:57 | was probably so, but I couldn't out how to do it. I |
|
141:01 | not figure out how to do a excuse me, could not have figure |
|
141:06 | how to do an isotropic a So um uh and now we |
|
141:11 | so I'm gonna uh uh talk about next. So here is the uh |
|
141:17 | that we had from before. Uh recognize the oh yeah, a vo |
|
141:26 | intercept great and character in terms of isotropic ideas jumps across the horizon of |
|
141:35 | is. So I think you're all with that, we analyzed it in |
|
141:40 | of this diagram here. And we the fact of uh uh by looking |
|
141:47 | the cross plot and seeing where these board um uh events, this is |
|
141:53 | uh intercept versus slope in the cross . And uh we saw these uh |
|
141:58 | things that look like uh on the on the cross pot, they look |
|
142:02 | like noise, but they all come the top of the structure. So |
|
142:05 | says they're not noise, they're telling a pattern and all that was very |
|
142:09 | good. And uh we, and by doing that, we learned how |
|
142:13 | use cross spots like this to find and to find fluid anomalies at your |
|
142:23 | greatly reduced risk. But there's something with this analysis. The slope of |
|
142:29 | green curve in here, which is taken from off to the side for |
|
142:34 | are no anomalous fluids in here and too steep. Something is wrong. |
|
142:40 | uh what I decide what I learned in 1981 is that if you do |
|
142:47 | same analysis with an anisotropic half but an an isotropic half space, |
|
142:52 | get a formula that looks pretty much same accept that in the gradient |
|
143:01 | And in the curvature team, you to find explicit and isotopic terms. |
|
143:07 | as soon as I saw it back 1980 I thought to myself, you |
|
143:12 | , this is big trouble because these here are gonna be small compared to |
|
143:18 | , but there's no one in the . All of these terms are |
|
143:23 | So this term which we always ignore be just as big as these other |
|
143:29 | terms could be, it could actually the sign to change the algebraic sign |
|
143:35 | the result if we include these So how to include that? Uh |
|
143:42 | , uh we can't solve uh uh problem like we did before because when |
|
143:47 | uh are trying to solve for these , which we, the vertical velocity |
|
143:52 | density, we got the uh the in anisotropy on the right hand |
|
143:58 | three equations in five unknown. So just to show you how crucial |
|
144:08 | is, we're gonna do an exercise now. So let's get out of |
|
144:14 | . So I'm going to stop And while I'm fumbling around here, |
|
144:20 | bring up the classroom exercises which I shared with you before via canvas. |
|
144:28 | what I'm gonna do right now is , I have to search around a |
|
144:37 | bit. Hold on a second. , I'm browsing for the file that |
|
144:58 | want and you should be browsing for , for the classroom exercises, trial |
|
145:10 | I provided to you for. it's called a VO exercise. |
|
145:20 | So, uh, now I am to show this, iii I don't |
|
145:27 | it. I cannot find it in en canvas. Oh, you don't |
|
145:32 | this. Uh, it might be I forgot to give it to |
|
145:36 | OK? I will check because this probably the most important part of this |
|
145:42 | course right here. And so uh uh um what I'm showing you here |
|
145:49 | an a vo exercise and uh I'm share with you, share the |
|
146:05 | I'm gonna have to stop sharing the screen. Can you see this um |
|
146:57 | the spreadsheet with the red and yellow and some curves at the bottom and |
|
147:03 | at the right? Yes. So let me show you what we |
|
147:08 | here and you can see my cursor around. OK. So never mind |
|
147:12 | tables to the right. They just the cursor. So what we're gonna |
|
147:16 | is some forward modeling of anisotropic A . And so we have uh uh |
|
147:21 | model parameters here in color. So got a red formation over a yellow |
|
147:27 | in the red formation. We got uh VP zero, which is gonna |
|
147:31 | you're gonna s uh set it to own parameters by uh using this |
|
147:36 | So you put your cursor on top the slider, depress the left key |
|
147:41 | then move it and you see the uh the number here changes. And |
|
147:46 | that's a different velocity that I OK. And the same thing with |
|
147:52 | ratio down here, same thing with , same thing with delta and |
|
147:57 | So what the first thing I'm gonna is I'm gonna put the uh the |
|
148:01 | delta equals zero here and I'm gonna epsilon equals zero. And I'm gonna |
|
148:09 | the same thing down here in the body. OK? Now you're gonna |
|
148:17 | down a little bit, scroll down little bit. And so the curves |
|
148:24 | here and that's a typical a VO that's exactly calculated from uh the model |
|
148:30 | I showed here. Um I shouldn't exact. So this is uh the |
|
148:35 | , the, the linearized approximation, A VO and this is what you |
|
148:39 | calculate for uh uh uh for a that you do. And uh so |
|
148:45 | can see that there are three values uh three curves plotted on top of |
|
148:50 | other. OK. So uh and I'm just gonna uh change uh |
|
148:55 | the, the velocity ratio up here you can see how the curve |
|
149:01 | OK? As soon as you let of the, of the slider, |
|
149:05 | the curve changes. OK. Now gonna put in here uh oh and |
|
149:09 | the way uh we got here. correct. Thank you. OK. |
|
149:33 | I got 0% here and this is um uh delta and epsilon. And |
|
149:39 | are the values which are uh reported this uh field from the, from |
|
149:45 | slider. OK. And down here the calculator value for A. |
|
149:51 | So now let's um iii I put here just oh And by the |
|
149:56 | these are, are um uh uh scale here is exactly given from |
|
150:05 | So, so it shows a 6% um uh P wave reflection coefficient. |
|
150:14 | are exact numbers. So let's put here a little bit of delta. |
|
150:18 | I'm gonna put in, so imagine is a shale above. So I'm |
|
150:22 | put in a little bit of delta . OK. So I put in |
|
150:31 | of delta. Suddenly we see different uh um different curves. And so |
|
150:39 | got this uh uh the uh the in dark blue, that's the same |
|
150:44 | as we had before. That's the curve. And here we see down |
|
150:49 | two curves on top of each One showing the uh uh the Avio |
|
150:55 | uh e exactly calculated properly. And other one calculated not quite right. |
|
151:02 | um Never mind the difference from And o off here to the |
|
151:08 | you can see the uh uh the gradient which is 2.9% for the isotropic |
|
151:16 | and 1% for the an isotropic. uh by putting in just a little |
|
151:23 | of delta right here, I changed gradient by a factor of three. |
|
151:30 | you see that changed the gradient by factor of three? But only changing |
|
151:36 | uh the delta but a little And if I, you know, |
|
151:43 | the epsilon, I'm gonna put in a little bit of epsilon make it |
|
151:51 | because usually epsilon is bigger than And now things are changing some |
|
151:57 | So that did not change the um gradient. If you look back at |
|
152:01 | formula, uh uh the gradient depends only on delta uh the jump in |
|
152:09 | in uh upper and lower. And I just now changed uh uh |
|
152:13 | And now you can see that these two anisotropic curves are a little |
|
152:18 | different. Uh So, uh since is short, I'm not gonna talk |
|
152:23 | these differences only. Uh but uh only about these differences only. And |
|
152:28 | main thing is that by um wi imagining just a little bit of |
|
152:36 | we changed the gradient by three times change in. And uh and, |
|
152:44 | , but we do this all the , we ignored Dal. And what |
|
152:49 | exercise proves to you that delta could very important depending on the other |
|
152:54 | So let me just uh change uh up here and I'm gonna uh change |
|
152:59 | of these down here. This is change the velocity ratio in the reflecting |
|
153:05 | . I'm just gonna move it And when I let go, uh |
|
153:11 | see that everything changed, of Uh um but the um it's what |
|
153:19 | see and you see, I didn't the anti I just changed the, |
|
153:23 | ratio of vertical velocities in one And again, there's major differences in |
|
153:29 | gradient and uh now they're both negative major changes in the gradient. Uh |
|
153:36 | uh uh that uh comes from uh uh all these terms work together and |
|
153:43 | relative size of the anti soy terms depends upon the uh the size of |
|
153:50 | isotropic terms. So I just changed isotropic term and you saw the |
|
153:56 | So this uh spreadsheet exercise. Uh very simple but it should uh give |
|
154:06 | uh grave concerns if you're interested in ll because you can see that the |
|
154:12 | that when we ignore an isotropy, could be making major major errors in |
|
154:22 | analysis because the anti soy term might just as big as the terms that |
|
154:28 | calculated. Yeah. So think about and talk about it with your other |
|
154:35 | in your company who are uh um vo expert. And think about |
|
154:39 | I will then I will upload this um spreadsheet to you tonight. |
|
154:46 | I should have done that before. And you should talk to your um |
|
154:51 | talk about this with your friends and uh guys, we could be making |
|
154:56 | big mistake by ignoring anti. So since time is short, I'm gonna |
|
155:03 | at that. And since I didn't uh give it to you, I |
|
155:06 | to spend some time uh allowing you play with this, but I will |
|
155:10 | it to you tonight and that you uh uh play with it and then |
|
155:18 | it to your colleagues and say what about this? And now I'm |
|
155:24 | go back to the, as a . Yeah. Are you with |
|
156:09 | OK. So this is where we at office. And so uh I |
|
156:15 | um discovered this effect back in 1981 before you all were born and uh |
|
156:26 | I could not figure out how to determine that parameter delta delta from the |
|
156:35 | . Easy to imagine what it could . And, and uh the uh |
|
156:41 | allows you to put in some plausible and immediately you're shocked by how large |
|
156:47 | contribution it makes, but I couldn't out how to evaluate it from the |
|
156:54 | . And so that's what I'm gonna you next. So uh uh he |
|
157:03 | our problem to uh analyze the A data while accounting for the anti and |
|
157:09 | issues are that logs have high spatial and accuracy. They don't measure |
|
157:16 | Why is that? Because they have array pads, only vertical ray pads |
|
157:25 | . Now the surface seismic travel times low special resolutions. So uh you |
|
157:31 | that um delta is uh in the out velocity and you can uh um |
|
157:38 | can deduce it by comparing with vertical from the VSP. But even if |
|
157:43 | did that, it's gonna have low resolution. So uh uh so you're |
|
157:50 | gonna be able to find a local in delta from surface seismic travel |
|
157:59 | You can find trends in delta, you can find uh uh you |
|
158:04 | uh uh of course, averaging of , but you can't find a local |
|
158:10 | in delta from surface travel times. when you look at surface amplitudes, |
|
158:17 | figure it's in there. But uh the amplitudes are affected by many |
|
158:22 | We talked about those uh of those you can imagine them uh before you |
|
158:27 | , there's uh uh uh geometrical there's a generation, there's uh uh |
|
158:35 | effect of trans uh transmission uh coefficients the overburden. There's uh uh uh |
|
158:43 | as with, there's a angle variation the source strength. So many things |
|
158:48 | you can think of uh affect the a and basically, it's simply not |
|
158:55 | to make deterministic corrections. So this was solved uh with a uh master's |
|
159:02 | here at the University of Houston um 10 years ago by now. And |
|
159:06 | this is the uh the um the you can uh if you're scribbling it |
|
159:12 | . Well, yeah, you have in in your files. But uh |
|
159:16 | you have to do is look up these names in the seg annual |
|
159:21 | Expanded abstract from 2013. And so is our, this is what we |
|
159:28 | thinking. Remember this convolutional description of propagation from less than four uh uh |
|
159:36 | our seismic data. Uh uh uh start right here. So here is |
|
159:41 | source strength uh uh parameter. Here's initial wavel, then it propagates down |
|
159:47 | some complicated uh operator that we don't what it is, then we have |
|
159:52 | reflection and then it propagates back upwards a complicated operation. And we have |
|
159:58 | summation over many different reflections. And the data then gets involved with an |
|
160:03 | response and uh something is going on the computer. And then on top |
|
160:08 | that, we have all this So uh most of these different effects |
|
160:13 | result in angle dependent variation of amplitude is offset the penetration that but most |
|
160:21 | these a vo effects accumulate gradually as go down a little bit and down |
|
160:29 | little bit more and down a little more. All of these things change |
|
160:34 | . And the only thing that changes is this one right here, their |
|
160:40 | . So when you go to f to Hampton Russell, they are going |
|
160:46 | uh compare the surface seismic data with uh uh reflect a synthetic reflectivity which |
|
160:54 | get from the law. And so surface agent data has units, you |
|
160:58 | , plus or minus 1000 or so . But the the the reflectivity are |
|
161:03 | probably less than plus or minus So what uh Hampson Russell is gonna |
|
161:09 | ? They're gonna normalize the surface seismic intercepting curvature to an isotropic synthetic reflectivity |
|
161:17 | which is constructed from logs and doesn't any of the propagation effects in there |
|
161:27 | these normalization. The slowly growing parts good acquisition and propagation effect. And |
|
161:33 | only part that varies rapidly with vertical time is the reflectivity. So what |
|
161:46 | Thompson did in 2013 was to realize what this uh these normalization factors, |
|
161:53 | functions of time, there are functions uh uh vertical travel time. And |
|
161:59 | so um uh that means they have , a fourier um spectrum and the |
|
162:07 | spectrum has uh high frequency parts and frequency parts. So if you don't |
|
162:12 | this, if you don't um uh e everything, all the surface uh |
|
162:18 | seismic through this uh synthetic gather, only um only use the slowly varying |
|
162:28 | , slowly running parts of the app factors use that only to normalize thereby |
|
162:38 | correct for the propagation acquisition effects leaving rapidly changing reflectivity unchanged. So here's |
|
162:46 | example of that, this is from uh master's thesis 11 years ago. |
|
162:53 | here we have a vertical arrival time milliseconds from 400 milliseconds to 700 mills |
|
163:00 | there, we had um uh 123456 reflectors. And so uh um we |
|
163:11 | uh and we have had a seismic uh at the midpoint right at that |
|
163:19 | right at that at the where the was held. And from the seismic |
|
163:24 | , we took all those seismic um uh received gradients and divided them by |
|
163:32 | factor of uh 10,000 so that we plot them on this graph. And |
|
163:37 | slopes looked like this for each of various uh reflectors. Understand how you |
|
163:44 | take a surface seismic ganner. Look each, you know, you flatten |
|
163:48 | and do everything. All the magic Hampson rustle uh gives stretch and squeeze |
|
163:53 | and everything. So it get to the well. Exactly. And then |
|
163:58 | look at the um um uh at uh the gradient for that seismic data |
|
164:09 | you normalize it to the isotropic reflectivity or for you, you normalize it |
|
164:16 | by 10,000. So you can put on this block and then uh look |
|
164:20 | the syn syn synthetic isotropic reflectivity slope the logs. So from the |
|
164:26 | you can calculate R two and R and R zero everything. And that's |
|
164:31 | red curve. And now if we this to this, we get |
|
164:37 | we learn nothing. All we and what Hampson muscle does. It, |
|
164:40 | teaches you how to model surface seismic terms of sy synthetic. Um And |
|
164:47 | terms of isotropic changes, we want find the local job in delta. |
|
164:53 | what we do is we low pass that normalization function and only normalize with |
|
165:00 | low pass fraction of the normalization So you see we did not get |
|
165:07 | normalization. Matter of fact, uh every one of these is wrong and |
|
165:11 | wrong by a little bit. You this one is wrong by 30%. |
|
165:16 | one is accurate. This one has opposite side. So these, we're |
|
165:24 | attribute these differences to uh uh the term of delta delta in the |
|
165:30 | which we could not measure in the . And we could not measure in |
|
165:34 | seismic. We're gonna estimate it by procedure and we just noticed that these |
|
165:39 | are not small as a um a percentages of the total. So now |
|
165:47 | do we know whether or not we're a reasonable answer? What are we |
|
165:52 | use for ground truth? Well, we're estimating here is delta delta. |
|
166:00 | no way for us to measure no matter what we do uh in |
|
166:06 | earth for delta delta. But what can measure is delta. And so |
|
166:13 | we're gonna do is we're going to that at the uh uh uh we're |
|
166:18 | assume a value for delta at one these layers and then add up all |
|
166:23 | delta deltas and get an absolute value delta itself. In each of the |
|
166:28 | , you know how that works. , and uh uh uh the, |
|
166:36 | , so that process is very Let me show you the next slide |
|
166:40 | . This is the deltas that we uh uh came up with. And |
|
166:46 | gonna start off here with delta equals right here. And you can see |
|
166:51 | the maximum here is about d equal 20% possible. But who knows whether |
|
166:59 | right or not, where did we up with this zero right here? |
|
167:03 | , you can see the gamma ray over here and uh uh uh right |
|
167:08 | is a uh as a player with gamma ray. So we uh gone |
|
167:14 | decades. What that means this is sandstone, but let's assume that in |
|
167:19 | sandstone layer, delta is zero. that's why we come up with this |
|
167:23 | here, then we add delta delta we add delta delta again. And |
|
167:27 | , and we come up with this . And so uh uh what you |
|
167:31 | here is that uh um uh that's calibration sandstone. And you see that |
|
167:38 | we've uh calculated uh uh high values delta, those correspond to high values |
|
167:45 | the gamma log, gamma ray So everybody knows that this means that |
|
167:49 | is here. And so that's what predicting uh with high values of |
|
167:54 | That means uh we're expecting high values delta in a shale. And sure |
|
167:59 | , we see a shale here. I si sent Miss Lynn into her |
|
168:05 | thesis uh examination with this. I it was pretty good. The other |
|
168:12 | here at the university who were much hard nosed than me. They said |
|
168:17 | ? You've only got two successful predictions you call that a success. I |
|
168:22 | , well, she said, you , she was answering their criticisms on |
|
168:26 | feet in front of the examination committee she said, well, you |
|
168:31 | I I only had this much log . So one of the professors here |
|
168:35 | named Professor Stewart, you probably know . And he said, well, |
|
168:40 | uh I know where this data set from. This comes from Canada, |
|
168:44 | from Canada. I know that We have lots more logs from that |
|
168:48 | field here in our files at the . You go into our database, |
|
168:53 | University of Houston, find some nearby from nearby wells in the same |
|
168:59 | same formation and add them in. that's what she did. And so |
|
169:03 | a nearby log, she was able extend the um uh the logging, |
|
169:09 | gamma ray logging uh higher up. you see, here's another sandstone and |
|
169:13 | shale she predicted low. And then said that's five out of five so |
|
169:18 | they gave her a pass. So uh uh uh graduated with distinction with |
|
169:23 | master's thesis in 2013. And I hoping she would come back for a |
|
169:31 | thesis but she said I'm tired of poor. I wanna go get a |
|
169:38 | . And so she went out and a job. Well, she made |
|
169:43 | money but she wasn't having a lot fun. So she, uh, |
|
169:47 | back to the university a few years and got a phd and now she |
|
169:51 | a good job somewhere else. I forgot, uh, where, |
|
169:55 | I'm sorry to report that she uh, do for her phd. |
|
170:00 | did not do an extension of this . On the other hand, what |
|
170:03 | means is an extension of this work still available for, uh, U |
|
170:09 | H students uh, today. And can think of lots of ways |
|
170:13 | uh, uh, do better. , this was, uh, |
|
170:20 | just a master's thesis and I think are several phd thesis, uh, |
|
170:26 | to be, uh, pursued along same lines. Well, the University |
|
170:33 | , uh, Houston has a betting this and we hope to find something |
|
170:37 | make some money out of that. I would have to say that so |
|
170:40 | . We haven't made any money. , so that is all, |
|
170:46 | uh, uh, to the point this study limited as it is shows |
|
170:53 | the, uh, the, uh, the plausible, um, |
|
171:03 | that the neglected effect on polar an coming from polar an iso could be |
|
171:09 | important. And this data set shows , in fact, in this |
|
171:13 | it, it was and, and is. And uh what I think |
|
171:18 | does is it, it uh goes doubt every single isotropic a analysis ever |
|
171:26 | . So every such isotropic Avio analysis in your company could be rethought and |
|
171:34 | it needs to be revised or maybe depending on uh what are local values |
|
171:39 | delta delta in that area. And this algorithm, you can figure it |
|
171:47 | and furthermore, you can probably um out improvements to this algorithm. So |
|
171:59 | leave that to your imagination. So do a little quiz here. Here's |
|
172:03 | quiz. Um because the anisotropic effect A VL was not recognized for so |
|
172:10 | , 30 years now. ABC or or all of the above it |
|
172:15 | Uh um So um uh li how about a, is that |
|
172:26 | Yeah, that's true. Uh uh B is this true when you look |
|
172:32 | seismic surface seismic amplitude, they contain variations, but those are due to |
|
172:39 | propagation and reflection. Is that I think it's true, professor. |
|
172:46 | true. And we just named it of them and you can think of |
|
172:49 | probably um uh uh versa. How you? This, the seismic altitudes |
|
172:55 | on you usually normalized using logs and normalization forces the agreement with isotropy. |
|
173:03 | can't learn anything about anisotropy. yeah. Uh uh uh uh uh |
|
173:10 | including I think, and uh uh sha software uh so that means I |
|
173:16 | read all of the above. So uh back to you le le |
|
173:21 | says uh anisotropy effects which coefficient um gradient curvature B and C are all |
|
173:30 | the above. What would you Mhm OK. And we talked about |
|
173:45 | , we didn't talk about C but uh we didn't talk about it much |
|
173:48 | we showed it. Yeah. So uh so it's B and C and |
|
173:52 | we didn't talk about C A curvature that's usually very noisy. Mhm |
|
173:59 | So um uh Carlos since the anisotropy usually weak when it's, you |
|
174:06 | measured as a rock property, uh effect on the AC O intercept is |
|
174:11 | small, true or false. I it's false professor. And why is |
|
174:18 | false? Because I mean, the can be uh yeah, you, |
|
174:22 | , you were actually showing that we have a significant significant impact in |
|
174:28 | in the response of the A what said is true, but you did |
|
174:32 | read carefully. See it says the Aviel intercept, you were thinking |
|
174:38 | the next quiz question, which is A vo gradient. So that was |
|
174:43 | trick question. So when you um uh do the uh the um |
|
174:51 | the alert for trick questions like OK. OK. Please. So |
|
174:59 | here's the next, even though um anti sci is usually its effect on |
|
175:05 | Avio curvature is usually large. Is true or false? I think it's |
|
175:17 | . Well, I'm gonna say it's because it says usually, and, |
|
175:21 | , you know, we don't have experience. The set on Avio curvature |
|
175:26 | usually very noisy. It's usually so . We don't even look at it |
|
175:29 | all. So, uh who, to say, what is the effect |
|
175:33 | anisotropy? I would say that this uh uh uh uh this is false |
|
175:39 | stated, although it might be we might learn this true. And |
|
175:44 | fact, there's another set of phd to learn how to use a vo |
|
175:49 | uh curvature. OK. Next as you uh uh if you have |
|
175:57 | estimate of delta as a function of for move out, you can use |
|
176:02 | to estimate the jump at the target . Is that true or false? |
|
176:12 | gonna say that's false because usually this is has low spatial resolution, |
|
176:18 | So you can't get a local jump a smoothly varying curve. Mhm |
|
176:24 | So uh uh uh from uh you , you just don't have the spatial |
|
176:31 | , you need to get D and we could, I would have solved |
|
176:34 | 3040 years ago. But uh uh , that's a false. OK. |
|
176:40 | , so much for polar and So now we're gonna deal with auth |
|
176:46 | symptom. And the first thing I is that uh you're probably familiar with |
|
176:54 | model called tilted polar on isotopy. here's a picture of a tilted mat |
|
177:03 | you can be sure that there beneath joints, there's um uh uh layers |
|
177:10 | uh the layers are uh uh uh originally polar an isotropic. And now |
|
177:16 | whole thing has been tilted. So common uh a model even today is |
|
177:23 | uh uh you assume tilted or an pain. But that's probably not plausible |
|
177:32 | you know, the same stresses which the dip probably opens up fractures. |
|
177:38 | so these uh here, here's a example, these fractures, these joints |
|
177:44 | parallel to the strike of the uh of the structure. That's not an |
|
177:50 | . These uh joints were opened up the same stresses which caused the structure |
|
177:56 | the same in the um in the place. If they were oriented at |
|
178:01 | random angle with respect to the you'd say that it's just an |
|
178:06 | but these are lined up parallel to strike. So what that shows is |
|
178:12 | uh uh probably a a common So these beds here have tilted or |
|
178:18 | symmetry. Why? Because these joints on the strike of the structure. |
|
178:26 | furthermore, so-called ht I horizontal transverse is never plausible in the sedimentary |
|
178:34 | And here's the reason hori horizontal transverse could come from a single set of |
|
178:43 | circular fractures in an otherwise isotropic And in fact, this was the |
|
178:49 | that we first had when we first thinking about is and isotropy uh uh |
|
178:55 | years ago. Uh But here's the , there are, there aren't any |
|
179:00 | fractures anywhere. In fact, the ones which we're most interested in are |
|
179:04 | shaped like we showed in the previous . Uh Furthermore, in the segmented |
|
179:10 | of the background is always an So both of these um assumptions are |
|
179:17 | . So uh we should um not a TIHT I assumptions today. Uh |
|
179:24 | were suitable when I first came into business, but not today. |
|
179:32 | it would be nice if we could fracture. But it turns out that |
|
179:38 | we can't, the most realistic approximation that of orthorhombic an isotropy. So |
|
179:46 | is an air photo of a, place in um uh Southwester United States |
|
179:51 | you can see the laying. you know, it's gonna be uh |
|
179:55 | a uh uh gonna have floss which with the incident angle and you can |
|
180:03 | the joints and you can see the other joints and there's a good |
|
180:07 | why um uh uh we commonly see joint sets like this. I don't |
|
180:15 | we have time today to talk about that is, we don't always see |
|
180:21 | , but we frequently see it. it's obvious in this uh in this |
|
180:26 | here. So in that case, uh the, the symmetry matrix is |
|
180:33 | complicated. This is what we had polar symmetry for un fractured shells and |
|
180:39 | bed sequences. Now, for uh Robic Symmetry with uh which shows fourth |
|
180:48 | sequences with either one or two vertical crack sets. If the Ortho to |
|
180:55 | other, then that's gonna lead to Robic Symmetry. And that has um |
|
181:00 | uh an anti stick. This matrix like this still has a lot |
|
181:06 | of zeros off here. But we nine independent quantities in here instead of |
|
181:15 | . So for orthorhombic formation, it the, the symmetry of a |
|
181:22 | So uh uh it has three different symmetry directions uh perpendicular to each |
|
181:30 | And uh notice in this graph, have right handed cord system, uh |
|
181:36 | lay your uh your pro your index along the X one direction, curl |
|
181:41 | over uh to make it into the two direction and look which way your |
|
181:46 | is pointing. It's the downward X direction. If you're using your right |
|
181:52 | , if you're using the left it would be the opposite. So |
|
181:58 | the equations for velocity in terms of two angles, which is the angle |
|
182:08 | incidence here and the an angle those equations are very complicated. If |
|
182:16 | uh uh the first time I, looked at them on the screen, |
|
182:19 | printed them out and it, it about 10 pages of output for one |
|
182:24 | for the, for the P wave . It was just amazingly complicated. |
|
182:31 | I gave up but my friend, , whose name you saw earlier uh |
|
182:37 | than me. And he didn't give . And he looked at those equations |
|
182:41 | thought about them and he realized that in the case where um if you |
|
182:50 | at only um propagation in this one it, it has the those complicated |
|
182:59 | reduced to the equations for polar anisotropy approximation. And that's true, not |
|
183:08 | for that this plan, but you know, a any plan parallel |
|
183:11 | this one has the same um um same polar anisotropic media. And we |
|
183:18 | how to handle all ANAs tropic And furthermore, if you looked at |
|
183:23 | o this other plant also true, it's a different uh uh uh a |
|
183:29 | um or an a striping medium for plant than for this one. So |
|
183:39 | um look to see what we uh uh let's parameterize that system for those |
|
183:46 | planes. So in the one plane have uh uh a P wave, |
|
183:51 | P wave velocity, vertical shear wave and three antihero parameter here. And |
|
183:58 | the other plane, you have the vertical uh uh P wave velocity and |
|
184:04 | a different shear wave velocity and uh more anisotropic parameters. And it looks |
|
184:12 | you have nine parameters in all you know, four plus four plus |
|
184:17 | makes nine what you need, but is smarter than that. And he |
|
184:21 | carefully at these at the, and saw that within this set of, |
|
184:26 | uh parameters. You uh don't have 13, excuse me. Uh You |
|
184:34 | have AC 12, here's AC here's AC 23, but you don't |
|
184:39 | AC 12 anywhere. So he made uh out of whole cloth, he |
|
184:44 | up another parameter which has a functional like this except it's got AC 12 |
|
184:50 | there. And so now we have set of nine parameters with uh which |
|
184:56 | all the hook in parameters in And we know how to uh handle |
|
185:02 | anisotropy. In that case, for , the P wave anisotropy has a |
|
185:08 | phone. It looks exactly like we before. Except that look here, |
|
185:14 | anti sucky parameters have an auth dependence them. So if you want to |
|
185:20 | what that is and this is what is. Here's the azimuthal dependence of |
|
185:26 | , it's got delta in one of planes, one of the vertical planes |
|
185:30 | the trigonometric factor and uh plus delta uh the other plant, another trig |
|
185:36 | factor, et cetera. And so here that this is the uh the |
|
185:43 | , this is an elliptical as an variation. So if we look at |
|
185:50 | offsets, we're gonna see only elliptical of delta. OK. So let |
|
186:00 | show you this in real data. is data from Lynn et all in |
|
186:06 | . This is a different len than one that I just mentioned that uh |
|
186:10 | L also a woman um was an colleague of mine back in the |
|
186:15 | And then she left Amaco and she up a company of her own called |
|
186:20 | Incorporated and you can hire them. uh I think you can just uh |
|
186:26 | Lyn Incorporated and you can hire them do specialized data processing even today. |
|
186:34 | , um for all this time, , she's made a living uh as |
|
186:38 | boutique processing house. And so this set comes from a study she did |
|
186:44 | 1999 for the, the US Department Energy and it's a wide Asus |
|
186:51 | And so uh uh look what she's . She's taken the data over |
|
186:56 | this is a time series of data she's made a narrow and it's wide |
|
187:02 | data with all asthma represented in this set. And she took uh uh |
|
187:08 | took a small uh asthma sector around degrees of asthma from zero from 0 |
|
187:15 | 10 degrees. And that's plotted here is uh about 20 degrees, 30 |
|
187:21 | and all around here to 180 And then of course, it wraps |
|
187:26 | the rest of the compass, wraps it because of the scalar theorem of |
|
187:32 | . And so I want you to on this um uh uh event down |
|
187:39 | and I want you to ignore the hyperbolic move out here and, and |
|
187:45 | on, she, she's applied uh short spread hyperbolic move out to this |
|
187:54 | only that she wasn't interested in And so uh you can see that |
|
188:01 | right here, the uh layers are . But if you go around the |
|
188:09 | , you see they get to be and less flat and here they're, |
|
188:13 | straight but they're not flat. And you keep on going around the circle |
|
188:18 | it gets better and better and then comes back here to flat. So |
|
188:23 | all of these are flattened with the uh same applied velocity function. And |
|
188:31 | see the difference between where it's most and where it's most incorrect, that's |
|
188:37 | degrees difference. That's the 90 degree in move out velocity as a function |
|
188:46 | as muscle uh variation right there in of your eyes. Nothing difficult was |
|
188:53 | to uh to show this. Uh know, if you do some kind |
|
188:58 | migration with uh uh anti sci beer , or whatever, it's always a |
|
189:04 | business and you never know where you have made a mistake. But |
|
189:08 | basically, all she did was sort data into subsets and apply an AIC |
|
189:14 | out. And so she couldn't possibly made a mistake. So because of |
|
189:21 | previous formula about um uh let me up here. Yeah, because of |
|
189:28 | formula here, we can use that to understand what happens between the |
|
189:38 | Now, let's uh uh think about . We're looking at um of the |
|
189:46 | of orthorhombic anisotropy on seismic. It's complicated, I would say, not |
|
189:54 | . But it's uh um uh it's when I was a boy back in |
|
190:01 | 19 forties. My father was finding in uh East Texas as a geophysicist |
|
190:10 | he didn't know about any of He didn't know about Ortho Robic |
|
190:16 | He didn't know about anti soy. didn't know about anything and he was |
|
190:20 | finding a lot of oil and it's our while to ask ourselves, how |
|
190:24 | those guys back in the day? they were so ignorant, they weren't |
|
190:29 | but they were ignorant. So many we know today, they didn't know |
|
190:36 | and yet they found a lot of . How could that be? |
|
190:39 | here's a partial answer back in that , he was doing one D size |
|
190:44 | , excuse me, two D size of a single line of acquisition. |
|
190:49 | furthermore, he was using short So he didn't know or care about |
|
190:56 | hyperbolic. Move up. Here's a of his line, his uh two |
|
191:04 | acquisition line right here, going at angle uh uh with respect to the |
|
191:09 | Ortho IIC symmetry. He knew nothing this, but he's assuming that the |
|
191:18 | are gonna go down flat line country East Texas. So he's assuming the |
|
191:24 | are, are uh uh uh the are gonna lie in the vertical plant |
|
191:29 | this acquisition line. Well, he know it, but you see the |
|
191:33 | rays were actually bouncing off to the and coming back into his receivers. |
|
191:39 | didn't know that all he knew was had uh uh to the light flat |
|
191:45 | uh geometry and he said, I for sure. Now the rays are |
|
191:49 | be lying in this same vertical plank that whatever images he was making were |
|
191:58 | , not too much, but some bouncing at the mid, near the |
|
192:03 | here. Um And he had uh hyperbolic move. So here's this hyperbolic |
|
192:15 | and he didn't know it. But there was a uh an anti parameter |
|
192:20 | which he didn't know about. And was the, the one which is |
|
192:25 | for this angle. So he didn't about this separation between anisotropy and uh |
|
192:32 | . He didn't know about the Asmal variation because it wasn't, it wasn't |
|
192:36 | his data and he used these numbers convert time and depth and he got |
|
192:44 | this time, but everybody got his . And so he didn't worry about |
|
192:48 | much uh found a lot of This is why guys in, in |
|
192:56 | era were so successful despite the fact they were so ignorant compared to what |
|
193:01 | know today. Now, let's think azimuthal variation of Avio uh uh move |
|
193:10 | , move out velocity, not move out velocity. Uh This is |
|
193:13 | problem that we, that I showed uh that Lynn was uh looking |
|
193:18 | So this is our one D tailor . And for uh if we have |
|
193:23 | variation, if we have uh uh 3d survey, uh oh uh land |
|
193:30 | , uh marine survey, whatever we're have uh uh offsets in the X |
|
193:35 | , offsets in the Y direction. uh Mr Taylor says that if we |
|
193:40 | first order Taylor expansion in small X small Y, you get three terms |
|
193:46 | this, you get an X square , you get a Y square term |
|
193:49 | you also get an xy term. that is the equation for any lips |
|
193:55 | in the uh in the plane of map. So you never get any |
|
194:02 | in the vertical plant. But you got any ellipse in the plane of |
|
194:05 | map if you restrict your yourself to offsets. And you can see there |
|
194:11 | 123 parameters to be determined here and , they've been given uh complicated names |
|
194:18 | , but you can see there's only parameters to, to determine. So |
|
194:23 | me show you um OK, some to illustrate this point. So this |
|
194:33 | a cartoon showing a single uh uh with hyperbolic move out from wide estimate |
|
194:39 | . And uh uh uh you can there's some uh a scatter about the |
|
194:44 | the best fit hyperbole here. And offsets are represented here uh uh |
|
194:49 | in this and the little pics here give the, the tops of the |
|
194:55 | peaks in the survival. And so looks like a bunch of scatter. |
|
195:01 | We did the best job we could this uh average. But if you |
|
195:05 | closely, you can see that the from the north and the south come |
|
195:12 | early and the arrivals from the east the west come in late. And |
|
195:18 | can you see that the uh uh discrepancies at far offsets are bigger than |
|
195:24 | discrepancies at their offsets. So these patterns in the data. And whenever |
|
195:29 | see patterns in the data, uh you know, your uh patterns in |
|
195:34 | residual, you're look, you you're looking at uh data, not |
|
195:39 | noise. So when you remove uh that hyperbola is flattened uh |
|
195:48 | then it looks like. So and this looks like it's been uh |
|
195:53 | doesn't it? But I move out uh NMO stretch, but it's |
|
195:57 | you see, uh it's the thighs basically of my um uh pointer. |
|
196:03 | Now I'm gonna go back and you see it's the same size here. |
|
196:07 | we haven't stretched anything. We that's just the way it looks. |
|
196:12 | more clear to the uh to the to the eyeball. Now that we've |
|
196:18 | the best fit hyperbole. No, good way to look at this kind |
|
196:28 | data is to look at the not by mixing all the offsets |
|
196:33 | So this is what we call an uh ordered gather. But uh uh |
|
196:42 | I'm gonna show you what happens if sort the gather by aims. |
|
196:47 | before we sort it by. So is from 0 to 100 and 80 |
|
196:51 | . So the, so before we it, we wanna limit the |
|
196:54 | So we have only uh uh say um uh uh 9000 m and 10,000 |
|
197:03 | offset between 8500 m off. So somewhere range limiting and then an animal |
|
197:12 | it. And then the uh you're likely to see a um uh |
|
197:16 | variation like this. These are the uh picks that you saw before, |
|
197:22 | off of the um uh top of peaks. And number one, you |
|
197:28 | that not all asthma are recorded, always happens that when you do a |
|
197:32 | asthma survey, some as are not . So here we have some blank |
|
197:39 | then the events here are not uh . And so all the uh we |
|
197:47 | then uh uh we, we see the, the peak right here is |
|
197:52 | degrees away from the trough. So can remove it with the uh with |
|
197:56 | parametric residual move algorithm. Um All do is uh uh estimate those uh |
|
198:04 | uh uh the amount of that ellipsoidal move out velocity. So we uh |
|
198:13 | the uh slow velocity, the fast and the, and the auth uh |
|
198:19 | the orientation of the ellipse. So me show you what this looks like |
|
198:23 | real data. This comes from Val , talked about that before val. |
|
198:32 | comes from um off to the side the uh the data quality is |
|
198:36 | This is the top of the, the reservoir right here. It's obviously |
|
198:42 | flattened very well. But if you closely, can you see this jitter |
|
198:49 | arrival times, nearby traces have uh uh have uh differences in arrival times |
|
198:56 | several uh milliseconds, maybe a 10 or more. Uh uh And then |
|
199:03 | changes back and forth in a, rapid variation. You can't imagine any |
|
199:10 | of, of one B structure which gonna uh do, do that. |
|
199:17 | so um this is a conventional gather , sorted and flat. So now |
|
199:25 | gonna do it as, as sorted limited and now we put it as |
|
199:31 | function of Asma. So number you see that there are some places |
|
199:35 | there's no data. And number there are places where there's early arrivals |
|
199:41 | then nearby places which are later So I thought this, this uh |
|
199:47 | ship, I'm gonna go back, , look at this place when I |
|
199:51 | off the white LA. And so heavy black um uh uh first arrival |
|
199:59 | in right here and then it comes later over here. And so um |
|
200:12 | what this the 90 degrees difference in . This is exactly what Lynn showed |
|
200:19 | her data set. Now, I'm show you as a mostly an isotropic |
|
200:25 | vo. So to make it we're gonna use only the two terms |
|
200:31 | , of a Avio gradient and uh 22 terms of, of Avio |
|
200:37 | We're gonna ignore this four third term here. The gradient term has all |
|
200:42 | features in it. It has a jump in uh uh vertical uh |
|
200:48 | VP it's got um uh a jump vertical shear modulus. This is modified |
|
200:56 | from what you're familiar with. And it's got a jump in delta as |
|
201:00 | function of asthma. And uh so is where the differences are from uh |
|
201:08 | what we saw before. So in of verbals, let's just uh um |
|
201:13 | all the uh um the terms which as smoothly verbal out of here and |
|
201:20 | them um uh into one term. inside here is a bunch of stuff |
|
201:25 | I don't want to show you But uh uh it's at a single |
|
201:30 | with this kind of as a little . So now let's look at some |
|
201:37 | um white as with data. And this is uh uh an amplitude data |
|
201:47 | uh uh this is the kind of curve that you're accustomed to looking |
|
201:51 | You see, it's got an intercept and curvature and you also see a |
|
201:57 | of scatter. So this is a three term a curve. So uh |
|
202:05 | here up here is another data point not part of the legend, that's |
|
202:10 | data point way out. Uh uh uh couldn't possibly think of any physical |
|
202:17 | for it to be there. So is a procedure in statistics called robust |
|
202:22 | fitting, which enables you to eliminate like this. But don't eliminate terms |
|
202:28 | this using uh a standard uh criteria deciding uh we're gonna just throw this |
|
202:35 | , but we're not gonna throw out one. So you ask your |
|
202:41 | the question is this variation of the or is it signal the way you |
|
202:49 | that question is you look for pattern the variation patterns in the residual. |
|
202:56 | one way to look for it is rep plot this data in terms of |
|
203:01 | . So now you see it has coherent variation with asthma, there's still |
|
203:07 | lot of scatter, but you can see a high zone and a low |
|
203:11 | here. And now let's look at same real data set. Well, |
|
203:15 | is real data from elsewhere in, uh vowel hall. But uh uh |
|
203:20 | look at the same vowel hall data that we just saw uh asthma, |
|
203:25 | range limit range with asthma. And you see that uh the uh uh |
|
203:33 | bright amplitudes are here and the dim are here 90 degrees separation. That's |
|
203:44 | , that's a, an, an variation in amplitude just like the elliptical |
|
203:51 | and move out that we saw a slides ago. So let me look |
|
203:57 | a map view here um of that effect. So this is about 505 |
|
204:04 | m by about 2000 m, something that. And you see in, |
|
204:11 | each pixel of the map, this a map view in each pixel, |
|
204:16 | a little arrow without an arrowhead. so the length of the arrow tells |
|
204:20 | the amount of Ainu variation ra and orientation of the arrow indicates the uh |
|
204:29 | hot direction. And so the same variation is uh shown in color. |
|
204:35 | uh all the the bi the big variations is in red, see down |
|
204:40 | in blue hardly in amplitude variation and scale is given over here. So |
|
204:44 | red, it's 200% variation. We're talking about small effects here, 200% |
|
204:53 | in Avio gradient. So maybe you see a right angle here or maybe |
|
204:59 | . Uh I always look for right in this kind of data and I |
|
205:03 | see them, maybe I'm fooling who knows. Uh But you can |
|
205:08 | that if you were doing a two survey along this line here, as |
|
205:12 | as you came into this area, would immediately realize that you had entered |
|
205:18 | Avio anomaly area just looking at the . But if you came at it |
|
205:26 | this direction, you wouldn't see a because all because uh these gradient is |
|
205:31 | vector quality, all the variation is this northwest direction. There's no, |
|
205:36 | no variation in the northeast direction. this, this line is perpendicular to |
|
205:42 | these other, these other lines So uh uh let uh this is |
|
205:48 | small study. Let's look at a study. This is about five kilometers |
|
205:53 | about uh 15 kilometers. And uh can see uh it's the same uh |
|
206:01 | bar basically, except that the maximum white instead of red. And also |
|
206:06 | can see a dark area in here no variation that's corresponds to these navy |
|
206:14 | um uh colors down here. and it means that where we didn't |
|
206:20 | the data or the data was so that we didn't believe it, we |
|
206:25 | colored it navy blue and that was clever thing to do because it took |
|
206:32 | away from our eyeball and then all colors that we can see are gonna |
|
206:36 | things that we do believe in. , I gotta tell you about these |
|
206:41 | uh red lines here. This is uh uh BPS, um permanent insulation |
|
206:50 | receivers. They buried 10,000 seismic receivers the sea floor about a meter down |
|
206:58 | the mud along these lines, connected all up to uh uh uh to |
|
207:03 | central platform which lies right here. then from there, they sent a |
|
207:08 | cable to the shore. And so they were a uh able to uh |
|
207:14 | these things, uh uh uh you , analyze them onshore. The reason |
|
207:21 | this expensive installation of c ocean bottom was so they could do cheap rapid |
|
207:34 | D reshoots. Think about this, you do an ordinary um uh four |
|
207:39 | survey, you do a 3d survey . That's your baseline survey. It |
|
207:44 | cost to five or $10 million and uh you come back about a year |
|
207:49 | and do another one. And it costs five or $10 million and you |
|
207:53 | never do a third one because it's expensive. But BP had the idea |
|
207:59 | if they spend uh um uh $20 up front to bury the receivers in |
|
208:06 | sea floor, then the seismic reshoots be very cheap, maybe $50,000 |
|
208:12 | Why is that? Because they they come back with a uh uh |
|
208:17 | of a large um uh uh seismic vessel with 10 kilometers of cable stretching |
|
208:24 | uh you know, on a small with only a source behind it because |
|
208:28 | receivers are already in place. And with that, that's a very good |
|
208:33 | . And they've done now over 20 and over this field uh uh using |
|
208:38 | clever idea. And uh furthermore, can see that the, the, |
|
208:44 | can see that platform is right here the middle and we, we only |
|
208:51 | instrumented half of the field, there's half of the field over here that |
|
208:55 | didn't instrument, obviously, the platform gonna be in the middle of the |
|
209:00 | . Uh uh Because when we did , it was expensive and we, |
|
209:04 | weren't sure whether it was gonna be . So then after we operated this |
|
209:09 | of the field for a while, we realized this is very useful. |
|
209:13 | we uh instrument the other half. that happened actually that after I left |
|
209:19 | , so I don't have any, data from there. Um So notice |
|
209:26 | values are only shown where the consonants the fit is more than 95% of |
|
209:32 | . So let the maximum values show than 100% difference in um avi O |
|
209:40 | as a myth damage. So let's in on the northwest corner. You |
|
209:43 | see areas here like in here where little pixel is calculated independently and uh |
|
209:50 | has a consistent direction, then there's no uh a narrow uh transition zone |
|
209:56 | then another area over here again, uniform with a different orientation. And |
|
210:03 | that's typical of uh of uh fractures my experience is that they're locally uniform |
|
210:12 | a narrow transition zone. And then else on the other side of the |
|
210:16 | zone, let's look at the southwest , the same corner thing, a |
|
210:22 | at some data that goes into So this shows the set of uh |
|
210:27 | a for a given point a given in that map. This shows the |
|
210:32 | of asthma and, or, and offsets that went into that. And |
|
210:36 | shows uh the uh Avio gradient. can see there's lots of scatter but |
|
210:43 | a, a well defined uh sinusoidal that goes through there here shows um |
|
210:54 | another case with lots of scatter. uh here, the scatter is so |
|
211:01 | that you don't believe the variation. is not a straight green line |
|
211:05 | This is a, a Sinusoidal um uh variation with a small uh a |
|
211:12 | variation, but we don't whatever it , we don't believe it because of |
|
211:16 | large value amount of scatter. And are standard statistical ways for you to |
|
211:22 | when you want to believe it or . This is the same data set |
|
211:28 | a different color bar. And we that to our colleague in Norway, |
|
211:33 | knowing what these colors meant. And sent us back this curve. The |
|
211:38 | thing about this curve is that he also the uh central platform with the |
|
211:45 | uh the boreholes. They go straight from the platform and then they deviate |
|
211:49 | all corners of the field. And see that every bore hole ends up |
|
211:54 | a patch of color and there are patches of color except for this one |
|
211:59 | don't have a borehole in them. what that means at the uh what |
|
212:05 | means uh you gotta realize at the of this bore hole, they perforate |
|
212:09 | lining and they, I inject um fluid or they produce uh hydrocarbons from |
|
212:16 | perforated part. So that when you the perforations in like, so that |
|
212:21 | that this, these colors were caused those operations in the field over the |
|
212:27 | 15 years, uh uh uh uh activity. And when we see them |
|
212:34 | in the uh uh surface seismic now take a look at this |
|
212:42 | Oh I should be operating uh by way, this is the sec. |
|
212:45 | is the second survey. This is first one, they tell mostly the |
|
212:48 | story. We'll come back to the between those in a second. You |
|
212:52 | this one missed it. It's very that uh when they put this well |
|
212:56 | here, they were expecting to produce a, an area around like. |
|
213:00 | instead all the changes came from What that means is that there's a |
|
213:05 | barrier in here somewhere which which prevents fluids from this side of that bore |
|
213:11 | from entering this uh borehole. uh what that means is there's an |
|
213:16 | uh uh there's a exploitation opportunity put another borehole right in here and |
|
213:22 | produce this all multimillion dollar uh uh to be made from high tech to |
|
213:30 | f So now let me show you difference between these two. So those |
|
213:35 | were taken uh uh with uh uh uh the time difference here between uh |
|
213:41 | three months time, between here and . And you can see that these |
|
213:44 | the differences, see the differences all in the same places where the uh |
|
213:50 | just grow over time. So that the uh one of the uh the |
|
213:55 | fidelity that you can measure these things high temporal resolution as well as vertical |
|
214:02 | . So folks, there are many topics uh uh that I uh would |
|
214:08 | to talk to you about, but have run out of time today. |
|
214:12 | so um I wanna uh uh uh leave you with this with this result |
|
214:21 | this study. Uh uh Today, have talked about three sorts of oo |
|
214:29 | uh effects from seismic anisotropy. All which come from the fact that the |
|
214:36 | and the kinds of rocks we explore is usually small, less than |
|
214:42 | Nonetheless, there are three types of you might see of that we call |
|
214:49 | second order effects, which are also of the order of 10% like in |
|
214:53 | move out velocity, they're not but they're small of the order of |
|
214:57 | . And there are large effects, effects. That comes for example, |
|
215:04 | ale from the same small anisotropy block , you get a 100% you can |
|
215:13 | a 100% effect difference in the A Grady and then there are zero order |
|
215:19 | which are completely new, completely um seen at all in isotropic seismic. |
|
215:25 | that includes, for example, sheer splitting, which we didn't have a |
|
215:29 | to talk about today. So with , what we're gonna, what um |
|
215:34 | we're gonna do is uh um all to uh uh close and uh since |
|
215:42 | is out, I have to hurry of here. Um But what I |
|
215:47 | do when I get home tonight, gonna send you all uh in the |
|
215:52 | exam and it will be in the of a, of a, of |
|
215:58 | attachment to an email. And so you do is uh uh you, |
|
216:03 | look at that and set it aside you have three hours to look at |
|
216:08 | and then you uh uh uh open attachment. My advice is to print |
|
216:15 | off on paper so that you can in your answers by hand and you |
|
216:19 | have to worry about the logistics of uh entering stuff by the computer. |
|
216:25 | then once you uh uh and it's be uh unlimited time, open |
|
216:30 | open notes. Um So it's a of your understanding, not of your |
|
216:37 | , not of your mathematical, not of your understanding. You will turn |
|
216:42 | this in before midnight on Wednesday, turn it into Utah when he collects |
|
216:48 | all, he will send them to . Um your eye um you will |
|
216:54 | unlimited time, open, book, notes, no consultation with anybody |
|
217:00 | So se uh, set aside some at least three hours when, |
|
217:06 | you, uh, will be undisturbed you did this. So, with |
|
217:11 | , uh, I should say there's a pleasure to teach you and I |
|
217:15 | you do well on the exam uh, I will grade the exams |
|
217:21 | them and get them back from your . Oh, stop sharing right |
|
217:27 | And uh, so for you folks remotely, uh I'm going to say |
|
217:33 | and thank you. I will um uh expect to hear from you via |
|
217:40 | by Wednesday from me later tonight after . Thank you. Thank you. |
|
217:46 | you. Bye |
|