© Distribution of this video is restricted by its owner
Transcript ×
Auto highlight
Font-size
00:00 Recording. OK. So backing up backing up, we uh uh we

00:23 with this conclusion because of the previous , we know that uh uh uh

00:33 pre velocity uh for a rock containing gas is the left oops I need

00:42 pointer. OK. The P velocity rock K and KG is less than

00:51 P velocity for rock, same rock containing Bryant. And uh uh of

00:57 , uh the, the text down , it reminds us that uh uh

01:02 , to calculate, to calculate this gas dependence here, uh You have

01:07 know about properties of hydrocarbons under high and temperature and that's all well known

01:13 uh from uh many uh years of work. And for example, uh

01:19 rock properties uh grew up here at University of Houston has been active in

01:25 for many years. They know how tell you uh what are the,

01:31 in compressibility of any mixture of And uh um uh uh also when

01:40 uh partially saturated with Bryant. So that understanding, uh what we want

01:45 do is uh uh point out that um understanding for velocities leads to anomalously

01:53 reflections and to anomalous A O. , uh the anomalous black reflections um

02:03 uh uh uh a, a predecessor a vo when I first came into

02:09 business. Uh The, the only direct detection that we had was what

02:14 call bright spots, which is what now call a VO as expressed in

02:20 stack. And o we were only to trace that down to the angular

02:26 with offset uh shortly after I arrived the business. And because of

02:31 we can, we can have uh 40 seismic surveillance. And uh of

02:36 , that this makes a big difference the economies of our exploiting these

02:43 So uh I want to uh uh uh briefly what we talked about before

02:52 that was before we had this knowledge the effects of fluids. So remember

02:57 talked about for normal set of inter , we have this situation for the

03:03 vo intercept and the gradient. And on laboratory data, we know that

03:08 term uh uh dominates over this And because of that, that uh

03:15 gradient term has an algebraic sign opposite that of uh the instant media.

03:21 that's positive, this is negative and versa. Also for an interface between

03:28 brine and a gas on uh this like that uh uh no changes

03:34 in mythology, just a change from to gas. Maybe we're talking about

03:39 uh the top of uh uh of partially filled reservoir, same mythology above

03:47 below the reflection in this case, you know, that's a special

03:51 an extreme case and any real case be a combination of this case and

03:56 previous case. So uh uh uh at this case only we have the

04:02 expression, very same expressions for intercepts . But now we know that this

04:08 is zero from poor elasticity. That's we learned. Uh we saw this

04:13 back in electron six. Now you're the middle of uh lesson eight.

04:17 now you understand that this comes from essentially in 1941. So in that

04:23 , the algebraic sign of uh this a zero. So this always has

04:27 same algebraic sign as this one. uh that is what we call anomous

04:33 vo. So uh as we said , the real world may be more

04:38 , but there's lots of um um now with Avi L and it all

04:44 depends uh this is the essential idea there's many complications um um that those

04:52 deal with. That's that course not course. So here is a,

04:56 couple of uh of um quiz It says uh let's see. Uh

05:06 to Gasman, the fluid dependence of launch shoulder mous is given by

05:11 Is this true or false? Let turn to um uh Meade for this

05:16 this one true or false. I it is true. Yes, it

05:28 true. But that's not what I you. Remember, I told

05:31 I, I gave you a funnel this with Ks on the left.

05:35 uh yeah, but even, but answer is correct, I taught you

05:40 difference with cases. But you you know that this m is equal

05:44 K plus four thirds mu. So just add in your mind uh plus

05:49 thirds mu and plus four thirds new . And that difference between the frame

05:56 and the uh uh moduli for sheer zero. So that putting in that

06:03 modification makes no difference on the right , it's still true. So,

06:12 good. Uh Let us uh look uh first question. Uh uh Number

06:20 , it says uh the assertion that classic gas monetary needs to be

06:25 that's the assertion that I told you that I published just last year requires

06:30 experimental confirmation. Um uh Lee would you say that's true or

06:37 Yeah, I think it's true. So I'm, I'm the one who

06:41 out the error. Uh But even say maybe that error is important and

06:45 not. Uh Yeah, uh we'll out was it good for experiments

06:51 So, question number three, uh goes to you Carlos using either classic

06:59 line theory or its refinement as discussed morning, the presence of gas in

07:04 pores can, can lead to significant in P wave velocity and impedance.

07:10 that true or false? Well, think it's true. Yeah, I

07:17 call that true. Yeah, that , and furthermore, uh that has

07:23 uh responsible for a lot of uh enormous amount of discoveries of hydrocarbons.

07:31 tell you this story here. Uh uh uh two parts of the

07:35 Uh my father was the one who uh bright spots back in the early

07:45 forties. Imagine that about the time was uh uh born, he was

07:52 bright spots and he was, he a field geophysicist. He was uh

07:58 the crew and the crew would move and uh uh uh the office crew

08:03 follow the field crew and uh he look at these primitive records coming to

08:11 , you know, just wiggles. didn't have any computers, he didn't

08:15 any workstations. Uh They would um the data on photographic film during the

08:22 . They would develop it overnight in morning, he would look at it

08:26 he would say, oh well, we need to uh uh modify the

08:33 plan for uh today's uh um you , acquisition will uh go this way

08:40 of that way. So uh that's kind of conditions he was uh uh

08:46 under uh by current san, you , hopelessly primitive. But he,

08:52 , he noticed that uh whenever he make a recommendation for drilling. If

08:58 recommendation came out of bright reflections instead dim reflections, uh, it was

09:04 likely to be successful. So he fabulously successful as an oil finder.

09:10 found, uh, he, uh, his record for,

09:14 recommendations was that, uh, when recommended a, a well, there

09:21 a 25% chance that, that well gonna be good. So today,

09:26 think 25%? That's terrible. Uh If, if we had that kind

09:31 work to today, we would all fired. And that's true for

09:34 But in his day, that was uh uh you know, very rare

09:39 find a, a geophysicist who could with 25% accuracy where to drill.

09:47 he, he became a famous guy Amal and that's why they hired me

09:52 years later. So he, uh that was his story of developing

09:58 of, of inventing bright spots. The, the next stage of,

10:03 , of that. Uh I, gotta tell you more about that.

10:07 , he uh he made these recommendations of course, he was ignored for

10:12 years by Amaco Management. Uh What said was uh uh uh you

10:19 Mr Thompson, uh uh your job just to uh examine the wiggles.

10:24 there's any important new ideas to come , that'll happen in our research center

10:29 Tulsa. If you'll just please follow recipe, follow the recommended procedures,

10:34 would appreciate that. So he was with that kind of a management brush

10:40 for years and years. Uh later uh discovered uh the same thing.

10:47 uh here's how we learned about Uh It was the early days of

10:53 exploration in the Gulf of Mexico and uh the way it works in

10:58 is uh those uh mineral rights for uh um parcels are owned by the

11:06 government. And uh the, the government uh puts up for auction

11:12 year. A few parcels and oil bid on the rights to explore and

11:18 in those parcels. And people began notice that on certain parcels uh a

11:24 the, the they're put up for and uh there's period for uh for

11:29 uh exploring and you can e either with proprietary equipment or explore with the

11:36 companies. In those days, there a lot more done by all companies

11:41 . Uh And it was two D uh short streamers o over the uh

11:47 in the Gulf of Mexico. And began to notice uncertain um prospects.

11:55 was laying down large bit bets, bits and so they would bid maybe

12:00 million dollars. It was a big when everybody else was bidding uh uh

12:05 30,000. So uh uh it was clear that the Mo Mobile knew exactly

12:13 they wanted to drill, they knew sort of secret that we didn't

12:17 So eventually that secret leaked out and was bright spots. The predecessor of

12:22 BO in those days, we didn't enough uh spread links to do a

12:29 uh uh measurement of the, of amplitude variation with offset all we we

12:35 at was the, the stack So, uh that's where bright spots

12:40 from. And uh so, uh was due to uh this um my

12:49 here. uh The, the gas theory is uh uh tells us that

12:58 you have gas in the forest that lead to significant reduction in P wave

13:02 and impedes hence brighter reflections. uh so uh I talked over you

13:11 Carl. So let me give this to you Carlos, it says uh

13:15 or false using either plastic gas mount , quartz, modern refinement discussed

13:21 the presence of gas in the phosphates lead to significant reduction in P wave

13:26 S wave velocity who are false. here, it says it specifically is

13:35 the, that's correct. That's the . Yeah, very good. So

13:40 gonna call this one false and, uh uh uh that, that would

13:44 correct. Uh It affects the P , not the S waves.

13:49 remember the S waves have in the density, they have the sheer

13:54 uh divided by uh the density and density will be affected by uh the

14:00 . But um uh this isn't about den density, this isn't about sheer

14:06 , it's about sheer velocity. And that's true. So, uh that's

14:13 good news here. And I'm thinking all thinking, well, what's the

14:17 , what's the big deal? What found out was that uh when we

14:20 from elasticity on homogeneous solids to poor on, in porous rocks, it's

14:28 making a, a minor difference. we learned is that um uh moduli

14:35 the density depend upon uh yeah, over the uh uh the heterogeneity that's

14:47 in Iraq. And in particular, sure dependence is given uh in a

14:52 way by vo and it's uh in and that's refined in 2023 72 years

15:02 . Um uh But they both lead the same conclusions here. With regard

15:08 exploration, we will learn about the . Uh whether or not that's an

15:13 refinement in the next few years. now I want to uh uh return

15:20 the idea of bo slow waves. I told you that Bo taught us

15:26 for porous rocks, we can expect kind of wave. So that's called

15:32 bo slow wave. So we got only P waves and sheer waves,

15:36 we have a mo bo slow And so uh uh that slow wave

15:43 because of the heterogeneity. And the he heterogeneity is we got truths versus

15:50 . Yeah. Right. So, way back in time, uh uh

15:54 uh Max Barn was a famous physicist in, in those days. And

15:59 asked himself the question of what happens you have uh uh a one d

16:05 uh object like this with beads, heavy beads uh separated by springs.

16:11 so I don't know why he would about that. This was uh a

16:16 in, in elasticity. And uh is uh mostly famous for his works

16:22 quantum mechanics. But here's a big and oh a classical physics which he

16:30 uh was interested in. And when found out that when uh as the

16:36 types of call it, he's got types of breeds. And you

16:39 he's got big ones, the heavy ones and lighter ones. And

16:44 the two types of beads move in that makes an ordinary move of

16:48 you just push on this from from the side. And uh all

16:53 beads are gonna be moving uh uh or less in phase. But when

16:57 move out of phase that makes another and it's slower. And he called

17:03 the optical mode because if these uh two beads are moving in opposite

17:09 that's uh and they're both elect electrically , of course, in his

17:14 And so that's gonna radiate uh um uh hello oh ra el electromagnetic

17:24 And so, uh and that's why called it the optical mode of

17:29 It's the same thing in rocks. , uh what B found in 1941

17:36 when the fluid and solid move in , that's just ordinary, sound like

17:40 been talking about. But when they out of phase, that's the new

17:44 of wave which B invented and realize the same basic physics that Warren found

17:51 13 years earlier. Now, uh We wanna ask where, where does

17:58 show up in our business? not too much because it's mostly uh

18:03 uh at high frequency. Uh If want to excite these things, you

18:08 to have high frequency, this means sonic band or ultrasonic band, probably

18:14 seismic band. Now, these were first detected uh in, in 1980

18:22 a, a Schlumberger guy uh uh . Do you know the legend of

18:27 Clona? No, no, I a very good experimentalist working up in

18:36 . Sure. He's now retired and uh he, after all, uh

18:41 imagine um 39 years between the prediction these waves and the uh detection by

18:50 uh it was a long time for Beau to wait frankly, I don't

18:55 if Beau lived long enough to see happen. Now, they travel very

19:01 100 m per second and with very attenuation uh uh QP not much uh

19:08 uh the, the uh the quality Q is less than one, we'll

19:14 more about the quality factors um in next lecture. So because of this

19:20 uh attenuation, if we, we we were to excite these ways uh

19:27 , uh it would be um oh would be attenuated very quickly because of

19:35 very low quality factor too. Uh we, we maybe need to think

19:46 it anyway. And seismic in the band because you can be sure that

19:50 every sedimentary interface, some ordinary energy converted into this kind of wave happens

19:57 in high frequency than low frequency, maybe some would be converted. And

20:01 constitutes an effective mode of attenuation. , it's gonna affect amplitudes in an

20:10 dependent way, which is not included Standard Avio theory. So we didn't

20:16 this uh uh uh in the core theory that we discussed earlier this

20:22 Uh uh because we assume poor uniform pressure. So the bo wave comes

20:28 the frequencies are high enough, the slow wave comes when you have frequencies

20:33 enough. So that uh the fluid is not uniform, but it's uh

20:39 in uh this part of the pore compared with that kind of a pore

20:44 . Uh We assumed earlier that that happen, we were operating in a

20:48 frequency and, but we cannot, the uh blood pressure does work,

20:57 then you're gonna get real slow waves this velocity, the velocity of those

21:03 is determined in part by the permeability the rock, you gotta be

21:08 How slow depends on the peril of rock. Now, when flow,

21:14 fluid flows locally, that means on grand scale during the passage of the

21:20 , this flow is called fluid It's a bit of a misnomer because

21:25 you think about squirting, you're thinking macroscopic movement of uh fluids uh from

21:33 part to the to the next. uh you know, the amplitudes are

21:37 low. So uh uh it might that we should use a different word

21:42 squirt because the fluid hardly moved at . But uh um when, when

21:47 does move, it has a big . It's responsible for the discrepancy that

21:54 showed earlier between ultrasonic data and gas theory. Remember this slide that I

22:00 earlier. So here is ultrasonic data uh in blue and the theory in

22:07 red showing clear discrepancy. And the for this discrepancy we now understand is

22:14 of fluid squirt at high frequency ultrasonic in the rock. And we uh

22:22 have to worry about these effects and seismic data because our frequencies in seismic

22:29 are so much lower than ultrasonic We uh we still have to worry

22:36 the other thing we talked about this . Now, when fluid flows

22:44 we have not only this stuff which just talked about, but it results

22:49 the attenuation ordinary ways which is lecture , which just happens to be the

22:56 topic. Ok. So let's, , let's move directly into that.

23:02 I will remind you that, after we leave tonight, uh,

23:06 afternoon, late this afternoon, you're go home and you're gonna be thinking

23:11 what we talked about and you're gonna down some questions. And before Friday

23:16 gonna send those questions to me on topic of poor elasticity and also on

23:21 next topic of attenuation. So I going to stop sharing here and I'm

23:31 to um the screen and I'm gonna up the next trial. They gonna

24:16 that in presentation or, and then going to um share that with you

24:31 . OK. Lesson nine A OK. So this is another one

24:38 those topics which is not uh included a standard course and ways and ray

24:44 because uh well, that's classical thinking we need to talk about attenuation because

24:51 see attenuation in our data all the . Have you look at your uh

24:58 workstation? You'll see that the uh coming in at long times have a

25:05 frequency content in the higher uh than than the data coming in at short

25:12 . The reason for that is that those long times uh means long way

25:18 , the uh uh uh the waves lost their high frequency. So now

25:25 gonna talk about that effect and uh what it means for us. So

25:33 see. So by the end of lesson, you'll understand how Books Law

25:41 to be modified to include a generation how this results in a wave equation

25:47 in includes a generation. Of if we're gonna change Hook's Law,

25:51 gonna change the wave equation. We're find plane wave solutions with the

25:57 We're going to uh uh find out that affects reflection. Remember uh uh

26:04 when everything we did with reflectivity um uh I know it generation assumed the

26:11 version of books law. Now here's , an A AAA topic which you

26:17 have um hm not seen coming uh turns out to be related to

26:27 So, dispersion is the fact that depend upon frequency. And now we

26:32 out according to this, that um that is connected to attenuation about

26:40 that's a bit of a surprise. . And then we're gonna talk about

26:47 of a generation. We didn't talk mechanisms of, of elasticity, did

26:54 ? But we're gonna be talking about of a situation. And then we're

