Документ взят из кэша поисковой машины. Адрес оригинального документа : http://geophys.geol.msu.ru/shevnin/publ/117.pdf
Дата изменения: Wed Oct 17 15:56:22 2012
Дата индексирования: Sat Feb 2 23:40:01 2013
Кодировка:
Bolshakov D.K., Modin I.N., Pervago E.V., Shevnin V.A.

Modeling and interpretation of azimuthal resistivity sounding over two-layered model with arbitrary - oriented anisotropy in each layer.
MSU, Geological faculty, 119899 Moscow, Russia. Fax: (007-095) -9394963, E-mail: sh@geophys.geol.msu.ru Proceedings of EAGE 60th Conference, Leipzig - 1998. P110.

The basic moment for creation of new algorithm was clear aware of half-space model limitation. Twolayered model with isotropic first layer (Bolshakov et al., 1997) and anisotropic basement has allowed to receive many new peculiarities of azimuthal soundings. The practice has shown that further development of model is required. Therefore an algorithm and program of modeling and interpretation of azimuthal resistivity sounding (ARS) data I over two-layered model with arbitrary - oriented anisotropy in each layer has been developed. The decisions for such model were known earlier (Gurevich, Sagina, 1977; Ping Li et al., 1997). Our purpose was in receiving fast and H effective computing algorithm, which would give precise results. The basis of the decision is spectral approach. The potential value is calculated as a sum of spectral harmonics C2n. At estimation of potential for various azimuths, its harmonics are calculated only once, that increases effectiveness of algorithm. The final expressions for an electrical field potential Fig.1. Two-layered anisotropic model look like:
U (r , , z ) =
Cn (r ) = 1 4 2 r 2 I1,l
n = -





C

2n

(r , z ) exp( -i 2n )
1 m Dn,0 (r , k ) + 2 K (k ) Dn,m (r , k ) b1(k ) m =1 1

-





En ( r , k ) exp(i nk ) dk ,
D (r , k ) =

En ( r , k ) =



n,m

1 + ms( r , k )
= n t , g = sin2 ( ) + 1
2

(

)2

1 + ms(r , k )

(

)2

+ ms(r , k )
cos2 (k - ),

n

g1r The final formulas are similar to the formulas for two-layered isotropic model and just as these allow to fulfill calculations with high accuracy. On the basis of harmonics for potential it is possible to receive harmonics for any component of electrical field and a for any array. The same formulas allow to fulfill transformation of azimuthal data for dipole equatorial array to linear one (D.K. Bolshakov et al., 1997), that is very useful for reduction of deep and near-surface inhomogeneities' influence, for comparison of results for different arrays and their interpretation. To study equivalence problem for the given model, two new parameters ef and Hef were entered. ef = 2 sin 2 ( ) + cos 2 ( ), H ef = 1 H , 1ef Thus, instead of three parameters (, , H) there is the influence of only two effective parameters (ef, Hef) or theoretical equivalence. At measurements on the earth surface it is impossible the estimate these three parameters separately. We calculated a series of two-layered models with an anisotropy in first or in the second layer. We reasoned that the task with anisotropy in both layers makes the model too complex. To study a large series of data the program was supplied with the system of visualization. It allows to see all data for any azimuth depending on spacing or to see all data for any spacing depending on azimuth. Except these the program shows diagrams a = f(r), harmonics of spectrum C2n = f(r), and values l, m, n for sounding curves, transformed to a uniform half-space. These transformations are estimated on azimuthal diagrams' spectra and their introduction is similar to introduction of apparent resistivity for isotropic media. The examples of accounts for two models are shown on fig.2. A number of anomalous or unusual phenomena were found out here. The most complicated forms of sounding curves are found for dipole array at its

cos2 ( ),

b(k ) = g sin2 (k - ) +

1
2

s( r , k ) =

2b1( k )H