27:02 uh uh point out that there's another another issue that we can call a

27:08 attenuation different from real generation, but looks the same if you look at

27:16 uh casually. So that's the program the rest of the, this

27:23 So most of what we are doing uh for the uh first eight lectures

27:29 been classic seismology equally suitable for exploration for investigations of the deep earth.

27:38 now we know that none of it truly suitable for exploration since it ignores

27:43 effects of attenuation. Now we see iteration every day in our data,

27:51 we normally we uh ignore it. we shouldn't ignore it. Maybe they

27:56 teach us something. So that's what says here and we see it every

28:01 . Uh the main way we see is uh because of the loss of

28:07 and also the loss of frequency at recording time. That is all the

28:13 are attenuated, but the high frequencies attenuated more than the low frequency.

28:19 um that's why you progressively uh don't the higher frequencies at longer reflection

28:28 Now, I'm, I wanna point to you that it's a good thing

28:31 uh sounds do die away. Otherwise the sounds ever made on earth would

28:36 be echoing on. Just imagine if didn't have a continuation, then all

28:42 sounds ever made on the surface of earth would be echoing around inside the

28:46 bouncing around. Uh you know, uh uh the stomping on the ground

28:51 the dinosaurs. Uh uh All the, the volcano of everything,

28:56 those noises would still be with us we didn't have attenuation. So uh

29:02 a real good thing that we have . It's also a good thing that

29:06 usually mild attenuation so that um uh can uh observe propagation of waves through

29:16 distances of several kilometers if the a were very high. So the things

29:22 out in 5 m, then we'd out of luck for exploring at least

29:27 the techniques that we typically use So that's what it says. We're

29:33 that the generation is usually weak. , um we're gonna see in

29:41 I think it was the sixth topic the previous list that we always get

29:47 when we have attenuation. But the is also usually weak. And so

29:53 uh we will talk more about that we get to the uh that part

29:58 the election. So first thing we're have to do is we're gonna have

30:03 modify Hook's Law. So this was Law as we uh uh saw it

30:11 we said that uh um with law says that stress is linearly related

30:18 strain. Here's the, the straight and the slope of that line is

30:22 the compliance. And if we did the other way, uh uh the

30:27 of the line would be called the . And let's look at it this

30:30 in terms of compliance. And uh according to hook, the stress

30:38 the strain without delay, you apply stress and immediately without any delay,

30:44 got a strain or vice versa, be the strain causes the stress.

30:51 it be that when you squeeze on rock or on a sponge or on

30:55 balloon. What you're applying is a and what you're feeling, pushing back

31:02 the resultant stress. So, cook not know or care about the answer

31:08 this question, but we should, we should care. We're gonna find

31:15 . Um Well, I, I , you know, already let me

31:18 you the question when you have a of uh seismic energy in the,

31:24 the earth. Uh what happens to energy? Does the energy disappear or

31:29 it get changed to some other So I'm gonna pose that question for

31:35 that. When we attenuate sound in earth, what happens to the energy

31:43 is its uh does the energy disappear does it get changed into some other

31:49 ? It, it, it changes energy doesn't disappear. Yeah, energy

31:54 not disappear. Energy is always And if you think that you've lost

31:58 energy, that means you haven't looked enough and you have to uh um

32:05 out what happened to the energy that think you're lost. And so now

32:09 let me uh um turn to uh le uh if the energy is changing

32:16 some other form, what form is changing into? Say it again?

32:42 . Yes, it's changing into Yeah. So uh uh and so

32:49 soon as she says, thermal, know we're talking about the second law

32:53 thermodynamics and the second law of the says, the, that the entropy

32:59 a closed system always increases. And that means that we uh as

33:04 as the uh when it goes through rock, uh uh uh the energy

33:11 not disappearing, some, most of is propagating, but some of it

33:16 uh increases the entropy in that In other words, it increases the

33:24 . So, um that showed and uh the second law of thermic

33:33 always with us, we never can the second law. It's uh uh

33:39 one of the fundamental laws of And so when Cook is pretending that

33:48 no attenuation here, he's ignoring the law. Uh And he uh he's

33:56 uh so, so that's a dangerous . You don't want to be in

34:00 position of ignoring one of the fundamental of the universe. So this picture

34:08 oversimplified. Of course, when you at real rocks that they behave more

34:13 less like this, when you cyclically uh squeeze a rock, uh think

34:21 uh um any frequency you want in laboratory, imagine squeezing it, un

34:26 it and so on. And when un squeezes, it doesn't come exactly

34:31 the uh uh the same path through strain space, it has uh

34:37 it, there's a little bit of here. This is called a,

34:41 hysteresis loop, a hysteresis loop. so time is going in the direction

34:48 the arrow. So it's going this uh in a diagram like this,

34:53 can't uh you gotta have arrows showing direction of the cycling because um um

35:02 she comes to the wrong conclusion. now here is the point of maximum

35:08 right here. So this is the active here. So this is the

35:12 of maximum stress. But the point maximum strain comes after you go around

35:18 corner and come and coming back. that's later. So what we can

35:27 is the stress always leads the So stress causes strain, not strain

35:34 stress. So I think I posed you the question on the first lecture

35:40 we were talking about hook, the caused strain or strain causes stress and

35:45 forget your answers. But I think was some confusion and I don't blame

35:51 for that because um uh we were uh uh adopting the assumptions of h

36:00 pretty obvious at the time, uh and proportional restraint when you think about

36:05 , that's a violation of the second of thermodynamics. And so uh real

36:11 are more like this. Now, I've shown you here is Annie Lifts

36:14 of course, real walks are not be exactly elliptical like that. But

36:19 was easy for me to draw using um facilities that are given to me

36:26 Mr Bill Gates. OK. we wanna incorporate this into the previous

36:35 that we had developed for hooking So this is what we're gonna

36:40 We're gonna press hook's law just like . This is just like we did

36:45 . But now the stiffness is It's all we have to do and

36:52 . Well, so, uh, , it could be the same with

36:55 . And also, uh, uh, it's amazing to me that

37:01 can recover from the enormous mistake that made earlier. We earlier, we

37:09 one of the basic laws of the . And now turns out that we

37:15 uh uh recover from that with a simple state that when we look at

37:22 Law like this, we got to that um the stiffness coefficients might be

37:34 . So you're gonna see in a how this leads to attenuation.

37:39 So for isotropic rocks, we for example, the in compressibility has

37:43 real part and an imaginary part same the sheer motos, same with the

37:49 modulus real part and an imaginary So now, normally we don't write

37:59 uh the uh you know, the modules in this way. Instead,

38:04 we do is we uh uh factor the real part and that left with

38:11 plus I times this ratio and this , we give the name one over

38:17 and for the P um uh for P wave, it will be a

38:22 of P. And so obviously, , it says uh uh QP is

38:30 one divided by this one, this divided by this. OK. So

38:35 could have defined the inverse of But uh conventionally, we define Q

38:41 the way for you to remember that um uh Q stands for quality.

38:48 the second law says an implication of second law is that Q should

38:53 always be positive. Now, the Q is the less attenuation we have

39:01 the bigger two is means the smaller number is and the uh uh the

39:06 the imaginary part is compared to the part. So um now, in

39:13 of velocity, what we, what we have? So that the P

39:17 velocity squared times the density is exactly to uh the lunch models which we

39:24 learned how to express it in, this way here. And since the

39:29 is real, we can uh uh sort of factor out the density and

39:34 that uh the square of the velocity equal to the real part of the

39:41 of the velocity plus one over Um Because all we did was divide

39:48 density. Now, since the attenuation weak, that is Q is

39:54 then uh we can just take the root of that and that brings in

39:58 factor of, of a half right . And in um the uh

40:05 the solutions for the plane wave, didn't have the velocity to the uh

40:10 the first power we had the, inverse of velocity that's the quantity which

40:16 in uh the phase factor of the wave oscillator factor uh that we showed

40:25 . So we need the inverse of VP. So again, because

40:30 is large, uh we can easily this one from this one simply by

40:36 the sign that's Taylor series arithmetic. we've been doing all this course and

40:42 similar for she. So now uh the difference of course is that this

40:51 is different from the other Q. 22 factors are independent properties of the

40:58 . And normally we have to determine uh from the data and normally they

41:05 uh well, uh OK. I'll it that normally we have to uh

41:11 to uh determine them from the data normally they're in the range of say

41:18 to a couple of 100 four rind , well consolidated, normally pressured

41:25 So if it got, if it's artful gas saturation, it's a different

41:29 . If it's uh a poorly it's a different story. If it's

41:34 over pressure, it's a different But for most of our rocks in

41:38 seismic band, we're gonna be expecting which are a lot bigger than one

41:47 probably bigger than 10 and maybe bigger 30 but less than 500.

41:54 we're gonna be expecting values in this . We will talk later this afternoon

42:03 about this case here and not bri it partially saturated. With gas.

42:09 that's gonna make you substantially lower. that might be interesting point for us

42:15 think about. So. No, consider a sandstone, pure sandstone with

42:24 quartz and um uh Ryan. yeah, for the mineral courts,

42:33 queue is very high. I think in, in intuitively. Uh uh

42:40 think you, you have in mind you have a, um and a

42:45 crystal quartz, uh and you uh it from the side, it's gonna

42:52 for a long time because it's such perfect uh uh material. Uh uh

43:00 uh it's also true that for the , uh if you, if you

43:04 the attenuation of sound um in in a uh in a container of

43:13 this, your brine, it's gonna a large number also 200. But

43:17 , you put these two together and the brine saturated S stuff, you

43:22 get something in between, you get which is lower. So,

43:27 uh the queue for the rock is an average of the queues for its

43:32 . Isn't that interesting? We, uh uh w when we were talking

43:38 um the density, for example, a mineral assemblance, we just found

43:43 was an average of the densities of the inconsistent mines. So uh that

43:53 fails completely when we talk about The reason for that is that the

44:03 factor Q uh uh depends upon a interaction, an interaction between the

44:09 the constituents. Interesting. So, here's AAA question for uh uh for

44:20 and I believe uh Carlos, it's turn. So, uh let's look

44:24 the question and notice down here we none of the above. So,

44:31 oh yeah, I'll, I'll ask uh uh Carlos about part A it

44:38 the theory of elasticity is easily extended a 10 to media by applying it

44:44 rocks. Um Is that true? see. I see. I think

44:55 not, it's not true professor because sort of a non state uh uh

45:03 to rocks. What does that Uh So uh uh uh I like

45:07 answer and I'm turning to Brisa, how about part B the theory of

45:14 is easily extended to a tentative media considering that the elastic businesses are

45:21 Is that tur false? You that is true. Yeah. Uh

45:38 gonna count that as true even though did not yet show you, I

45:42 yet show you how those complex uh is, make it for a

45:48 So that was a tough question because was asking you to uh uh understand

45:54 yourself the implications of what I taught . See, I, I taught

45:59 about the complex distances but I didn't you how that leads to attenuation,

46:06 you figured it out on your own you're understanding, at least you're beginning

46:11 understand you'll understand a lot better um in a few minutes. But that's

46:18 kind of question I like uh a that doesn't uh rely on your memory

46:26 relies on your understanding, your physical and you pass the test good for

46:31 . So turning to Lily, this better be uh uh uh better be

46:36 false because we already found one. true. Uh uh So, uh

46:42 I want you to tell me why false. Says the theory of elasticity

46:47 easily extended to a tenant to medium inside the earth, considering the effects

46:52 high pressure. Yeah. So that's we're expecting that to be not

46:58 But um why? Yeah, nothing said about pressure is implying that you

47:10 uh that makes a situation. So , that, that this is uh

47:15 off the point. OK. So we are going to turn, did

47:20 hear something? Now we're gonna uh uh look at this business here.

47:32 uh Answer B and we're gonna show why Merce it is right there.

47:37 we're gonna do the quasi easic wave . OK. So here's the wave

47:44 as we defined it before. This the vector wave equation uh uh for

47:49 , the uh he wave displacement and got the P wave velocity in

47:55 And if you think about how we this, we never assumed that,

48:00 the P wave velocity is real. can, you can go back to

48:05 uh the second lecture and uh where drive this wave ation at no

48:11 did we ever assume that that's, a real number? So, uh

48:17 uh and now we're just gonna allow to be complex, why not?

48:22 furthermore, look at this, we assumed that it was constant. So

48:27 we're also gonna allow to depend upon . How about that? And you

48:34 verify both of these statements, go to our derivation of the wave

48:39 Uh We never did say whether it's or not, we did never did

48:43 that it's real or not. So now, we're gonna assume that uh

48:46 might be complex with a real part an imaginary part because we didn't assume

48:56 was uh uh real. We can accept that as uh uh a

49:02 still work uh Because of what we before and we don't have to modify

49:08 . So now let's look at some to that. OK. So let's

49:13 a solution as before. But now have a complex velocity which is also

49:19 dependent. So here's our N wave . Yeah, it's got, it's

49:24 uh it's, it's a vector displacement a function of time space and

49:30 And we, it's got an amplitude which is a function of frequency and

49:36 it's got this oscillator factor. And uh you don't see the velocity anywhere

49:42 . What you see is omega T or minus a dot X and a

49:48 uh uh has, can be written this term with three different components of

49:54 uh the uh the wave vector The length of it is given by

49:59 sum of the squares square root of sum of the squares. And that's

50:05 by uh Omega over VP. And we've got VP complex. So uh

50:12 simplicity, let's consider one D propagation positive Z direction. That means that

50:19 can, uh we can drop off , the plus that we only have

50:22 minus here. And uh we got uh uh A minus Z over

50:29 And uh previous slide, we um spelled that out and spell that out

50:36 terms of the real part of Uh And uh uh uh uh and

50:46 over QP, also time is a part of VP. And then um

51:03 separate out uh uh the, the part. So here is um uh

51:13 there's only real numbers in this part the uh of the phase. And

51:19 we still have our imaginary uh uh here. So this part here is

51:24 make for oscillations and this part here have any, I, you see

51:31 we have an I here and I where that and makes minus one.

51:36 there's no I in this part. that's why we separated it out.

51:42 this is frankly only the uh uh well, uh the whole wave is

51:48 in the Z direction. So of , this depends only on Z.

51:54 so this part is gonna lead to , but this part is gonna lead

51:59 attenuation because as Z gets bigger, whole thing gets smaller because of that

52:06 . So, and you can do similar thing for sure. And the

52:11 difference is that you have uh uh some subscripts s here. Now for

52:18 mode, since the velocity is equal the wavelength times the frequency divided by

52:25 pi we can rewrite that in this . So that the wavelength um is

52:31 explicit and I kind of like this now uh this is uh it's um

52:39 I look at something like this, always make sure that we have an

52:43 which has no dimensions to it. here we have uh dimensions of

52:48 And here we uh we also have of length in the denominator Q has

52:54 dimensions I as a no dimension. the whole thing has no dimensions.

52:59 is the attenuation factor only. And what you can see is that uh

53:05 whenever each cycle for each cycle of , uh um it goes one

53:16 So when we, when, when go from uh uh Z to Z

53:20 delta Z sit on, when we from Z to Z plus lambda,

53:27 uh uh that means an, an factor of uh uh of one of

53:37 I'm saying it wrong, sing it . When we propagate over an interval

53:43 Z, the interval is equal to wavelength clamor. Then this uh canceled

53:53 schedule and left with either the IP Q. And so for each such

53:59 , we're gonna reduce the amplitude by factor of uh e to the minus

54:05 over 50 that's about 94%. So a good thing to, to uh

54:13 that as waves are propagating along, are reducing in amplitude by the same

54:22 independent of frequency. See uh when we introduce the wavelength instead of

54:27 frequency, uh uh the frequency is longer explicit. So for every wavelength

54:34 propagation that a wave uh uh it's losing a few percent of its

54:42 um of the tude per cycle. , of course, we have,

54:52 you have higher frequencies that means shorter and you get more cycles per

54:57 So the higher frequencies attenuate more uh a different distance down, for

55:03 from surface down to the target. when we see that the, the

55:09 frequencies are dying out, it's because executing more cycle, not necessarily because

55:18 is less than QP that might be . But the, the main thing

55:23 that um uh shear waves execute Now, uh we specifically said we

55:31 propagating in the distance in the direction positive Z. So if we just

55:37 things out this way we can uh uh propagation in the direction of negative

55:42 . And also we can uh uh in our four year decompositions. We

55:47 include negative frequencies. Uh We put uh uh uh absolute value signs so

55:56 we were getting a AAA decrease in factor because we know that the uh

56:05 two factor are, and the real of velocity are both positive. But

56:10 , it often happens that we're sloppy uh we don't um uh I put

56:17 those um uh we, we don't in those absolute values. So I

56:30 this is uh uh you know, course, how we go about uh

56:33 uh determining velocity and when you determine in those ways that uh you're all

56:40 with, um you're familiar with the that you have been taught to use

56:47 calculating the real part of the So how would we calculate the imaginary

56:52 of the velocity? In other how would we calculate two?

56:57 let's just measure the cube by measuring loss of high frequencies. So,

57:03 we're going to, you know, the same expression that we had before

57:08 this continual fact continuity factor up in , put it inside of a bracket

57:15 call the whole thing uh uh uh amplitude as a function of frequency and

57:23 length of propagation. So this one have any Z in there. This

57:28 does. And then the oscillator factor here. Now we form the ratio

57:34 spectral component we before after and before propagated distant Zs le let me back

57:42 with this uh you take this thing uh so you look at your,

57:46 your data, look at all the uh your work or take a single

57:49 from your workstation and, and look a um um a reflection event and

57:57 a, a window that off and , take a spectrum of that wavelet

58:02 the window surrounding that reflection event. you take another one, you,

58:07 , you look deeper and look at time and find another a a prominent

58:13 further down longer time, make a uh spectrum for that. And then

58:22 uh for a given frequency, then going to uh excuse me, uh

58:28 , the first spectrum that you find this one, the, the,

58:33 second spectrum from a lower reflector is one you form this ratio of these

58:39 spectral components? Excuse me. Professor , I have a question regarding those

58:45 that you said you take the spectral . So does it matter or how

58:52 the, the way that we choose windows affect the result? Like it's

58:58 too a a and so it, never as simple as I just

59:04 you're always gonna say, well, I gonna look at this uh um

59:11 uh is a small uh uh peak between or am I just gonna look

59:15 the big peaks. Uh uh that's matter of judgment. And so you

59:18 it this way and you do it way and see how the answer changes

59:23 respect to your choices of the right? So, uh uh uh

59:30 might be sitting there at your workstation you say to your buddy, I

59:33 , how long should this window And he'll say, oh, it

59:36 be 50 milliseconds. So you try milliseconds. And then because you're a

59:42 person, you also try 60 milliseconds 40 milliseconds and see how things

59:50 Um When you change the arbitrary decision the window L right? You wanna

59:59 an answer that doesn't depend upon the , right? You wanna have an

60:04 that tells you about the medium, about uh uh your uh the way

60:10 handling. So you, you uh it several ways and decide for yourself

60:15 uh for that particular situation. Um uh are the best parameters to

60:23 You folks are always uh adjusting parameters in all of the processing steps.

60:30 furthermore, well, uh let me delay that uh comment. I'll,

60:36 go furthermore, uh in, in couple of seconds, we're gonna

60:40 form this ratio here as a function several different um frequencies. And because

60:48 uh uh this ratio is gonna depend uh uh when you cancel out the

60:55 solitary fact, when you take when look at the spectrum, you're no

61:00 looking at the wavel, you're looking the frequency domain. And you see

61:04 , there's no, um there's no terms here, there's no complex terms

61:10 . Everything here is real. You uh uh you determine the uh the

61:17 of uh depth difference between uh uh reflection and that reflection de determine that

61:25 um uh by uh uh estimating And so, uh by making this

61:33 thing on the right, we're going uh infer what this is. And

61:37 we're gonna infer what is uh um what, what is the cue?

61:48 now if you, if you choose , if you choose two reflectors,

61:54 are too close together, that is you have a small delta Z,

61:59 might get into trouble, uh You'll out for yourself that if you do

62:04 uh for a respected group, uh you do it for um reflection events

62:12 are too close together, you get answers. So you try it,

62:17 uh to say, OK, I forget this uh event, I'm gonna

62:20 down, pick a lower event and you do that, you're gonna lose

62:26 . So, uh uh uh when do this with a large value

62:31 of a delta Z between the upper and the lower reflector, um you're

62:39 uh only determine the average value of in between those two reflectors and the

62:46 the distance, the more heterogeneity is inside that average, the smaller the

62:52 z uh if it gets to be Z is too small, you find

62:57 uh uh nonsense answers. So this is a parameter that you have

63:03 select um uh with your good And um uh so like that,

63:10 having chosen a number of different you know, when you're looking at

63:14 data, you can't choose any delta that you want, you want to

63:18 a delta Z which uh corresponds to different prominent reflectors. And if you

63:24 to, to uh uh do something , you'll find um uh nonsense

63:31 So this is the kind of a that you're normally gonna get as a

63:36 of frequency by whoever drew this cartoon 12345678 frequencies and found that the generally

63:48 ratio is generally decreasing with some scatter the straight line. And uh we

63:56 the best fit straight line through And that, that's the slope of

64:00 uh logarithm of the ratio. And the Q value in by measuring the

64:07 of this line that's measuring how fast losing the high frequencies. Um That's

64:18 uh depend uh that's gonna give you an estimate for the average value of

64:23 inside here. That's an average over interval delta eight and also an average

64:31 the size of band because you're you're calculating the spectra only from here

64:38 here. That's the maximum frequency uh uh that you had available in the

64:43 . This is the minimum frequency. if you were able to look for

64:48 frequencies outside that you might find that uh uh instead of a straight

64:56 you're gonna get a curve. um normally that's not a problem.

65:03 , uh there's lots of scout, it's hard to see curvature in this

65:09 . So you normally, what you're do is settle for a straight line

65:13 and get uh thereby an average value Q over the interval and over the

65:20 of B. So now let me to li li for this quiz

65:30 it says the complex velocity leads to because high squared equals minus one.

65:38 that true? That's false. I the statement is true, but

65:42 it doesn't say anything about complex right? OK. So uh

65:50 it says if you true or high frequencies extenuated more than low frequencies

65:56 Q as the function of frequency is for high frequencies, is that

66:04 I think it's the other way Well, actually, we, we

66:10 say anything about the frequency dependence of . Let's back up here here,

66:15 have a Q which is frequency We have a, a straight line

66:20 . And uh the QE is in formula for that straight line. And

66:25 uh if there is a frequency dependence we can't see it because of the

66:32 . So all we can see is average over this band. And we

66:37 say anything about this uh uh frequency of Q. Maybe it's in there

66:41 maybe not. But I can tell that in most cases, when you're

66:45 at Q, at seismic, at seismic data, normally, you

66:51 do better than to assume um a value for Q independent of frequency.

66:58 band is not wide enough. If had a band that went from five

67:03 to 1000 Hertz, maybe we would some uh frequency dependence of QE.

67:10 normally, we don't have that for the S band, we don't

67:14 it for the sonic band, we have it because in the sonic

67:18 we might have 500 Hertz to 2000 . Again, that's a small uh

67:23 range of frequencies. And so normally can't see a difference in frequency in

67:30 because of the, within the sonic because uh the band is not wide

67:37 . However, it is possible that Sonic cube could be different than the

67:43 cue because those two bands are And furthermore, you can say the

67:48 thing about with the ultrasonic band. normally the ultrasonic band by itself is

67:55 wide enough for you to measure um dependence of Q inside that band.

68:03 when you compare the numbers for the band with the numbers for the Sonic

68:08 and the siding band, you might see differences because, uh,

68:13 uh, now you're looking at cross . Now this is, uh,

68:25 , the next question and, I remember, uh, what Carlos

68:29 is, he said, he thinks the other way around. Uh,

68:32 so this is, uh, the way around. Uh, but I

68:35 the answer is the same that normally don't, uh, well, we

68:41 see uh the sequence depends on And I will be showing you uh

68:48 this afternoon, I will be showing um what we should be expecting for

68:56 frequency dependence of Q. OK. this one, I think I

69:03 I, I just wanted to comment , I am in the slide 22

69:07 have a comment that says that since frequencies execute more cycles per meter,

69:14 higher frequencies at any way more over given distance. Yeah, that's

69:20 But that's because of the uh wavelength uh when you have um uh uh

69:29 when you have high frequencies, they more wavelengths. So the A U

69:34 giving you the amount of um of energy loss per wavelength in that

69:42 So it's not happen because the wave executed more uh cycles, not because

69:50 key is different. So sorry, back to the, to that question

69:56 I was looking at those slides again I was confused. So here,

69:59 trick in the question is that it that that high frequency sin no

70:06 But it's not because it's depend, the Q, it's not of the

70:11 dependency of the frequency. That's It, it's, it,

70:16 it's because those high frequencies attenuate more execute more cycles. Mhm OK,

70:27 . OK. So uh uh um one goes to mesa, it says

70:32 frequencies that generate more than low frequencies higher frequencies execute more cycles than lower

70:39 . Yeah, that's true. That's we just said. Yeah.

70:43 So, uh uh so now we how complex stiffness is make for attenuation

70:52 as it propagates. Now, let's about the next thing what happens to

71:00 when we consider that the rocks are attenuating rocks. So we already said

71:08 when you have a, a sediment gas in it that's highly attenuated.

71:13 shoes are low for these gasses Uh uh uh uh uh uh uh

71:20 , uh you uh you know, that um gas saturated sediments are

71:26 they have slow velocity. That's the part. And now um uh we're

71:31 about the attenuator part, the, imaginary part of the I won't uh

71:41 stiffness. And now I'm telling you partially saturated rocks have high,

71:48 have low Q as well as low parts of the velocity. And the

71:54 is, well, can we use uh can we use that fact to

71:59 for gas saturated rock. Wouldn't that neat? No. Uh The next

72:05 , usually the gas reservoirs are so that there's not too much loss of

72:10 frequency due to a two wave propagation them. That is if it's only

72:16 m thick. Uh So the heat going down the back is only going

72:22 m. If it's vertical propagation, that's not enough to notice the loss

72:27 high frequencies due to the gas in layer. But should we give up

72:39 this line of inquire or should we thinking about it? The next question

72:44 here said, is there effect of on the reflectivity itself? And we're

72:50 about reflection coefficient, not the two propagation. OK. So let's look

72:57 this here is the normal reflection It's a, a jump in um

73:03 impedance divided by uh uh twice the competence, you can separate it out

73:09 the uh uh density part and the part, the density part is

73:14 So we can leave that alone. the velocity part now has uh um

73:20 real part and an immersion part. uh when we say delta VP,

73:25 means uh uh incident subtracted from uh reflecting. And so here is the

73:34 the VP two minus VP one, spelling it out in terms of uh

73:41 real parts and imaginary parts. And down in the denominator, we have

73:46 same thing with the temp plus instead a minus. So uh we introduce

73:52 notation, uh that means we're we're factoring out of each one of

73:57 terms like this, the real part then the imaginary part is proportional to

74:02 1/2, 2. And where the comes from, comes from the fact

74:08 this is uh uh the, the as is the square root of the

74:14 . So that brings in a one factor here. And so here we

74:18 for the uh reflecting body here, have for the instant body here,

74:22 have the sum of them. And what we wanna do is collect the

74:27 and imaginary parts. OK. So there's no approximation here. And um

74:34 , let me see what else I'm do here. Um Oh yeah.

74:45 I'm gonna express this part here as the average and this part here as

74:52 the difference of the real parts. , uh before when we were analyzing

75:01 , we didn't have any terms like one, this one, we,

75:04 didn't have any imaginary terms in in the reflectivity coefficient. And so

75:10 delta V over vs that we saw , in that case, uh uh

75:15 we recognize that it was a, change across the reflecting horizon of the

75:20 part of the velocity. And now an imaginary part that we hadn't been

75:25 about before. Yeah. So I'm to make the approximation here that we

75:39 neglect this part here if the cues large and cues are bigger than uh

75:44 10 or so, we're gonna ignore part, but we'll remember that,

75:50 we ignored it. So, uh we ever come to a case where

75:54 not true, we can go back , and include that part. And

76:01 I wanna do my favorite trick here to use a tailor expansion. We're

76:06 ignore this part. And up we're gonna use a tailor expansion to

76:11 uh separate out the, the linear in uh the cues. So when

76:16 do that, uh we deduce that this approximation that the reflection coefficient is

76:25 we had before with an imaginary which depends upon the jump in cue

76:32 the reflecting horizon. And down here have uh the product of the two

76:39 . OK. So in words, can say that when we look at

76:47 linearized an elastic plan reflection car that's . So uh uh we're gonna do

76:55 a vo problem with linear the plan . We don't want to deal with

77:01 exact not sole equations, we'll deal the linearized expressions. And here we

77:08 did the uh normal incidence term and norm at uh non normal incidence uh

77:14 there's a corresponding um generalization from real complex. I don't want to show

77:20 that here. I just wanna talk this real um uh huh oh The

77:31 this normal ancestor, say it it's got a real part which is

77:36 by uh uh the uh the, real part of the jump pen impedance

77:43 by the real part of the uh impedance. That's pretty much what we

77:48 before. But now we have another on here because now we realize that

77:54 rocks that we're reflecting off of uh are continuing. And so when we

78:03 this complex number in the reflectivity, means that the reflected wave has a

78:10 shape an infinite way. Wow, we had the reflection reflected wavelength looked

78:18 like the incident wavel with a different . And now we see that also

78:23 has a different shape because of this contribution to the normal of reflection of

78:34 . So now how uh how important this? So let's consider a case

78:40 the, the real part is very . So um and I can,

78:46 can, we can ignore this So we're looking at weak reflections here

78:52 we left uh uh uh with uh only this imaginary term. And so

78:59 means that the uh the phase which uh uh that's gonna mean that the

79:06 wavelength off of such an interface is shifted by 90 degrees because the real

79:14 is zero and the imaginary part is zero in this approximation. And um

79:21 uh so now we should ask ourselves OK. So we have a,

79:26 weak reflection, it is phase So if it's, if it's uh

79:32 in at a close to zero it's coming out at close to 90

79:36 phase shifting. Uh And so how is that? Uh Well, so

79:42 just uh do it. Uh put numbers, let's set the, uh

79:48 uh the cue for the uh for incoming uh um uh I for the

79:59 media, I'm gonna set that set to 50. So that will be

80:05 , for example, a Brian saturated might have a Q number about like

80:12 , but it's the top reflection off a, a reservoir. And

80:17 the reservoir is a gas sand for gas sand. That's a much smaller

80:22 . Let's set it here at five five. OK. So then just

80:27 in those numbers, we find that reflection coefficient is 4.5% imaginary.

80:35 that's not such a small number. is not a small number. We've

80:41 , we've set, we assume that thing is really small at less than

80:45 and a lot less than 4.5% So, we might expect to find

80:51 significant reflection coefficient in this situation phase . Wow. Now, in this

81:12 , the incoming wave reflects off the of the gas sand reservoir. So

81:19 never propagated into the reservoir at The reflected P wave did not go

81:24 the uh uh reservoir at all. it did, it would have lost

81:30 frequency here because the, uh, Q factor in the sand is so

81:35 but it doesn't, uh, go , it reflects off the top.

81:39 has not traveled through the itinerary through reservoir, gas res reservoir, which

81:46 highly a generating. So it has lost any high frequencies. Wow.

81:52 this is a way to detect it's a way to detect gas,

82:00 gas reservoirs or even thick gas Yeah, let's, let's never mind

82:06 thinness. It's a way to detect reservoirs independent of anything else that we

82:15 talked about. We haven't done any vo here only uh normal instance,

82:21 . And we found that if uh uh in this limited case where the

82:26 part is negligible, the imaginary part lead to a significant uh phase shift

82:33 the wavel without the loss of any frequencies at all. So when you're