Model 1 Model 2 orientation near 45° from A C strike direction. The visible M M number of layers on fig.2, B, D differs from true. The 1 50 1 10 1.8 right and left asymptotes 2 50 1.8 10 1 fall outside the limits of l, D B m values. On fig.3 the values of apparent anisotropy a appear more than true greatest and smallest ones. The values a for D - array on fig.3 were recalculated to E F equivalent linear array R/h (D.K. Bolshakov et al., Fig.3. a graphs for two models from 1997) to compare with Schlumberger array Dipole equatorial array fig.2 linear array's data (E). 50 m From comparison of a graphs for two models with an anisotropy in the first or in second layer, it is visible, that phenomena of Fig.2. Sounding curves and azimuthal amplification (for model 2) and decrease of anisotropy (for model 1) diagrams for two anisotropic models are similar. Case 1 gives the change of ellipse ,Ohm.m C A orientation for E-array (disappearance of anisotropy Z, m 50 0 0 10 paradox). D-array for both models gives smaller 450 30 0 20 30 deviations of apparent and true parameters and has no 20 0 ARS 45 40 change in ellipse orientation. 10 0 On fig.4 an example of interpretation of 500 600 X, m 0 100 200 300 400 R, m 6 0 azimuthal resistivity soundings, received in a place 3 571 0 2 03 0 5 0 10 0 Patil (Ukraine, Crimea) above two-layered earth with ,Ohm.m 50 0 anisotropy in both layers (A), is presented. In the first B D 30 0 layer the bedding (and anisotropy) orientation is close 0° 20 0 to horizontal, while in the second layer these are close to vertical. Due to inclined earth surface cut the 10 0 first layer, its bedding plane appears inclined to the 5 0 earth surface. Interpretation or fitting experimental 90° and theoretical data fulfills and checks on all types of 3 0 R=30 m diagrams; that is on azimuthal diagrams for any 2 0 field data spacing (B), on sounding curves for any azimuth R, m 1 0 model data (along and across the strike on example D) and on 3 571 0 2 03 0 5 0 10 0 spectral characteristics, in particular on harmonics, a Fig.4. Example of practical data interpretation. graphs and l, m, n diagrams (C). The best fitting is achieved when of the first layer equals 4 and dip angle of its bedding =5°. These results are close to geological estimations. Conclusions The decision for a, potential and electric field's calculation for many arrays above two-layered model with arbitrary oriented anisotropy in each layer has been received. Principles of theoretical equivalence of (, , H) parameters for such model are formulated. Some anomalous and unusual effects on sounding curves and =f (r) diagrams were revealed for twolayered anisotropic models. This program is suitable for interpretation of experimental azimuthal resistivity sounding data. References Bolshakov D.K., Modin I.N., Pervago E.V., Shevnin V.A. Anisotropy effects investigations by resistivity method in some inhomogeneous media. Proceedings of EAEG 57th Annual Meeting, Glasgow, May 28-June 2, 1995. P034. D.K. Bolshakov, I.N. Modin, E.V.Pervago, V.A.Shevnin. The new approach to the analysis of the azimuthal resistivity data over anisotropic media. Proceedings of 59th EAGE Conference, Geneva - 1997. P139. Ju.M.Gurevich, O.V.Sagina. Electrical field from point source of current lowered in two layered anisotropic model // Razvedochnaja geophysika. 1977. Issue. 74. P. 37-45. (In Russian). Ping Li and Norm Uren. Analytical solution for direct current electrical potentials in a two layer anisotropic earth medium. Proceedings of 59th EAGE Conference, Geneva - 1997. F021. Schlumberger

Model

1st layer

2nd layer

100

100



1

N

50



1

50

1

0

, m

20 10 5



, m

2

N

20 10 5



1

L

90



a

2M

0

a



2

R, m

2

L

90

R, m

0 .5 1

2

5

10

20

50100200

0 .5 1

2

5

10

20

50 10020

200 100 50

200 100 50



0
1 N



1

20
1 L

Dipole equatorial array

20 10



0



2

N

20 10

30

ef

2M

20

2

L

, m

30

, m

5

5

90



2

40

60

a

a

2

90

2

0.5

0.5

40

50

1

50 60

1

0.2

R, m
200

0.2

R, m
200

0.5

1

2

5

10

20

50

100

0.5

1

2

5

10

20

50

100

a

L

M

N

a

hm.m 0O 20

1 00