82:42 at data, I always challenge students remember this fact about reflectivity. The

82:51 is gonna be complex because the, rocks are janitors you can count

82:58 And next, next question is uh that gonna be a big enough effect

83:03 uh uh see it in my And we just showed that it might

83:10 , you could put in different numbers find different answers. But it,

83:14 possible that you could find a reflection efficient uh in this scenario in a

83:24 between two layers which don't have a uh we, we, we still

83:34 a, a AAA strong impedance a strong real impedance change. The

83:42 could come from the contrast in cues from the contrast and impedance. If

83:48 true, the reflected wavelength is gonna back up where the uh an appreciable

83:54 and it's not, it's not um , right as it goes back

83:58 it never went through the highly attenuated . So it's gonna go back up

84:04 the same attenuation that came down. it's gonna get back to the,

84:09 to the receivers. And how are gonna notice it? Well, it's

84:12 have a different shape than the Uh uh uh because uh the,

84:19 , the highly attenuating reservoir that it from the top of put in

84:25 um, a phase shifted reflection car . No. So is this a

84:35 way to find gas? Well, answer is in principle, it is

84:39 maybe not in practice. And in , nobody has ever published uh a

84:45 saying I discovered this gas reservoir because the uh uh effects on, on

84:52 attenuation on reflectivity. And we didn't it through a vo or anything like

84:57 . We discovered it through uh the uh the complex nature of the reflection

85:06 and the effects on way, what nobody has ever published a paper like

85:14 . Nonetheless, I think it might true. Why do why do uh

85:19 why do we have uh this ambiguous today? When what I just said

85:32 pretty clear, we should expect phase wavelengths, phase shifted reflections at normal

85:39 off of gas saturated reservoir. We expect that. So how come it's

85:45 a big deal? Well, here's reason why it's not a big

85:48 It's because other things might cause a effect. For example, there might

85:53 an interference from a nearby reflector if have a nearby reflector with a uh

85:58 uh uh normal instance, real car , a real normal instance,

86:05 car fission of opposite sign that's gonna uh uh uh uh give a reflected

86:15 superposition of those two top of those reflections, which is gonna look like

86:20 phase shifted wavel talk that good. uh That would uh look a lot

86:28 the one that we just uh attributed uh uh uh a situation.

86:33 the way you figure you, the you um uh sort that out is

86:38 look at the uh offset dependence and on like that, uh uh you

86:43 keep this in mind while you're looking real data and see that, that

86:47 you uh uh ever see a, puzzling um oh reflection event, maybe

86:55 due to the complex nature of reflection fish, just keep this in

87:02 You might be the one who finds way to use this to discover

87:08 you would be become, I can you that you would instantly become famous

87:13 you discover this, that this effect reliably u utilize to find uh oil

87:21 gas in the subsurface. It would just as much a revolution as the

87:26 of a bo So um if you to be famous, keep this in

87:32 when you're looking at data, maybe uh maybe you'll be the one who

87:36 it. So, uh I think one goes to Lily uh true or

87:48 . It says, of course, ref reflect the reflectivity must be complex

87:53 the density is complex, that's false we never said the density is

87:58 Only the stiffness is complex good for . OK. Carlos. Uh It

88:04 , of course, the reflectivity must complex since the stiffness is complex.

88:09 that true or false? Hm I think it's false. Now,

88:21 didn't we just say uh uh I you might not have been uh following

88:27 the previous discussion. Uh I'm gonna this is true since the stiffness is

88:36 , we've got to have a reflectivity . Uh uh Car uh I'm gonna

88:41 you to go back, go back the lecture material and you'll see where

88:45 uh we showed this explicitly. Yeah. OK. So Rosa,

88:51 one goes here says true or A large Q contrast at a reflecting

88:57 can produce a phase shifting reflection without loss of high frequencies. Is that

89:06 ? Yeah, that's true. And me that's what's remarkable about this whole

89:09 is we didn't lose any high frequencies we never propagated down through the Ainu

89:15 media. Uh, we, uh course lost some high frequencies, uh

89:21 propagating down through the overbred, but was normal. Uh We didn't propagate

89:27 this top reflection. Um Yeah. we didn't propagate through the highly Agenor

89:35 at all. What we got was phase shifted reflection with the same frequency

89:43 . We, we affected the, phase spectrum, not the power,

89:47 the uh not the power spectrum of the outcome upcoming movie.

89:55 So um they actually come to this which I would say this topic is

90:09 uh an obvious one. So what says here is there's a close con

90:18 link between attenuation and dispersion. So let's look at uh this simple

90:26 , we have an impulsive source at r equal zero and it's gonna radiate

90:31 in all directions, unbounded medium. , no attenuation or dispersion or,

90:38 anything. So when you take this of uh uh of the impulse response

90:47 can do a four year decomposition uh uh of this time function, you

90:54 what that is, that's uh zero negative times infinite for uh zero times

91:01 zero for a positive times. So infinite spike at equals zero, that's

91:06 time function. It's a peculiar one it can be for it, but

91:10 just like any time function, it be really decomposed um um into uh

91:18 uh all these different four year And uh an elementary result of uh

91:25 of that four analysis which you learn the first couple weeks of any

91:31 And uh complex analysis is that this here is flat. It is,

91:38 a one. And how did we that's one or we do the inverse

91:42 a uh transform? We, we an integral instead of integral over

91:48 we have an integral over time, same factor here. But it's got

91:52 minus sign. See that minus there's no minus sign up here.

91:57 then what we're doing is we're uh that delta function that spike that I

92:03 before. And if you uh uh not hard to perform this in integral

92:09 that's a one independent of frequency. that's number is a one right

92:15 So let's put this uh uh analysis uh what we previously called the N

92:24 wave equation. So we found a for that previously. And if you

92:30 go back in your natural finance, solution looks like this. It has

92:34 oscillator part where uh uh uh uh but it uh it's oscillating in a

92:42 direction. So we don't, we need a vector product here. We

92:46 have a K times R here and going uh uh uh upwards. Uh

92:51 got a minus sign here. That uh it's going outwards as towards increasing

92:58 and because it's got a uh a term in there, it's got a

93:02 over R out there in front. uh uh then these constants here to

93:07 it uh uh um uh this one has to be the uh since this

93:15 a pressure pulse, this has got be a pressure uh have, have

93:19 dimensions of a pressure also. So uh make it look simpler by writing

93:25 way. So here's the fourier amplitude the source. And uh if the

93:31 is impulsive, that's gonna be AAA number. Uh OK. So now

93:38 gonna do the inverse transform. So is the uh the uh solution as

93:45 function of frequency. We're gonna find , the uh solution as a function

93:49 time by doing the inverse four you . And so all we do is

93:56 uh integrate. Yeah, uh we're integrating over uh the frequency.

94:14 And here's the oil factor. And you walk through this uh interval,

94:21 find that uh uh this, it doesn't take much work at

94:27 Just sort of um casual inspection of . Lets you see that this integral

94:33 , in fact the uh uh the, the fourier recomposes position,

94:39 put all these frequencies together with one here for the uh for the uh

94:46 the amplitude. But a one right that gives you the uh the a

94:50 function uh uh arriving not at T zero, but at T equals minus

94:57 over V. It's just an expanding . So thi this makes perfect

95:05 we can fire off a shot. And the origin and what you get

95:09 is AAA sphere um uh um um expanding sphere expanding outwards, decreasing it

95:19 amplitude with um uh according to one our, our, that's just the

95:24 spreading factor, just an expanding So that's um uh just working

95:31 the uh the source problem for it then backwards and finding that, that

95:36 get an expanding a shell of sound out from our soil. So that

95:45 for elastic medium. So now, assume that the medium has a ation

95:50 it. Let's assume it's a constant for all frequencies. And we're gonna

95:55 no dispersion whatsoever. So, putting assumptions into the solution, we know

96:02 that solution is here is um uh oscillator part. Uh And now we

96:09 an additional part here coming from the that two is uh is not infinite

96:16 is some sort of large um finite . And you see this part here

96:22 no um I in it got an over here, but we don't have

96:27 I here. And so this part gonna be leading to uh uh uh

96:33 decay as, as our increases. part here is gonna lead to smaller

96:39 smaller amplitudes. Now, we do 48 in uh in inverse transform.

96:45 uh uh uh we have an additional out here coming from the Q.

96:50 is what we did before. we have a new term out

96:53 You work through these intervals and you out that the answer looks like

97:00 And in graphical form, it looks this. We do not get an

97:05 uh uh an infinite impulse, we a finite impulse. Uh uh And

97:11 is the arrival time here of, , of this impulse, but it

97:15 forward in time and backwards in So this shows that the energy begins

97:22 arrive well before the arrival time. this is uh uh the arrival time

97:28 call it is the radius divided by real part of the velocity.

97:33 so see uh uh uh the energy arriving way here early and then most

97:42 it is arriving when you expect, some of it is arriving early.

97:46 does not make sense. Furthermore, symmetrical about the peak, it's zero

97:53 . Whereas real data is always close minimum phase. So these two conclusions

98:02 that analysis do not make physical So what we conclude is that one

98:08 the assumptions that we made to get this point must be invalid.

98:14 we only assume a linear wave equation a uniform medium but in homogeneous

98:21 right, it has a source turn , I mean by N homogeneous

98:26 Uh But it's linear and I don't anybody wants to argue with a linear

98:32 equation. We assume constant Q. that one is problematic. Maybe Q

98:38 not really constant. That's what, what we assume and we assume constant

98:45 . And so maybe that's not In fact, both of the,

98:51 these live are, are physically I think nobody is surprised if I

98:57 you that for real walk Q is be dependent on frequency. And you're

99:05 aware that for real walks the velocity depends upon frequency. So um uh

99:14 that these are implausible, we can back to the lab and uh do

99:18 lot more laboratory experiments and they, we're gonna prove then that both of

99:22 are mathematically impossible. So the conclusion , yeah, because of the second

99:29 , Q must be finite, we're have a generation, it's not gonna

99:36 infinite anywhere. I mean, Q not gonna be infinite anywhere. That

99:41 that a generation is not gonna be anywhere we're gonna have Q. And

99:47 , uh uh uh uh we conclude she gotta be frequency dependent and velocity

99:54 also be frequency dependent. So that uh uh came as a big surprise

100:02 geophysicist. And it was, it about the time that I was coming

100:07 this business and I was new And , so I was not so surprised

100:12 others, but people who had been it a long time said,

100:16 I've never thought about that before. , um, you know, that's

100:21 . So you can s see more this discussion and the textbook by Akan

100:27 , which we talked about before. even though this is, um,

100:34 uh these statements are generally true, also true that when we have a

100:39 bandwidth, it's a convenient uh a convenient approximation to assume that both

100:48 and velocity are constant within our sighted . And we can make that operationally

100:57 assumption. Even though we know that we look at other bands, we

101:01 find other values for Q and other guys from velocity. That's what it

101:09 here. When we consider other we can see it, but we're

101:13 not gonna be able to see any defensive Q or um velocity within the

101:21 B. No, we're gonna shortly up the issue of mechanisms of a

101:34 . It's gonna be a messy But for a wide variety of different

101:41 , this the relationship that we just concluded between attenuation and dispersion can be

101:48 in this way if you have a different frequencies. Um uh yeah,

101:57 uh sentence one and two take the and the ratio differs from one by

102:05 which depends upon one over Q and logarithm of the ratio of those two

102:12 . And so if you want to about um uh angular frequencies, uh

102:17 the same and uh uh uh uh frequencies is the same. You,

102:22 uh uh uh that doesn't make any difference we're talking about uh uh

102:30 cha this is the sort of change can expect to see um for almost

102:40 mechanism. So we're gonna um we're to postpone the issue of mechanism c

102:47 of mechanisms for some uh uh uh bit of time. And we're going

102:53 uh use this expression to understand what expect to find in our data.

103:02 just put in some numbers say, the frequency one at the top of

103:05 band and frequency two at the bottom the band and take two equals

103:10 And so you ask yourself uh uh is the, the ratio of velocities

103:16 the top of the band compared to bottom of the band? Because we

103:20 50 here in the denominator and we the logarithm of uh 10 here.

103:27 the ratio of the frequencies that comes to be 1.5%. So, uh

103:33 you might not think that uh you uh detect differences in velocity between uh

103:41 top of the band and the bottom the band at the level of

103:46 If it were 10% you might be to see it. But uh

103:51 Probably not. So, so this why we normally don't um uh worry

103:58 , oh uh no, non cons and about uh velocity dispersion when we're

104:09 restricting our attention to a single band data, like the seismic band,

104:16 . Not much you put in your number, you'll find that their small

104:21 . However, if you have partially sediments, same bandwidth put in here

104:27 instead of two equals five instead of equals 50. And that puts a

104:32 here instead of a 50 suddenly we're about 15%. So you might be

104:38 to see that. Now, it be that um uh this partially century

104:47 is so thin that you can't estimate uh velocity such across such a thin

104:53 accurately. That's a different issue. the resolution issue that we talked about

104:58 . Uh But you, you see uh uh this is a powerful

105:02 You can find that you can figure the sort of, of, of

105:08 you can expect to find over a B by this simple formula. But

105:16 let's compare different bands. So here have top up, the F one

105:20 uh in the sonic band 2000 Hertz one is in the middle of the

105:24 band. So we formed this ratio , we're gonna be applying this for

105:29 saturated segments. So we got a down here and now we have 22%

105:36 . So that means that you should to find differences in se velocity compared

105:44 sonic velocity of the order of 20% of dispersion. OK. Now,

105:54 what we said earlier about, Um uh about the difference between uh

106:06 Backus average velocities and ray theory Remember we said that that difference is

106:11 to friendly multiples and that if uh if you have lots of thin

106:17 which we always have that we're expecting uh low velocities because of the friendly

106:27 effect, that's entirely separate from this in the front line multiple effect.

106:34 didn't have any energy being converted to . We only had um arrivals superposing

106:43 itself which were delayed by various delays the various thin beds, which you

106:49 see in the sonic bands data, can't see them in the uh in

106:54 seismic data. But that makes an attenuation which is also uh uh in

107:02 uh uh uh uh it also is to uh decrease the velocities. The

107:10 uh band velocities are gonna be lower the sonic band velocities and it can

107:17 a similar amount. So uh uh expect to find exact matches between um

107:26 sonic and seismic velocities because of these separate effects and some of the other

107:33 that we talked about before. Just remind you the, the sic velocities

107:38 measuring um velocity near to the the C velocities are measuring velocities way

107:46 from the ball hole, you a, a kilometer, two

107:49 10 kilometers away from the borehole. so uh it might be that there's

107:55 and homogeneity going up. And so effect can also mess up the comparison

108:01 sonic and se velocities. So, and she uh uh things I can

108:10 really complicated uh because of uh of effects which we ignored with our first

108:18 through the theory. Now, this here, this is the difference between

108:24 velocity that at one frequency and a at. Well, another frequency let's

108:31 that the, the two frequencies are together. In that case, uh

108:36 uh these um uh uh in that , this ratio looks like a

108:52 So we can rearrange this formula so we see that the inverse of cube

108:57 proportional to the derivative of velocity with to frequency, see how we do

109:04 . So the um this number is be about one, I'm saying I

109:27 . If we assume these two frequencies similar, then this algorithm it can

109:34 uh we can express uh this in of uh uh uh uh v of

109:41 one equals V of omega two plus V. And this one omega one

109:48 one equals omega two equals omega one delta omega. Um And uh rearrange

109:56 in terms of differentials and this term into this. It's the logarithmic

110:04 The inverse of Q is proportional to logarithmic duties, which means uh the

110:10 dimensionalized derivative here, we have the omega and here we have omega over

110:15 taken together those uh things are called logarithmic derivative of V with respect to

110:23 multiplied by pi. So that's the form of this relationship here. Just

110:31 losing Taylor's air. So let's do little quiz here. And uh uh

110:39 it says here, true or I think uh I think Rocky,

110:43 one is uh uh Carlos. This uh coming to you. It says

110:48 the generation goes along with dispersion, follows that a higher value of Q

110:55 a larger dispersion. And now to that, let's just go back

111:01 Yeah, larger value of two means is smaller. And so this ratio

111:07 gonna be closer to one. So means less to this person.

111:13 OK. So this brings us two . But yeah, can we

111:20 can we back up three slides is I have a question about that when

111:25 mentioned that, you know, another , please. Mhm I think this

111:32 one. Yeah, when you yeah, that one. Yeah,

111:37 said that partially saturated segments. It that the segments are saturated with gas

111:44 this, this number of qeq equals five. It's, yeah,

111:50 that's typical for gas saturator that I I, that's not quite clear is

111:56 so thanks for pointing out? Um uh So here a mechanism of

112:05 . So this is a good place stop. Uh uh For a

112:08 Let's take a quick break here and back in 15 minutes to talk about

112:14 of attenuation. Now you're talking Yes, let's get started again.

112:20 up where we left off with mechanisms attenuation. Um uh You know,

112:26 don't know whether everybody is back but let's assume that everybody is back

112:33 we'll continue with mechanisms of a I think it's really interesting that we

112:39 talked about mechanisms of elasticity. Did uh uh talked about elasticity and,

112:45 how uh you apply a stress and get a strand, all that.

112:49 We never did mention anything like uh tiny particles, uh atomic particles with

112:56 between them. Uh We didn't say anything about the medium except what was

113:04 implicit in Oaks law uh proportional restrain since it's a ainu of we have

113:14 complex coefficients. But in there, no discussion of uh what's going on

113:21 the microscopic level. That's kind of . But now we have to say

113:27 it is that is causing this No Hook's law assumed perfect elasticity.

113:38 we want to do better. So gonna write down here below a more

113:44 consti more general constitutive equation, all equations for the standard linear solid.

113:53 you can see that it's got in derivative terms and it's got uh

113:58 a derivative of the stress with respect time and a derivative of the strain

114:04 the spec to time. And uh we ignore those terms, we get

114:11 Law. But this is a generalization H's law. Uh When we say

114:16 , a constitutive equation means we're talking equations which uh which uh describe the

114:25 of the mature. And so this one, this is more general than

114:30 uh that hook uh uh would have hook would have looked at this and

114:36 , I don't want those terms in , get rid of those terms.

114:40 in the 20th century people realized that apply the ideas of elasticity and poor

114:49 to more complicated systems, we have have a more general law. And

114:55 uh um this is uh uh uh express in which is, is sufficiently

115:08 to include and to include many different uh mechanisms. And if you look

115:18 here, we have a derivative of with respect of stress with respect of

115:23 and also some kind of a parameter . And so that uh is a

115:29 , it's gonna have the dimensions of , right? Because uh uh uh

115:34 term here has to have the same as this term here, gotta have

115:39 dimensions of stress. And so then gonna put in here uh uh the

115:44 sort of thing on the strain we're gonna have a characteristic relax relaxation

115:49 for strain at hooks law is just special case. Now, what does

115:56 look like in um uh in, graphical form? If we are

116:05 Uh Yeah, if we draw the , you'll, you'll recognize this as

116:13 modulus, the square of the uh velocity. And this has is uh

116:19 uh as a function of frequency down . And so this is uh um

116:27 the phase velocity in the relaxed So I'll tell you what that means

116:33 a second here is the un relaxed . And you see the curve shows

116:37 its functional frequency starts low and it goes through an increase and ends

116:43 high and flat. And uh before look at the other curve, let's

116:49 again, uh let's look more closely the uh at the frequency uh scale

116:56 here. This is non dimensional frequency the actual frequency is scaled by the

117:03 root of these two times. So see these are the same two times

117:07 define on the previous slide. And uh this uh square root has the

117:13 has the physical dimensions of time. so uh that time uh uh multiplied

117:21 omega means this is a non dimensional of frequency here. And so when

117:26 frequency is one, that's where the is increasing and not only at one

117:33 say between uh uh 0.1 and Uh uh that's where the uh

117:39 the velocity begins to uh deviate from the limiting ace here and heads up

117:48 ramp and then ends up flattening out 10 here. And as the same

117:54 the uh we call that upper upper uh condition, the un relaxed

118:04 . And we'll explain later what that . Uh uh Let's look at this

118:09 curve here. That is the attenuation a function of frequency and that peaks

118:17 uh non dimensional frequency equals one. it's uh it's a definitely a function

118:24 frequency, right? We said we, we wouldn't be surprised if

118:29 frequency, if the Q dependent on uh here is uh here's the case

118:36 that's explicit frequency dependent Q. And is this uh uh um maximum

118:43 of inverse QE maximum attenuation minimum of Q is given by this expression

118:49 where these TS are the same as uh uh T you saw. Uh

118:58 . And um so now let's uh say that uh uh put on here

119:04 geophysical numbers. We, we generally that the seismic band is down here

119:08 the relaxed regime. The ultrasonic band up here in the unreason band is

119:16 things are changing. And so this here is a good characterization I think

119:26 sedimentary rocks, if you were uh a geophysicist, you might choose some

119:33 , a way to characterize this. the uh uh the extenuation reaches the

119:43 , it's max in the middle of band where the, where the velocity

119:48 changing. So now let's talk about what we mean by relaxed and

119:54 So a good way to think about is um when we're talking about

120:00 the principal mechanism for a generation of wave is that uh the movement of

120:09 at the grain scale as the P is going through. So as the

120:14 wave is going through, consider any uh any mass element, say the

120:24 of your fist uh which has lots grains in it. Uh But it's

120:30 compared to the wavelength of the cell . So I, so we can

120:36 of the uh the uh the fluid in there to be uh uniform,

120:44 at high frequency, it's not Uh uh If we uh uh uh

120:51 think now not about seismic waves but ultrasonic waves going through that same

120:58 the wavelengths are much shorter wavelengths are smaller than your fist. Now,

121:04 wavelengths are of the order of a or so. Uh And but it

121:10 comprises many grains. And in uh that uh if you have a mass

121:20 uh the size of a centimeter and many grains and a complicated core

121:27 And as the wave goes through, configure, it's going to be squeezing

121:34 water from one part of the pore . To another because the uh the

121:40 parts of the pore space which are , don't deform as much as the

121:47 of the pore space which are like . So the crack like portion of

121:51 pore space is gonna be squeezed in the wave goes through, squeezing water

121:57 the equate pore space. And those idealize descriptions. But I think you

122:02 the idea because the pore space has complicated shape, fluid is gonna move

122:09 and as it does, so it's have viscous losses inside the fluid that's

122:15 convert the um uh elastic energy to uh like uh the was describing.

122:23 so there's a, there's a, attenuation loss and while it, while

122:33 doing that, uh uh I should that uh uh it, the fluid

122:40 to do that a as the uh the fluid is uh as the wave

122:44 passing through. If the frequency is high, it's not enough time for

122:49 fluid to move from one part of rock to another one. So we

122:53 that the a band, the, um the mechanism. In this

123:04 food court hasn't had time to happen the wave has come and gone.

123:09 that's un relaxed at the other end the spectrum down here uh at the

123:13 uh the lower end of the that's when uh those frequencies are so

123:19 that uh uh pressure equalization has already . Uh while the wave is going

123:27 and at intermediate frequencies that's partially So at low, at low enough

123:37 , we say that mechanism of attenuation relaxed. So you think,

123:47 this is pretty easy. Uh I'm to uh do a bunch of experiments

123:51 a bunch of frequencies and phrase and determine these curves. And therefore,

123:57 gonna determine what are the values of sub uh tau and T sub

124:04 Well, it turns out to be more complicated than that because real walks

124:09 have a single dispersive mechanism like but they have several. And so

124:15 you have several dispersion mechanisms, all the same total, um all with

124:21 same uh maximum uh for a uh it, it's gonna look like

124:28 like this. You can see that is a superposition of many bell shaped

124:32 like that. And when you superpose this uh or uh uh several curves

124:41 this, and you're also gonna have superpose several curves like this. And

124:45 you get then an extended ramp and velocity and uh uh you only get

124:52 have constant velocity when you've uh got frequencies which are higher than uh the

125:02 the maximum frequency for the uh the which is represented by this little tique

125:11 . And you can say, say same thing about, you know,

125:16 uh about the low frequency end. for rocks, um the general conclusion

125:23 , let, let me say here , at so frequencies, the dominant

125:29 mechanism of um of uh attenuation comes the scattering of waves off of the

125:39 grains. The fluid has not had to move uh uh at ultrasonic

125:47 So you can consider the ultrasonic frequencies the limit, they don't move at

125:52 . And so the pore pressure is be different within the cracks and within

125:57 pores as the food has not had to pressure equalize. But still,

126:03 will be uh scattering off of the high frequency. So that's the mechanism

126:13 scattering uh uh uh uh in the band in the sonic and the seismic

126:22 . The uh the dominant mechanism of is uh fluid squirt. Let me

126:33 that again. I I said it in the sonic band. The dominant

126:38 mechanism is smooth squirt. And when get down to the sonic band,

126:43 we hope is that smooth squirt has happened. It's all gone away and

126:47 have uniform pore pressure at uh uh the grain scale uh uh for uh

126:55 seismic waves. Now, it's still frequencies. You can imagine that the

127:00 might flow over long distances, more a grain scale, it might flow

127:06 uh uh the, the distances uh to uh a wavelength. So that's

127:13 long way for the fluid to, flow. And so that means

127:18 very low frequencies. We don't see sorts of low frequencies in our kinds

127:23 data. So um oh I should uh uh at um I should say

127:44 when we have sonic band frequencies, frequencies have long wavelengths. So whenever

127:55 have a long wavelength that of propagating the earth, it's always almost always

128:01 to encounter heterogeneity, heterogeneity of uh rock types. Like if you have

128:10 wave, which is 100 m you know, there's gonna be lots

128:14 different Ortho in uh uh contained within cycle of such a wave. So

128:22 that case, it's a bit um really difficult to say. Um what

128:33 is the mechanism of attenuation at sonic because of the heterogeneity, which is

128:44 gonna be affecting those waves also. um the uh there is gonna be

128:50 effect of apparent attenuation also in seismic uh band, we're gonna have this

128:59 of friendly multiples that we talked about seriously. So that's why this statement

129:05 extends the fluid start mechanism down into uh into the uh seismic band.

129:13 It's uh it, it's a bit whether fluid skirt stops it at 10

129:19 or maybe it stops at 100 Hertz somehow, but it's definitely present in

129:24 sonic band. So uh at ultrasonic , this is what I said um

129:32 scattering. Now, let's think about seismic b, the strength of the

129:40 squirt me mechanism depends upon what kind food we're talking about. If it's

129:47 with oil, then the squirt mechanism enhanced because of the greater viscosity of

129:52 compared to Bryan. Furthermore, if pore space is partly saturated with

130:00 then the squirt mechanism is greatly enhanced if you have gas in the pore

130:06 , basically that move, that gives the uh fluid room to move

130:12 And so uh that mechanism is greatly . So in these situations, uh

130:17 the, the, the quality factor P waves is low, but for

130:24 waves, it's normal. So I'm take you back to um uh to

130:30 hall and to show you what a state of the art image looked like

130:36 Val Hall in 1994. So this uh about the time that I was

130:45 the oil industry uh previous to I had been a professor in the

130:49 University of New York and I joined Amaco, I think in 1995.

130:58 I didn't know these guys. Uh I didn't know at that time,

131:01 didn't know this guy uh barred since . He's become a very good

131:05 And by the way, he no works for VP. Uh But at

131:09 time, he was working for Amao Norway and he was in charge of

131:15 this field here. And it's in southern part of the Norwegian North

131:22 It's near where the boundaries, the of Norway and Denmark and the United

131:33 in the middle of the North This is not too far from that

131:38 . And so their first images look this. You can see uh uh

131:45 fidelity of this image. You have confidence what's in here. Uh Because

131:57 is a time section, not um not a uh depth section. We

132:03 didn't know how to do depth migration that time. I'd say, I

132:07 say we were learning how to do migration. It was very common to

132:11 at um sections in uh as, time migrated sections. And so the

132:19 uh uh is this uh structural depth is this caused by a late arrival

132:27 because of low velocity up here? an open question when you look at

132:34 and you can see up here that uh uh reflection seem to be going

132:39 across, there's no stroke, there's depression here. And you see this

132:44 1400 milliseconds down. So uh the uh reservoir is gonna turn out

132:51 be down here at 2800. So is halfway down. And above

132:56 the ex the data quality was What kind of, of um imaging

133:01 this? This was done with the of the art at that time,

133:06 was called ad mo stack dip move stack. And it's uh uh a

133:13 of what we showed earlier about uh stacking traces with uh uh dick move

133:21 . So all those early ideas were enough for simple um substructure. But

133:28 we explored more and more, we more and more complicated um uh structures

133:34 the subsurface. And this is an . And so we're gonna need better

133:39 algorithm than DMO to see this. you know, we're not gonna wait

133:45 that happens. We're gonna drill it . So when, when Amaco was

133:51 that, I think this was uh in the two D era. I

133:54 it was two D data. And Amao uh was uh uh doing that

134:09 expression in the early part of the nineties, got this image in

134:16 And you know, some manager is say we're gonna drill that sucker,

134:21 they had no idea what was down . Some kind of a monster down

134:25 that it might be, you severely overpressure. Nobody knew what they

134:30 gonna get. Some adventuresome interpreter had a line through there and some lines

134:37 there. But nobody was willing to whether that was a structural depth or

134:43 , a velocity push down. Nobody . So you can uh try to

134:50 yourself as the uh as the staff who's going out on the drill

134:57 And so uh uh uh uh before leave the office, uh uh

135:02 you look over this image and the says, uh OK, Ralph,

135:07 gonna be our ge expert on the ship when, when we drill into

135:13 . So, the first thing you do is you wanna go back home

135:16 make sure your medical insurance is fully up and you, you kiss your

135:22 and kids tenderly goodbye and you say , I hope to see you in

135:27 weeks time and then you uh uh then you go up out and

135:32 So that's what they did. And discover that all of this was velocity

135:38 out and um uh attenuation due to uh you know, uh uh just

135:48 to, to Hewa attenuation here. the structures were really structures that looked

135:54 this across here, mild structure, like this. So uh they came

136:02 uh and they said, honey, home and uh uh um uh it

136:06 a million barrel AAA billion barrel oil . Great discovery by Amako in the

136:14 , early nineties based on this kind data. So the top of the

136:20 here is here. And by by the way off to the

136:23 uh it's good imaging simple structures and only the crest of the energy is

136:29 well, OK, it's messed Well, we have come to realize

136:33 the uh well, uh that this data quality comes from what we call

136:39 ga a cloud of gas. this is the peak, the,

136:43 actual reservoir looks like this sort of , a gentle dome. This is

136:48 crest of the dome right in You can't see it on this image

136:54 at the crest gasses leaked up out the crest of the dome over G

136:59 time collecting here and the overburden, over here, but only here and

137:04 the crest and not up here. the uh uh the reflectors go through

137:10 fine up here. So the starting here on down the um the gas

137:16 collected in there and it does two . It uh makes the wave,

137:20 P waves slow down and, and attenuates them. So it says slows

137:27 waves down and it's and retain. uh this is not the only place

137:32 this in the North Sea, there's be 20 places like this and one

137:39 them is owned by um mhm Before tell you that part of the

137:52 I'll, I'll, I'll complete that one of them nearby about 20 miles

137:56 was owned by conical and it looked similar and it was also a major

138:01 reservoir hiding beneath a cloud of gas that. So, um this is

138:07 picture of what we uh uh uh was happening. We had P waves

138:13 down uh and outside the gas cloud up through the gas cloud, uh

138:19 gas cloud. So it slows the uh the P wave down and it

138:23 lowers the quality. So the waves not only slowed down but they're

138:30 So, um, over, um, uh uh uh I said

138:47 miles away, there's a, a field, um, owned by Cono

138:51 called Eco Fsk. Um Norwegian by the way, Val is a

138:58 name, uh, that refers to , the, the, the home

139:03 the gods. So that when gods die in battle, they are taken

139:09 , uh, what we might call . It's called Val. And uh

139:13 that's uh in Norse mythology. That's it's called. Uh Sometime you might

139:23 to go to the Val Hall bar the rice campus. It's a,

139:28 hidden bar. Nobody knows where it . It has uh uh no signs

139:34 indicate where it is. But if find the right student say, take

139:40 to the Val Hall bar, you see that. OK. So,

139:46 Ko had a well, like had field like this, but they didn't

139:50 what to do with it either, maybe 50 miles away. State Oil

139:56 a field. So sta oil is Norwegian State Oil Company. And uh

140:02 they had the bright idea that let's um explore for, oh, uh

140:15 see if we can get a better of our reservoir that looked pretty much

140:21 this one over and they, they uh uh this field was called uh

140:29 hm So on the tip of my , but I can't bring it off

140:40 tip of my toelle, that's what field is called. Toleen. Another

140:45 name, don't know what, what means, but their field looks a

140:49 like this. They had the same of image data quality, but they

140:53 the bright idea let us explore for with sheer waves. Why? It's

141:01 if this is a shear wave going instead of a P wave, it's

141:04 compressing the uh rocks as it goes , it's shearing the rocks sideways.

141:10 gonna come down here and reflect and back. And so when it goes

141:13 the cloud of gas, it's gonna shearing the cloud back and forth,

141:18 compressing it in any way. So sheer wave doesn't care what is the

141:24 in the pore space, just like talked on the previous uh lecture lecture

141:30 eight about uh or elasticity. We that uh uh the nature of the

141:38 has a very minimal effect on sheer . It affects the density, but

141:42 does not affect the shear modulus. let's, they said, let's do

141:48 wave exploration and get a good shear image of our field at Toma

141:55 Well, the problem is that is there is a uh a ocean uh

142:01 there and you can't send sheer way the ocean. So they said,

142:07 , let's invent ocean bottom seismic receivers put the receivers on the sea

142:15 then we'll have our sources generating um P waves in the ocean water,

142:23 like before both P waves will go through the water, they'll convert to

142:28 waves at the sea floor. They'll down into the sea floor and uh

142:35 uh uh up through, uh down the uh c far down to the

142:41 , reflect, come back through the cloud, unaffected and received by our

142:46 component receivers easy. And so they it and they had success and they

142:58 their results at a meeting of the Ex uh Geophysical Society. And they

143:04 great applause for that successful imaging at through the gas cloud. What she

143:13 , I was not there, but colleague was there and he came back

143:17 told me all about it. He that's great stuff. And so shortly

143:22 our uh our friends at the chemical in Norway said, hey, let's

143:28 that. And uh and so um they said, will you help us

143:34 that? Well, we, we know how to do uh ocean bottom

143:38 uh side. So we, we , I was working in Houston at

143:42 at that time, you know, AMMO Corp office is out in West

143:47 uh about 10 miles from where we here. So we agreed to

143:53 we agreed to help the Amer Norway acquire and process um uh energy you're

144:02 uh a converted cheer. And so the process of that, we learned

144:08 important we learn that when the heat goes down through the water and hits

144:15 sea floor, most of it doesn't to, to um shear waves at

144:21 point. Most of it can uh as a P wave, the shear

144:25 only like 5% something like like it's converted to uh shear waves.

144:32 so by the time those shear waves back uh from the deep reflector,

144:38 they're all dissipated, they, they just too low amplitude to uh see

144:43 , you know, they, they uh they, they um a as

144:49 propagate, they lose energy faster than P wave because they're making more cycles

144:54 meter than the sea than the P are. And so they uh attenuate

145:01 more in ordinary rocks like we see , ordinary rocks, they attenuate more

145:06 P waves. So when those sheer came back, we detected them,

145:11 very, very weak it, but were lucky. Uh the, the

145:17 uh helped us in another way because found out that if we ignore that

145:22 and if we consider P waves which convert of the sea floor, but

145:27 all the way down to the reflector convert to S waves right here coming

145:34 through the uh uh red to to the gas cloud as a sheer

145:39 converted at the reflector that one we see. So here's the, here's

145:44 diagram of that. So here's our wave coming down. You see,

145:47 conversion point is different than the We talked about that a few lectures

145:53 , coming up as a sheer wave transversely to the propagation direction, it

145:59 through the sheer wave and it doesn't care what fluid is there, uh

146:04 or oil or gas or what it's continue. So it comes up

146:08 And so that makes uh uh uh wave arrival And it, it,

146:13 all the possibilities for converted wave conversion . This is what we call the

146:19 wave. That's the one which which converts upon reflection. That's this

146:25 here. So we made images at hall using that idea. And here

146:30 the images, this one came This was our first image, I

146:37 say it was very naive imaging, at least it recognized the uh uh

146:46 of converted waves. It's coming back than the uh uh P wave

146:52 See, let's see, this is , oh this is converted to P

146:56 times. So the, this is a, a time section. Uh

147:00 is uh 3000 milliseconds. Uh So converted to P wave time. So

147:06 looks normal to us. And sure , we see a low dome across

147:11 . And so maybe there's some issues uh uh with better imaging techniques developed

147:20 then, since 1996 we see that absolutely continuous across here. So all

147:28 did was normal move out processing. was in charge of the team that

147:33 this processing and I am not a uh I am not an expert in

147:40 wave imaging. That's why they, don't have me to be teaching the

147:44 in imaging. They have Professor he is a much more expert than

147:49 in imaging. All I knew how do was remove the move out swimming

147:55 layers of these are pretty much flat um I stacked them up according to

148:02 . That's all I knew. But recognized that the image point was uh

148:08 over here instead of here at the . And of course, we recognize

148:13 the uh wave coming up is gonna coming up as a sheer wave.

148:21 um you can see here that the this overburden here, here's to where

148:28 gas cloud is. Nothing really unusual , not very good imaging, but

148:33 not a disaster like we showed Here is the uh the top of

148:37 gas clouds coming in here. And uh um below, this is the

148:44 cloud, it's bad but not I would say. And so all

148:49 these um um um our boreholes posted this image by uh barred after we

148:58 the image. And he said these four of the boreholes that we have

149:04 they uh they validate exactly what this is. And so when we showed

149:09 image, uh um uh to uh went over there to Norway, I

149:16 I went alone, went over there Norway and uh um uh convened every

149:22 together to look at our um uh novel processing of converted wave data from

149:31 . And to show that to the people in Norway who were trying to

149:36 a living by producing here, they never seen an image like this

149:42 What they saw was a terrible image I showed you five slides ago.

149:48 I can tell you they stood and when they, when they saw this

149:52 because now they knew they could manage field properly knowing how it was shaped

149:58 uh in true space in, in depth. I'll tell you another uh

150:07 incident that comes from this. Uh told you about um I think I

150:13 you about this um of the image the uh the incident where the,

150:19 young man who had processed uh the data for our acquisition contractor. He

150:26 it got different results because he didn't the intricacies of converted wave propagation.

150:33 I talked about that some time uh several elections ago. Uh Now

150:39 wanna tell you more uh one more incident there. Uh A friend of

150:46 in Houston took the same data set he uh worked with it and worked

150:53 it and worked with it. He to me a year later and he

150:57 Leon, I've tried to uh uh this data set oh with 80 different

151:06 of velocity, 80 different velocity None of them showed any uh uh

151:13 like this one did. And you got this one on your first

151:18 So that all that has made a impression on me. He was thinking

151:23 this data set was pretty much like data sets that he had seen.

151:28 just gonna tweak his procedures and he's solve the uh flattening problem. And

151:35 , you know, he's gonna uh low velocities, uh you know,

151:39 the crest of the reservoir, it work and it didn't work and it

151:43 work, it didn't work. And 80 tries he gave up, we

151:49 success on our first try because it a new type of data. It

151:53 converted web date, it's not web date. It was converted web

151:58 with many different issues associated with And we, we discovered those and

152:04 dealt with them and we came up this image on our first trial.

152:09 you should do that too. Uh you look at uh uh a new

152:15 , but with the old type of , uh you're gonna be able to

152:19 that problem by tweaking the old But if you're looking, if you

152:24 a new type of data that you seen before, you should look at

152:30 with open eyes and figure this thing offer issues which were not even part

152:38 my previous thinking because it's a new data. So that's what we did

152:43 1996 and had this success. by the way, when we have

152:52 success of the, just the fact we have a successful image here tells

152:58 that the problem is attenuation in the . It tells us that the reason

153:06 didn't get um the heat wave images because we had no P wave energy

153:13 after it was attenuated by the cloud gas. Before, if you have

153:21 poor imaging with P and good imaging convert with C waves, that essentially

153:27 that the problem with the P waves coming from uh the gas CL.

153:33 so now that says, OK, we know there's lots of gas

153:36 Uh Is this maybe um commercial area we attend in it, intentionally tried

153:44 reduce the gas here? Um Would be worth our while? We think

153:51 answer is no, we think the of gas here is at the level

153:56 1% and probably not commercially attractive to to produce it. But who

154:03 maybe um maybe the they would decide the future that it's uh worthwhile.

154:13 I should tell you that in the couple of years BP has sold this

154:17 to another company where they have only uh ownership of that company and that

154:25 is uh trying to get the last billion barrels, few barrels out of

154:30 . I should say that when this was discovered, it was a 1

154:35 barrel field. Now, it's recognized be a 5 billion barrel field because

154:40 can, we can see it better modern techniques including converted wave techniques and

154:46 better imaging techniques. Um Of NMO processing NMO imaging would be considered

154:54 extremely naive today. But in uh in 1996 it was only a bit

155:02 and the fact that it was successful didn't well its parts everybody but uh

155:08 not the uh the surprise came from solving the the sea wave problems,

155:14 from solving the imaging problems. um let's have a quiz. Let

155:27 turn to you Verda uh true or . It says for segment rocks,

155:34 dominant mechanism of continuation in the se in P waves is fluid squirt.

155:41 that true or false? It is . Yeah, that's true. Uh

155:46 you Lily, for S rocks, dominant mechanism of a generation in the

155:52 band for P waves is enhanced if is present in the pore space who

155:57 false. That's also true that that's abundantly true at val Hall. So

156:04 Carlos for a generative sedimentary rock, waves might, may give images,

156:11 give images which is superior to those B waves, true or false.

156:16 would say it's true. And I gave you an example. And so

156:20 are many other examples. Uh And uh uh it turns out that converted

156:28 are useful for a number of uh solve a number of imaging problems of

156:34 this one is the most prominent, uh uh the gas cloud problem.

156:38 are others which are uh less um common, you know. But uh

156:46 working with us now for 20 we can say that uh converted wave

156:51 are useful mainly in the context where have overburdened gas. So let us

157:00 turn our attention to apparent attenuation. remind you what we mean by

157:06 This is a situation which does not converting elastic energy to heat. This

157:14 a apparent a generation and remind you this slide coming back from uh lecture

157:21 that as we were deriving here uh the wave equation, uh We uh

157:28 found that uh if the, if medium is uniform, we can take

157:35 stiffness tensor outside of this operator, it over here and then what's left

157:41 gonna turn out to be um uh the wave equation I'll back up

157:50 Um um This is not really the equation you can see here that there

157:54 three derivatives with the spec to space of two, you know that uh

158:00 but uh uh if you go back lecture three, you'll see how we

158:05 uh uh found the wave equation scalar equation with only uh one not three

158:13 of 27 term, but one equation for the uh uh scalar potential.

158:21 And that's all the story back in lecture three. So uh also from

158:28 same lecture, we realized that if the medium is non uniform, we

158:34 that problem by uh including another term the uh result. Uh If this

158:40 depends upon space, then we're going have um AAA derivative of that thing

158:47 respect to space. In addition to term, this term is gonna lead

158:52 uh uh the wave equation. And is a new term which comes from

158:58 non uniformity of the media. So now we're gonna see how this leads

159:03 apparent attenuation. OK. For a isotropic medium with vertical P wave

159:11 this becomes like this. So here have only the uh uh the,

159:17 vertical component of displacement. So we that lowercase W. So that's here

159:23 that's here and it's also here. you see this is only one derivative

159:27 respect to depth. And over here a derivative with respect to depth of

159:33 longitudinal modulus. And so this looks if, if we ignore this,

159:39 this is a zero right here, a locally uniform medium. So we

159:44 the wave equation here. But uh now we're gonna consider cases where this

159:50 not zero. So this is the equation and this term is going to

159:55 to apparent attenuation. So let's assume we have plane wave solution. And

160:02 think you're familiar with this and we only uh a K three component of

160:07 wave vector here because uh uh we're propagating in the vertical direction. So

160:15 put that into the equation of not the equation, not the wave

160:19 . We put that into the equation motion, which has this additional term

160:23 here. And um uh so, from two derivatives with respect to

160:30 we get uh a minus uh we I omega squared and from two derivatives

160:36 to K three, we get minus three squared. That's what we saw

160:41 the uh uh uh from the wave . And then we have this additional

160:48 where um uh we have um a was back to Z and only one

160:56 of um um of the displaced uh . So uh simplifying that we get

161:05 quadratic equation for the vertical uh component the wave vector. And it has

161:12 term in there which we didn't see . So uh let's assume that the

161:21 uh the um the stiffness uh the stiffness M is real that way

161:28 there's gonna be no true attenuation at . OK. So whatever comes out

161:34 this, whatever attenuation comes out of is gonna be apparent attenuation. So

161:40 everybody here knows how to solve a equation. I think you remember that

161:45 high school days. And it looks uh this term here, close to

161:50 the square root of this term Now, where does this come

161:58 This comes, I'm gonna back This comes from right here. That's

162:06 ordinary solution to a quadratic equation is take the coefficient of the linear

162:13 Here's the, here's the second order , here's the linear term, here's

162:18 uh the term with no unknown at in it. And so uh uh

162:24 answer to the um to the solution the quad, that quadratic equation is

162:31 here. You see it's um Even though this thing is assumed to

162:37 real, we get uh complex. Now, in this expression, the

162:44 and the stiffness are properties of the not of the wave. Where do

162:48 see the properties of the wave Well, the frequency, that's part

162:51 the wave, that's, that's not of the uh uh uh of

162:59 that's not part of the media, part of the wave. And so

163:03 the wave is specified by the frequency by the choice of algebraic sign in

163:09 second term, you see, I a typo here, I need to

163:12 back and, and remove that type . And so we're gonna consider propagation

163:18 the plus Z direction. So we're um let the plus. Yeah.

163:27 uh let's now um use our old to tailor approximation. So, in

163:34 case, we're gonna repla we're gonna this um square root with this term

163:42 , which has the same minus sign you see right here and everything

163:48 that's what you see here. But you go, uh what you,

163:50 new is the one half that comes the fact that this is a square

163:55 or one half. OK. So gonna average this over a wavelength.

164:04 so um uh um uh thi this is the average here, that's uh

164:13 swap places. So this is the here. That's the real because you

164:17 this is the real part here, is the imaginary part. So I'm

164:21 move the imaginary part over here. uh we, we, this is

164:26 , the real part of the average of the wave factor as it goes

164:32 this a gen homogeneous sequence. And is the imaginary part. So this

164:40 the friendly multiple delay that we talked before. If you look uh uh

164:46 at uh look back what we did . You see that we had uh

164:51 uh uh a real part of the factor depending on uh uh the square

164:56 the uh of the reflectivity. And is the same friendly multiples delay written

165:05 different language, different uh terms. again, I see I made a

165:10 out here. I forgot to have , a closed quote. Right

165:18 So now the plane wave solution looks this. This is our definition of

165:25 plane wave solution got an uppercase W front and the lowercase W means the

165:31 . And this is the amplitude. so we're gonna separate out the real

165:36 parts and the imaginary part is over . And you see that because of

165:44 , I multiplied by this, we have uh high squared minus one

165:51 here. So this is decreasing as um uh as the wave propagates

166:01 So this is not like a surface , a sur surface wave oscillates as

166:06 goes along um horizontally and it decays the amplitude de decays away vertically with

166:14 term looks more or less like But here we have decrease in amplitude

166:23 it goes down. So that's making attenuation. So remember this is

166:34 that's real. Everything is real here this term, but we do get

166:40 situation because of this factor with the sign here. Now let's consider a

166:45 where the stiffness shows a trend. uh uh and that means that uh

166:51 is gonna be greater than zero. this means that this thing is gonna

166:58 greater than zero exponential decay of amplitude of frequency, no true attenuation because

167:08 is real. OK. Now, assume suppose the wave is moving the

167:17 way coming back up. So in propagation in the minus Z direction,

167:22 gonna select the minus sign, same uh uh uh they have the same

167:29 for the uh imaginary term and uh the same thing but the uh you

167:34 , uh variation is small and get of the square root. And uh

167:40 so, um now, at the case, same case where the stiffness

167:50 a positive tr the apparent A generation still um um the same as we

167:58 before, but now Z is the wave is coming back up.

168:04 now, uh um yeah, um the trend is positive, this part

168:10 positive, but as the wave goes smaller Z here, the amplitude is

168:16 exponentially. Wow. So if it down at lost amplitude, but as

168:22 comes back up to the same sequence rocks, it recovers the, the

168:29 isn't that strange. So this apparent cancels itself out for a reflection problem

168:38 A VSP which is, you one way propagation uh that would be

168:44 would be a different scenario. But this case, the apparent cancellation cancels

168:53 . That's interesting. I'll bet that as a surprise to came as a

168:56 to me when I was working through . So now let's consider a case

169:01 the, the stiffness fluctuates, but has no tr so then the imaginary

169:06 of P three vanishes because this, part uh uh averages to zero because

169:12 , there's some places where in increases some places where it decreases, no

169:17 at all. So that means that part is zero. Now, let's

169:25 at the real part, the real is frequency dependent with a nonzero part

169:30 because we're squaring this driven for downgoing , it looks like this and then

169:37 square is nonzero. OK. uh uh that means we're having the

169:43 multiple delay effect when we have uh fluctuating stiffness with no tramp. It's

169:56 friendly multiple delay because of that minus here. Now, let's assume,

170:02 let's be more specific, instead of saying fluctuates, let's say there's a

170:07 uh variation in the stiffness. And let's describe the stiffness as a constant

170:14 of stiffness with uh uh uh uh attitude, uh a verbal part which

170:22 oscillates uh uh with a bit thickness H over two. So that's a

170:27 to describe a, a sinusoidal variation um a nonzero average in the

170:37 So then the stiffest derivative looks like where this part goes away. And

170:42 have only the uh the uh the the, the delta M, that's

170:47 delta M here. And so the looks like that and then the wave

170:52 looks like. So, um uh so this part here, um excuse

170:59 , oh yeah. Uh the way it has only uh see that's uh

171:20 the wave factor is looking like. yeah. So I noticed that we

171:25 the uh the frequency here. And at high frequency, that term goes

171:33 zero because this uh uh this frequency large. And so we're left with

171:40 this term here. That's the average uh at the average of the

171:47 This is the ray theory average and is the average of the slowness of

171:51 high frequencies. That's what we found lecture three for the ray theory.

171:59 result. And so um no, this is the high frequency wave

172:10 So uh recognizing that the velocity is to the wave vector by this relationship

172:17 , we find that the velocity is by the average over the layers of

172:22 inverse of velocity. In other slowness. And then of course,

172:28 have to take the inverse of This is the ra ra result which

172:32 found previously. So I'm gonna show back here, remember this slide here

172:37 we derive the theory by uh uh derive the velocity in a coarse layer

172:43 sequence, we call it coarse layer the layers are thick compared to the

172:49 wavelength. So this is the high limit and we found the same

172:56 Now, going back to uh back this situation, I'm gonna go back

173:02 slide back on, slide back So we're looking at this. Now

173:06 gonna look at the low frequency situation this number is low, here's low

173:16 . And uh um um uh we not gonna neglect this because this is

173:22 low number and the low frequency phase is low uh because of that.

173:33 , um a second here, Um This is the low frequency wave

173:43 , see. Uh And so we're gonna neglect this. We're gonna keep

173:47 , we're not gonna let it go uh a frequency go to zero,

173:51 we're gonna say it's low. So still in there. And uh because

173:57 this term, the low frequency velocity um um less than the high frequency

174:09 as the friend, the same friendly multiple delay. So now that's all

174:17 the propagation. How about the Well, you know that the uh

174:25 that uh we have the, the for attenuation is related to dispersion by

174:32 formula here. So let's apply that the previous situation here. Uh We're

174:39 take the derivative of the low frequency with respect to frequency and we're gonna

174:47 get nothing here because the uh the frequency part does not depend on the

174:53 , but there is a frequency dependence in here. So, working through

174:58 derivative, that derivative is given by to the apparent two factor from the

175:11 representation on page 53. So uh have your uh your notes go back

175:16 page 53. You'll see that the uh the apparent attenuation is related

175:23 the velocity dispersion by this term And so we get friendly multiple attenuation

175:31 well as friendly multiple delay, higher decay more rapidly than low frequencies

175:37 And you know that this doesn't come converting plastic energy to heat. It

175:45 from superposing many, many friendly multiples each other with different time delays.

175:53 they uh uh those different time delays the effect of killing out the high

176:01 . And that's what we, we here for the, uh that's what

176:05 found here for the velocity. The is uh is slower than the high

176:12 , high frequency velocity. And the the amplitudes are also down because of

176:18 . So when we get our display seismic arrivals on our workstation, we

176:30 that the uh uh uh uh at reflection times, we see that

176:37 the frequency content is less and that's to a combination of two attenuation and

176:45 attenuation. This is the apparent attenuation . Oh um La Cross um consider

176:59 quiz statement of uh uh begins with equation of motion, not the equation

177:06 not the wave equation, but the of motion with that additional term leads

177:12 in the case of cyclical bedding ABC A all of the above.

177:18 uh uh number one, let me to you uh Bria, how about

177:24 ? Uh uh is this true lower at higher frequencies? It is

177:31 Hm Who is that true? what the friendly multiples are gonna make

177:38 a delay. So the, the multiples are going to slow down

177:45 uh, um, uh the low waves. So it's gonna be at

177:50 frequencies at lower velocities. So this is false, you've got fueled by

177:59 trickery in the question. So, , le le how about Lee?

178:04 this one true? That one is . So, uh, so this

178:11 better be false. Um, uh sure enough, uh, that one

178:14 false. The uh just the fact we have cyclical betting does not necessarily

178:21 we have attenuation unless the beds themselves a generation in them. The,

178:29 um the apparent uh vellos of high due to the cyclical bedding leads to

178:41 current degeneration without any loss of energy heat. It does lose high frequencies

178:47 it propagates. It doesn't lose any that to heat. It just loses

178:51 to the superposition of the many, friendly multiples on top of each other

178:57 verbal uh um verbal delays that makes attenuation even when there is no true

179:06 . So, as conclusions, uh have learned that a generation is an

179:12 part of any realistic seismology. It the relative loss of high frequencies.

179:20 always accompanied by dispersion and its velocity upon frequency, many different physical

179:31 But in the seismic toson band, most important is fluid squirt. Here's

179:39 you probably didn't see coming at all you have a generation differences at reflectors

179:45 causes a phase shift of the reflected . Now, this might be important

179:52 because we know that for partially saturated below the top of the reservoir,

179:59 gonna make the and a generation difference abo uh for P waves, so

180:06 above the reservoir, we're gonna have attenuation inside the reservoir, we're gonna

180:12 high attenuation. But for that reflected , it doesn't go through from the

180:17 reflected wave from the top doesn't cause loss of high frequency as it says

180:25 , because that wave didn't go through highly reflective, highly attenuated medium,

180:31 goes back up to the normally attenuated . And so, uh even so

180:37 can detect that difference in uh uh uh attenuation at the reflecting horizon with

180:45 resolution, I might say because of phase shift. And so that could

180:51 uh a, a big economic We could use this fact, find

180:55 lot of uh gas which we didn't before. If we look more carefully

181:03 uh historical data, I think we do that. I, I think

181:06 is worth uh uh uh a student um thesis, not sure if it

181:14 be worth it for uh Schlumberger to versa to work on that or for

181:23 uh uh Carlo's company to pay in work for that because these two people

181:27 highly paid but uh uh uh students paid hardly anything. So uh uh

181:35 would be this is a great project a student to figure out whether uh

181:40 kind of uh reflection off the top a oh yeah, gas saturated

181:48 Does that constitute um an expiration clue not? If the answer is

181:55 then we can try maybe make a of money. That person will become

182:00 if the answer is no, because confusion from other effects like uh a

182:06 thin bed, nearby thin bed If, if you can't reliably um

182:13 use this phase shift on reflection, it's still an interesting thesis. And

182:27 , uh when we have lots of uh layers like we always do,

182:31 causes a pa uh uh uh parent . And uh uh uh uh in

182:39 case where the uh the uh the cycle call, we can see very

182:44 that uh that leads to apparent attenuation well as uh real dispersion.

182:51 So I'm sure you have lots of along these lines. But um we

183:01 enough time now, I'm gonna ask to hold your questions, send them

183:05 me by email for next Friday. I'll, I'll tell you what's gonna

183:12 . Now, we're gonna uh we 40 minutes left. We can begin

183:17 discussion of anisotropy. That's a big . I give five day courses,

183:24 full day courses on that. You get a half a day. So

183:28 see that it's gonna be a full a day, next Friday afternoon.

183:36 the end of that day, you not be experts in seismic uh wave

183:43 in case of seismic uh anti but you will know some important

183:49 And so I hope it, it you to learn more about anisotropy because

183:56 ? Because the rocks are an isotropic especially these days when we're looking at

184:02 with angle, we have to consider and we have, of course,

184:08 have variations with angle of incidence and have variations with angle of asthma of

184:15 propagation. And all of these are which were completely um off the radar

184:23 I joined Amao in 1979. But was lucky and I found uh a

184:32 the first week I was with The very first data set that I

184:39 showed the effects of anisotropy. And of my background, I could see

184:46 when others couldn't see it. Other who were much more expert than me

184:51 all of the topics of exploration to , maybe they had been working for

184:58 for 10 years. They were real class experts, but they had not

185:03 their minds pre prepared to think Maybe those subserous rocks are Amish

185:10 When we look at it, we find it and um come on before

185:21 uh before about 1979 or so, 78. Only a few people in

185:29 uh the world thought about anisotropy and people were largely regarded as cranks,

185:36 was paying any attention to them. so now anisotropy is a mainstream

185:42 Uh Lots of people know about We have lots of experts around the

185:47 uh uh and uh um mainstream. you will just get your toe dipped

185:56 this in the next um lecture uh will happen on Friday. Also on

186:07 , you can anticipate receiving from me examination and here's the way it's gonna

186:15 . I'm gonna send you the exam email. I don't like to do

186:19 because, um, I like for to have it as a piece of

186:25 copy in your hand, which is an envelope. So you don't see

186:30 and you don't see it until you to take the exam. But

186:38 uh uh the exam is gonna be book unlimited time. So you can

186:45 your notes, you can have you can have anything open while you

186:49 it. Take as much time as want to do this exam. I'm

186:54 try to make it where I think can do it in three hours.

187:01 usually, well, most people take than three hours, but I do

187:04 best to make it a three hour . And in court, you

187:08 the students could rush through it. might make their first pass in two

187:13 . But then they say, I have some more time, let

187:15 check this and let me check And so that checking process usually ends

187:20 to be more than three hours. when you take the exam set aside

187:25 time for yourself where you're gonna be and you're not gonna have other people

187:33 around bothering you with other things and have a block of time,

187:39 a generous block of time for you work on this. At that

187:43 you open the exam from your email until this time. Uh uh I

187:49 tell you by the way, this is gonna be in um um uh

187:54 attachment to the email. And so uh I think it will, it

187:58 be a PDF attachment. And um the smart thing for you to do

188:07 uh uh as soon as you open up, soon as you open the

188:11 , print off the attachment and then with it. Uh And uh with

188:17 and pencil, there will be plenty space on the paper for you to

188:21 whatever you need to do to solve questions. And then when you're

188:27 uh you scan it and send it to not, not send it to

188:32 , send it to uh Utah. when he gets all three of

188:37 he's gonna send them to me at time. So you can choose your

188:41 . The uh the due date I will be on the Wednesday following the

188:46 . So a Wednesday 10 days from . That's when it will be due

188:50 the end of Wednesday, midnight on Wednesday. But we're doing this by

188:55 because Carlos is over there in, , um, uh, Colombia

189:01 uh, uh, it, we've to send him the exam by

189:05 So that means we don't want it be an unfair, uh, advantage

189:10 disadvantage that he has. So we're do it all by email, uh

189:15 to le le here at the And as a bonus, uh uh

189:20 pres it does not have to drive to from Richmond. She'll get it

189:24 her email. Everybody will handle it the same way after dawn. Right

189:29 there. I took 3.5 hours. I took 2.5 hours, whatever.

189:35 uh, you, uh, wanna take, you take it,

189:39 10 hours if you want, write off and there'll be a space on

189:42 exam to say what, how much you took. That's for my information

189:47 there's no points taken or removed from , uh, the scores because of

189:55 time. I just, I'm interested the time. Send it back to

189:59 by email when he has all three them. He will send it to

190:06 . That's gonna happen after class on sometime Friday evening. I'll send it

190:12 all of you. I know your address. So I'll just send to

190:17 so that leaves us with half an to begin the discussion of anti

190:22 So, what I'm gonna do is gonna stop sharing here and I'm going

190:29 bring up the powerpoint and then I'm watch and I start to be,

190:51 sure did. I'm gonna put myself presentation mode and then I'm going to

191:01 this and I'm going to bring back Zoom session and I'm gonna share my

191:11 and there is less than 10 We're almost done here and share

191:19 So professor, we just don't have in the canvas. That's right.

191:27 not in the compass, but it be um uh uh tonight or uh

191:33 or maybe tomorrow. I, I put it in canvas tomorrow. And

191:38 , um um sorry, uh uh have to just watch the screen here

191:43 anything here. And uh again, having difficulty getting this to me full

191:53 . Um uh Utah help me. do I can't remember what it is

191:58 you do to make this full You go to that one.

192:04 That's it. So everybody sees this less than 10. OK. So

192:10 uh we do, we have uh time left before quitting time today to

192:15 to begin this topic. So, the uh uh the course objectives uh

192:23 . Are you gonna uh be able explain to your friends? Um uh

192:29 are the common classes of anisotropic rocks how to find the wave equation.

192:36 Remember that if it's uh if it's , then it's the same in all

192:41 . No problem. But if it's isotropic, it's gonna be different in

192:45 directions. So that is uh that why we have different classes of anisotropy

192:53 it says here. And then we're find a wave equation. We're gonna

192:57 what we did before uh to find wave equation and then we're gonna find

193:09 and then we're gonna find simpler So if we have, we're gonna

193:16 that if the anisotropy is weak, the solution, the expression for the

193:21 as a function of angle is gonna a lot simpler than these exact

193:27 As a matter of fact, when see these exact solutions, you will

193:30 quite unhappy, but you will be when you see these simple solutions.

193:35 in fact, this is a kind geophysical approximation that we often make.

193:40 recognize that most of the oil that's been found has been found by ignoring

193:46 , right? So uh uh it be AAA small effect. And so

193:55 assume that it's small but non that's what this is. And then

193:59 find that even though it's small, gonna have uh uh a noticeable effects

194:05 our seismic data. And we're gonna uh uh a particular particular attention to

194:15 effects of an I start on the Avio problem. And then uh

194:22 basically, all of this is gonna for P waves. But then we're

194:26 gonna sh show how it affects cheer and it's gonna be a fundamental

194:34 And the reason is because we're gonna that in anisotropic media, there are

194:39 different, share ways, still only P wave but uh two different sho

194:45 . So that's a fundamental difference. then because uh C wave is a

194:50 of PNS, uh there's gonna be effects there. And that's the program

194:56 can see it's a full program. uh It takes me five full days

195:01 instruction to explain all this uh You get half a day. So

195:07 uh let's see what we can learn high points. OK. So everything

195:14 done, the first nine lectures have classic discussion of seismology equally suitable for

195:22 or the deep interior. But none it is truly suitable for either one

195:29 previous discussion has ignored anti our friends do, who are interested in the

195:35 earth. They uh mostly ignore not all of them anymore. Um

195:42 Anisotropy is getting to be a more topic in um uh academic geophysics.

195:50 shouldn't say it that way. I say uh geophysics which is motivated by

195:56 about the earth in exploration Geophysics. motivated by uh we wanna have a

196:03 driven geophysics who motivated by finding by geophysics to solve socially important problems.

196:16 , you know, this of that hydrocarbons are normally found in se

196:20 Iraq, such rocks are normally an . Now, the anti syncopy comes

196:27 these different features. It comes from bed layering. You can see uh

196:31 in this outcrop, many thin And so the the seismic wavelength is

196:38 stretch from the top of this cliff the bottom of this cliff. And

196:43 uh we can say that in uh this layer in this zone, it

196:48 uh propagates as though it's uniform. gonna be reflecting off of boundaries like

196:57 . Uh But uh maybe this is , but you can see the

197:01 many layers in there. And so layers are gonna mean that this uh

197:07 velocity is traveling with an average velocity there. But it's obviously gonna be

197:12 different average if it's a if it's vertically compared to obliquely compared to

197:20 Now, among these layers, there's be shales and the shales are gonna

197:24 intrinsically an iso and then furthermore, gonna be other features in the rock

197:32 are s have a small scale compared the seismic wavelength with the preferred

197:39 And so you can see this joint here. It's it uh it has

197:44 fractured the layers here to here. looks like it didn't penetrate into here

197:50 maybe it did uh did, didn't here. So it's limited top and

197:55 by um uh pathology considerations, but can be very long as you,

198:03 , if we were able to dig into this rock, we would find

198:07 this uh, uh joint goes into rock quite a long ways. How

198:12 we know that? Well, look here, here's another one and another

198:16 probably oriented in the same way. look here is half of one bet

198:21 didn't see this. This is half a fracture and the other half has

198:25 into the lock. This happens to in Ireland. So this is a

198:30 on a lake and you can see here, same thing and these are

198:34 uh are more or less parallel So these are small scale structures smaller

198:40 a seismic wavelength with a preferred And these are gonna lead to uh

198:48 musical variations in in uh and uh . Obviously, if they had just

198:55 fractured layers, then all the asthma would be the same. But when

199:00 have uh uh joints in here, obviously gonna make Lauth and Iar.

199:08 , mother Nature tells us that all masses possess a fabric. So let's

199:15 at this set of hand scents. is a set here. You can

199:18 the layers in here. And so can, you can see obviously that

199:22 gonna have uh uh velocities which are vertically than horizontally just because the fact

199:28 the layers are all horizontal. And this is a small scale version of

199:34 laying that we saw in the left in the previous figure. So let's

199:42 at this. This is a piece metamorphic rock and it has layers

199:46 But these layers are not necessarily As a matter of fact, these

199:51 are uh these are horizontal because they but during the time when this sedimentary

199:58 was being created, gravity was always downwards. And so the layers are

200:04 uh horizontal layers. Gravity is not big issue here. These layers which

200:10 see in here are um perpendicular to uh a gradient and not a gradient

200:20 the pressure but a gradient in So as this mor metamorphic rock

200:25 it made uh different layers perpendicular to gradient of temperature. Furthermore, you

200:32 see the uh the outside shape of . This is a paperweight by the

200:37 and the craftsman who made this paperweight on there, some smooth sides here

200:42 here and here and here. And did that with his machinery is uh

200:50 AAA Craftsman's machinery. And uh uh outer faces are not necessarily related to

201:01 inner structure of the sample. But over here at this crystal. This

201:10 looks like it's homogeneous, matter of , you can see right through it

201:16 you uh uh the it it's got crystal faces here and here and here

201:21 those are determined by the internal distribution the atoms, the atoms here are

201:28 up in uh in uh atomic all lined up with the cells are

201:35 sh uh shaped in this same shape as the, the external uh uh

201:42 . So the the external phases of thing are determined by the internal

201:52 And so uh when I uh come the lecture tomorrow, uh on

201:58 I will bring with me this piece rock. And uh we will see

202:03 least those of us here in Houston see that uh the uh we can

202:09 optical an iso every it by looking through this crystal. Now, here's

202:16 piece of wood. And so this not a rock but uh also as

202:21 uh uh uh uh tree grew, grew, gravity was always pointing

202:27 And so as the tree grew, arranged its internal structure, the the

202:33 of uh the, the wood here arranged so that it's strong in this

202:39 , not so strong in this Uh You can see the, the

202:44 uh the branches came out here. here those are the hor branches which

202:48 gonna lead to horizontal branching, but trunk of the tree is vertical like

202:53 . And uh it's stronger in this than in uh the horizontal directions.

202:59 you've ever worked with wood, you about how you need to respect the

203:03 of the wood as you work as you work with it. And

203:07 can guarantee you that uh uh if have a house, if you live

203:12 a house where the wood frame, always arranged so that the wood is

203:16 out of the tree so that the strong direction of the tree is

203:22 in the strong direction of the lumber hold up the the to hold up

203:28 him. Uh house better. here's a piece of sandstone and this

203:35 , in fact, Baria sandstone, talked about that several time before you

203:39 at it. And you do not with your eyeball, any preferred any

203:44 structure or maybe you can maybe rarely the grains in here, but they're

203:49 randomly oriented. So this rock is except that when you measure it,

203:55 not really isotropic. If you measure uh uh with uh send uh he

204:01 down this way, they travel faster P waves across this. The reason

204:08 because this core was taken out of ground where the coring tool is a

204:14 core. And so the sides of are free and it's a tall

204:19 So the, the stress, whatever was in the rock mass before this

204:24 uh uh taken out of the rock . It has released the stress on

204:30 sides, but it hasn't released the on the top. So that means

204:34 cracks have all open inside this sandstone they are generally cracks with um um

204:46 are not radio cracks, they are ranks. So if you think of

204:52 as uh flat cracks, the the flat sides are always perpendicular or

204:58 perpendicular to this external radius here. uh if you look down on the

205:06 from the top, uh you should see the top edges of the

205:09 you don't see the flat sides of crack. So since the cracks have

205:15 small scale structure with preferred orientation that leads to anisotropy. And by the

205:21 , we verify everything I said here the craft by squeezing the rock.

205:26 when we squeeze it back to its pressure, the anisotropy goes away,

205:31 proves that it's due to uh You can't see them with your

205:38 you can't even see them under a uh because you make two dimensional slices

205:43 your microscope. So you don't see crafts, but you see them with

205:48 sound waves vertically as opposed to Now, when you look at any

205:56 , here's a typical outcrop and it like these rocks are homogeneous. And

206:02 see here, doesn't that look to a homogeneous uh the shales? Uh

206:08 should back up here. You see cliffs here where the rock has uh

206:13 is not shale, but it uh enough calcite cement to hold itself up

206:20 a vertical cliff. But over here see uh uh the rock has

206:25 So if you go, if you away all this stuff here, you'll

206:30 a shale here. Uh uh And reason I know that is because it

206:35 in this pattern with uh uh uh we call this an alluvial sand,

206:42 contrast, this is not shall because has a vertical cliff. So um

206:49 that's what it says here, behind one of these slope, sloping debris

206:55 that's a shale. So uh uh at this section, most of this

207:00 is shale. Now, I call a typical outcrop. Of course,

207:05 is not typical. This is one the most spectacular outcrops in the

207:09 It's a mile from the top to . This is the Grand Canyon.

207:14 , let's look behind uh this uh slope, find the uh the un

207:20 shale and we find um um that looks on a small scale like

207:27 So you can see grains of but most of it is grains of

207:33 and the grains of clay are shaped platelets. And during the sedimentary

207:38 these platelets fell out of the ocean onto the sea floor and mostly

207:46 very few of them landed on their . And then the there were mixed

207:50 a few equipped grains of quartz or Planica here and here. But you

207:56 tell from the orientation of the clay that the original um I don't know

208:03 of gravity is that maybe when you at this, maybe this is really

208:15 originally up in, in in this , maybe this was the top and

208:19 was the bottom. But I'm sure none of you thought that this direction

208:23 up because you can see the orientation a Lady uh Grant. And so

208:28 is gonna mean uh intrinsic anisotropy much the crystal that we saw. All

208:37 these uh uh uh grains of course crystals. And so the the crystal

208:44 anisotropy of this uh clay platelet also to the anti isotopy. Here's a

208:50 of quartz, I think this is , not sure either quartz or

208:54 either way this grain is anisotropic in . But when you look at all

209:02 grains, they are uh all uh oriented so that they yeah,

209:09 The anti oxy internal anisotropy of these does not of these equant grains does

209:16 contribute to the anti oc of the because they're all ran pointed because their

209:24 shape is more or less the same all directions. You can see that

209:29 and here however, that, that platelets are thin and flat and the

209:34 sides are always uh horizontal. And , each of them is an anisotropic

209:41 itself intrinsic an itself. So this is an isotropic and like we said

209:49 , most of this section is I think it's true that something like

209:54 of all the sedimentary rocks in the are shale and a small percentage are

210:02 and a small percentage are carbonates. . So now how are we going

210:10 deal with these anti sat topic We can't use the isotropic um wave

210:20 . I think I misspoke that we, we cannot use the

210:25 We request, we're gonna need the wave equation. And so let's try

210:34 separate out P waves and sheer waves assuming that this stiffness constant is

210:44 Uh I just now called it a constant. This is really a stiffness

210:50 three by three by three by And you know that we enjoyed representing

210:58 on a piece of paper as a by six matrix. So all the

211:04 in that six by six matrix what uh uh is uh represented by the

211:10 uh to me, let me back up. All the information in this

211:16 rank tensor is represented in the second matrix. So that second uh

211:26 we can plot on a piece of . And I'm gonna show you some

211:30 and think about it since it's has two indices, we can't even think

211:35 this because it has four indices. , if you remember that, uh

211:40 we looked at that six by six for isotropic rocks, there were only

211:46 uh there was a uh a lot zeros in there, there were three

211:52 and three muses and three lamin out that, we get two P,

211:59 types of waves two waves and sheer , we never find uh a wave

212:04 is propagating according to LAMBDA, we a wave which is propagating according to

212:10 and one propagating according to view, are respectively PNS waves. Although um

212:17 the LAMBDA is just another concept which in there and is useful for uh

212:24 some uh rock mechanics problems but not anything involving wave propagation. No.

212:34 did we get to PNS waves? , we had that simple stiffness matrix

212:40 , and uh from that stiffness we made a simple isotropic stiffness

212:46 And we put that into this equation . And we used Helmholtz uh uh

212:53 uh which says that whenever you have like this quantity um side guinea um

213:05 , what hes theorem says is that uh whenever you have a vector field

213:12 varies with XY and Z, you always separate it into two parts.

213:17 what hub Holt's term says. Go and check it. We have two

213:21 , the part which is uh uh uh divergence free and a part which

213:28 curl free, the curl free part the scalar of uh of an existing

213:34 for the scalar potential. And this still applies even though this thing is

213:42 isot we solved uh uh this equation it's isotropic. But now we need

213:52 be more realistic. OK. So us consider the simplest case of

214:00 So this is the simplest case of it's gonna apply to thin un

214:06 fractured shales and un fractured, thin sequence. So that's gonna lead to

214:12 what we call polar anti sore because has a pole of symmetry and all

214:18 horizontal uh directions are equivalent. And because as these uh uh these rocks

214:26 formed, gravity was always pointing down the downward direction to be a preferred

214:33 under that um that assumption where the is polar anisotropic, it has a

214:43 symmetry, you can rotate uh uh uh uh around this vertical axis and

214:51 get the same properties independent of That's what we have assumed here as

214:56 the simplest special case. So right , we've assumed axis one and two

215:05 uh soon to be equivalent. And the way, that's the same reason

215:09 we have uh these two elements, elements down here for the uh uh

215:17 mous. And this element is gonna a controlling P wave, vertical P

215:24 uh uh velocities. This one is be in controlling horizontal by velocities.

215:32 gonna need these others for oblique he velocities, right? Most of most

215:38 our data is gonna be uh traveling . So we're gonna need these three

215:44 then we have two different share ways here. Um uh We'll talk about

215:51 , the shear wave complications, then get to shear wave propagation on

215:59 Now count them up how many um how, how many uh independent parameters

216:06 we have? Well, we have true 345 parameters. So I,

216:16 regret to tell you that the simplest of an isotopy doesn't have just three

216:25 numbers. It has five different So um I I misspoke just

216:34 there is a case of anisotropy which only three different numbers scattered among uh

216:42 the uh elements here. So that's bad. It uh going from two

216:48 two to that other case with That's not bad, but that's not

216:53 case which we ever encounter in That's the case of cubic anisotropy.

216:59 we do have cubic crystals. If uh as you walk down the hall

217:03 class, look at in the uh cabinet there. I think we have

217:07 cubic crystals in that cabinet. And they have very simple anti. But

217:14 in geophysics, this is the simplest that we have. And it turned

217:19 that this one governs the vertical P velocity. This one governs the horizontal

217:25 wave velocity. We have two different moduli with two different shear waves.

217:31 then we have this one here which repeated here, but it's like a

217:36 parameter and it's not gonna uh clo , actually, um this one is

217:44 show up in wave propagation, it gonna show up in wave propagation.

217:49 gonna name it C 13. we're not gonna name it after lame

217:55 Lama didn't know anything about that, only knew about isotopic. So now

218:01 is the simplest case, we're gonna more complicated cases. And by the

218:07 of the day and of the evening Friday, you will be able to

218:12 manage this case and further more uh , more realistic cases. Why do

218:19 say more realistic? Because these shas usually not frac uh not the these

218:26 usually have fractures in them with preferred . That means the shales have lower

218:33 than is indicated here. Uh But is gonna be for un fractured shas

218:38 unraced tin bed sequence. We're gonna this matrix as a first step towards

218:47 an isotropic wave of now. Too . Even for the simplest case,

218:52 quad, it has gonna have a P wave and two quasi estimate

218:58 but they are coupled together. That the curl free part of the field

219:04 he gave us. That's not a . So we need to have a

219:09 math math mathematical idea. But here's beginnings of that better mathematical idea.

219:16 gonna go back to the full tensor of motion and we're gonna have it

219:23 gonna be the homogeneous equation of So uh it's gonna be this is

219:28 be a homogeneous with respect to no spatial dependence on that. But

219:35 not gonna be the isotropic special case we looked at before. It's gonna

219:41 , it's gonna have these kinds of in it. So we'll take that

219:48 on Friday. Uh, uh, we'll be able to read ahead,

219:53 this is a good place for us stop for. Now. Um,

219:57 can see what's, uh, what's go ahead. We're gonna guess we're

220:01 assume, uh, plane waves and , uh, this assumption into here

220:06 see if that works. And, , you will see, um uh

220:11 gets, it's not gonna be as uh simple and straightforward as we found

220:17 rocks. That's the reason why we this discussion of anisotropy until after you've

220:25 your minds around uh uh the uh consequences of assuming is arre I think

220:34 were a lot of complications that we at here uh with uh um isotropic

220:41 . And so now we're gonna see different kinds of complications arising from an

220:48 . And of course, we're driven do that because the rocks are in

220:51 to drink. If the rocks were , we would happily leave this

220:56 Well, let's leave it for the uh for there for today. And

221:00 I'm gonna stop sharing this and I'm uh say goodbye. Uh uh You

221:06 send me some questions arising from the today. Uh during the week

221:11 I will post onto canvas the 10th . I didn't do it before because

221:17 didn't know we would need it But, you know, we didn't

221:21 it very much. You will have shortly. Um, and,

221:25 so tomorrow morning after church you'll have chance to, uh, to read

221:32 . Ok. So I think that is a good time for, for

221:36 to say goodbye and leave you uh, for Fri until Friday at

221:41 .

-
+