CN102608666A - Fast and accurate depth inversion method of transient electromagnetic data - Google Patents

Fast and accurate depth inversion method of transient electromagnetic data Download PDF

Info

Publication number
CN102608666A
CN102608666A CN2012100797683A CN201210079768A CN102608666A CN 102608666 A CN102608666 A CN 102608666A CN 2012100797683 A CN2012100797683 A CN 2012100797683A CN 201210079768 A CN201210079768 A CN 201210079768A CN 102608666 A CN102608666 A CN 102608666A
Authority
CN
China
Prior art keywords
depth
layer
thickness
resistivity
transient electromagnetic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2012100797683A
Other languages
Chinese (zh)
Inventor
范涛
王信文
程建远
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN2012100797683A priority Critical patent/CN102608666A/en
Publication of CN102608666A publication Critical patent/CN102608666A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a fast and accurate depth inversion method of transient electromagnetic data. According to the delay time equivalent simple relationship in a periodic and transient electromagnetic method of a magnetotelluric method, a Bostick inversion method of the magnetotelluric method is adopted for building an initial model, the Bostick depth is taken, and finally, the thickness and the electric resistivity of each electric layer are sequentially adjusted according to differences between measured values and apparent resistivity theoretical fitted values calculated through inversion. The invention has the advantages that the method is simple and convenient, the operation is easy, good processing effects are realized on theoretical data and actual measurement data, the fast inversion can be realized to obtain the depth conforming to the geologic structure, and the data processing speed and the explanation precision of the transient electromagnetic method are greatly improved, so the great development and the application of the method technology are promoted.

Description

The transient electromagnetic data are degree of depth inversion method fast and accurately
Technical field
The invention belongs to the geophysical exploration technology field, relate to a kind of transient electromagnetic data degree of depth inversion method fast and accurately.
Background technology
Transient electromagnetic method is called for short TEM, is a kind of electromagnetism class method of exploration that is widely used in fields such as resource exploration, engineering geophysics, hydrology exploration in recent years.This method popular depth computing method at present mainly contains skin depth estimation algorithm, equivalent conductive plane method and smoke ring theoretical method etc.; These methods are calculated easy; But result of calculation often has than large deviation with actual geological condition, can not well reflect subsurface geological structure; And the generalized inverse inverting, computational accuracy is very high, but because computing method are complicated, calculated amount is very big, so the difficult realization of inverting, and speed is slower, and computing machine is had relatively high expectations, and is difficult to satisfy the needs of actual engineering
In fact; Though transient electromagnetic method belongs to TDEM; But it and frequency domain electromagnetic methods still have getting in touch of countless ties, in the derivation of many fundamental formulars, have all adopted elder generation in frequency field, to discuss; Convert the means of time domain at last into through Fourier transform, so the maturation method that can use for reference legacy frequencies territory electromagnetic method is explained in the processing of TEM data.
Magnetotelluric method is called for short MT; It is electromagnetic field horizontal component through the observation natural variation; Convert electromagnetic field signal to apparent resistivity curve and phase curve, the resistivity on each stratum and a kind of frequency domain electromagnetic method of one-tenth-value thickness 1/10 are tried to achieve in inverting then, and its degree of depth inversion method is simple; Reliable results; And the cycle of this method with have time delay among the TEM simple relation of equity mutually, can easily the TEM data-switching be become the MT data, utilize the MT inversion method to ask for and meet the architectonic degree of depth.
Summary of the invention
The technical matters that the present invention solved provides a kind of transient electromagnetic data that solve the problem that the inverting degree of depth that exists in the existing transient electromagnetic data processing technique is difficult to match with actual geological condition degree of depth inversion method fast and accurately.For solving above-mentioned technical matters, the technical scheme that the present invention takes:
A kind of transient electromagnetic data are degree of depth inversion method fast and accurately, and its special character is: realize through following steps:
(1), according to the time delay in the cycle in the magnetotelluric method and the transient electromagnetic method of the simple relation of equity mutually
Figure 2012100797683100002DEST_PATH_IMAGE001
Can be with the approximate actual measurement apparent resistivity curve that converts magnetotelluric method into of the apparent resistivity curve that transient electromagnetic method calculates;
Figure 361553DEST_PATH_IMAGE002
is the cycle in the magnetotelluric method in the formula, and
Figure 872169DEST_PATH_IMAGE004
is the time delay in the transient electromagnetic method;
(2), adopt Bostick inversion method body plan initial model in the magnetotelluric method, get the Bostic degree of depth
Figure 2012100797683100002DEST_PATH_IMAGE005
Be the resistivity
Figure 2012100797683100002DEST_PATH_IMAGE007
on the initial inverting degree of depth on the corresponding periodic point
Figure 715490DEST_PATH_IMAGE006
, the respective depth point, can convert by the differential of apparent resistivity (
Figure 173016DEST_PATH_IMAGE008
) curve
Figure 2012100797683100002DEST_PATH_IMAGE009
；
(3), adjust each electrically thickness and resistivity of layer successively according to the difference of the apparent resistivity theoretical fitting value of just drilling calculating and measured value.
Step in the above-mentioned inversion method (3) is adjusted each electrically thickness and resistivity of layer through following steps:
A) just drilling the theoretical fitting value of the apparent resistivity that calculates
Figure 393913DEST_PATH_IMAGE010
inferior model; Thought according to your process of iteration of Gauss's one Saden; Here with the 1st layer to
Figure 731222DEST_PATH_IMAGE012
layer already correction result join the process of finding the solution
Figure 857310DEST_PATH_IMAGE011
, promptly
Figure 869259DEST_PATH_IMAGE014
Figure 2012100797683100002DEST_PATH_IMAGE015
functional in the formula for just drilling;
B) refutation process is pressed the resistivity and the thickness of two formula corrections
Figure 2012100797683100002DEST_PATH_IMAGE017
layer
Figure 334876DEST_PATH_IMAGE018
In two formulas,
Figure 639824DEST_PATH_IMAGE020
is the resistivity correction factor;
Figure 2012100797683100002DEST_PATH_IMAGE021
is the poor of the measured value of apparent resistivity and the theoretical value of just drilling calculating; if
Figure 722050DEST_PATH_IMAGE022
; Then individual electrical layer resistivity will suitably increase a bit; I.e.
Figure 2012100797683100002DEST_PATH_IMAGE023
; if
Figure 529917DEST_PATH_IMAGE024
; Then individual electrical layer resistivity will suitably reduce a bit, i.e.
Figure 2012100797683100002DEST_PATH_IMAGE025
;
Figure 732545DEST_PATH_IMAGE026
is the thickness correction factor, handles respectively according to the mutual relationship between
Figure 2012100797683100002DEST_PATH_IMAGE027
,
Figure 584832DEST_PATH_IMAGE028
,
Figure 2012100797683100002DEST_PATH_IMAGE029
:
C) according to a) and b) to each electrically layer circulate, promptly
Figure 290620DEST_PATH_IMAGE017
circulates from 1 to N;
D) according to a), b) and c)
Figure 360207DEST_PATH_IMAGE030
circulated; When the difference of the measured value of apparent resistivity and the theoretical fitting value of just drilling calculating when overall relative error
Figure 2012100797683100002DEST_PATH_IMAGE031
is less than preassigned number
Figure 902178DEST_PATH_IMAGE032
on all periodic points, promptly
To stop circulation.
Mutual relationship between above-mentioned ,
Figure 389889DEST_PATH_IMAGE028
,
Figure 630378DEST_PATH_IMAGE029
is for when
Figure 440071DEST_PATH_IMAGE034
and
Figure 2012100797683100002DEST_PATH_IMAGE035
; if
Figure 525838DEST_PATH_IMAGE022
;
Figure 222530DEST_PATH_IMAGE036
; Promptly the suitable attenuate of individual electrical layer thickness thickens and levels thickness is corresponding; if
Figure 665330DEST_PATH_IMAGE024
;
Figure DEST_PATH_IMAGE037
; Promptly
Figure 663110DEST_PATH_IMAGE017
individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness.
Mutual relationship between above-mentioned , , is for as
Figure 664248DEST_PATH_IMAGE038
time; if
Figure 335531DEST_PATH_IMAGE022
; ; Then this layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down; if
Figure 823145DEST_PATH_IMAGE024
;
Figure 829147DEST_PATH_IMAGE039
; This layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer.
Mutual relationship between above-mentioned
Figure 325987DEST_PATH_IMAGE027
,
Figure 615892DEST_PATH_IMAGE028
,
Figure 274406DEST_PATH_IMAGE029
is for when
Figure 767705DEST_PATH_IMAGE040
and
Figure DEST_PATH_IMAGE041
;
Figure 537077DEST_PATH_IMAGE022
;
Figure 182954DEST_PATH_IMAGE037
; Then
Figure 277948DEST_PATH_IMAGE042
individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness; if ;
Figure 97186DEST_PATH_IMAGE036
; The suitable attenuate of
Figure 830524DEST_PATH_IMAGE017
individual electrical layer thickness thickens and levels thickness is corresponding.
Mutual relationship between above-mentioned
Figure 830841DEST_PATH_IMAGE027
,
Figure 298732DEST_PATH_IMAGE028
,
Figure 675487DEST_PATH_IMAGE029
is for as
Figure DEST_PATH_IMAGE043
time; if
Figure 358272DEST_PATH_IMAGE022
;
Figure 795069DEST_PATH_IMAGE039
; Then this layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer; if ; ; This layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down.
The resistivity correction factor
Figure 836035DEST_PATH_IMAGE020
, the thickness correction factor
Figure 319100DEST_PATH_IMAGE026
, a pre-specified number
Figure 371370DEST_PATH_IMAGE032
value of the relative difference between the measured value more than three orders of magnitude of the number.
Above-mentioned resistivity correction factor
Figure 948982DEST_PATH_IMAGE020
is 0.001.
Above-mentioned thickness correction factor
Figure 403097DEST_PATH_IMAGE026
is 0.001.
Above-mentioned preassigned number
Figure 290019DEST_PATH_IMAGE032
is 0.001.
Compared with prior art; The inventive method is easy, and no matter easy operating is for gross data or measured data; Good treatment effect is all arranged; Can fast inversion draw and meet the architectonic degree of depth, improve transient electromagnetic method Data Processing speed greatly and explain precision, thereby promote the development significantly and the application of this method and technology.
Description of drawings
Fig. 1 is a general flow chart of the present invention;
Fig. 2 is that theoretical model of the present invention calculates design sketch;
Fig. 3 is that theoretical model of the present invention calculates design sketch;
Fig. 4 is a measured data design sketch of the present invention;
Fig. 5 is a measured data design sketch of the present invention;
Fig. 6 is a measured data design sketch of the present invention;
Fig. 7 is a measured data design sketch of the present invention;
Fig. 8 is a measured data design sketch of the present invention;
Fig. 9 is a measured data design sketch of the present invention.
Embodiment
Below in conjunction with embodiment and Figure of description the present invention is elaborated.
Referring to Fig. 1, the present invention realizes through following steps:
(1), according to the time delay in the cycle in the magnetotelluric method and the transient electromagnetic method of the simple relation of equity mutually
Figure 954219DEST_PATH_IMAGE001
Can be with the approximate actual measurement apparent resistivity curve that converts magnetotelluric method into of the apparent resistivity curve that transient electromagnetic method calculates;
Figure 820675DEST_PATH_IMAGE002
is the cycle in the magnetotelluric method in the formula, and
Figure 831094DEST_PATH_IMAGE004
is the time delay in the transient electromagnetic method;
(2), adopt Bostick inversion method body plan initial model in the magnetotelluric method, get the Bostic degree of depth
Be the resistivity
Figure 336658DEST_PATH_IMAGE007
on the initial inverting degree of depth on the corresponding periodic point
Figure 151665DEST_PATH_IMAGE006
, the respective depth point, can convert by the differential of apparent resistivity (
Figure 499787DEST_PATH_IMAGE008
) curve
Figure 462932DEST_PATH_IMAGE009
；
(3), adjust each electrically thickness and resistivity of layer successively according to the difference of the apparent resistivity theoretical fitting value of just drilling calculating and measured value.
Step in the above-mentioned inversion method (3) is adjusted each electrically thickness and resistivity of layer through following steps:
A) just drilling the theoretical fitting value
Figure 575562DEST_PATH_IMAGE011
of the apparent resistivity that calculates
Figure 101724DEST_PATH_IMAGE010
inferior model; Thought according to your process of iteration of Gauss's one Saden; Here with the 1st layer to
Figure 593196DEST_PATH_IMAGE012
layer already correction result join the process of finding the solution
Figure 9134DEST_PATH_IMAGE011
, promptly
Figure 261353DEST_PATH_IMAGE014
Figure 992549DEST_PATH_IMAGE015
functional in the formula for just drilling;
B) refutation process is pressed the resistivity and the thickness of two formula corrections layer
Figure 818870DEST_PATH_IMAGE018
Figure 883778DEST_PATH_IMAGE019
In two formulas;
Figure 876005DEST_PATH_IMAGE020
is the resistivity correction factor, gets the less number of an absolute value usually as 0.001; is the poor of the measured value of apparent resistivity and the theoretical value of just drilling calculating; if
Figure 452535DEST_PATH_IMAGE022
; Then individual electrical layer resistivity will suitably increase a bit; I.e.
Figure 449758DEST_PATH_IMAGE023
; if
Figure 988186DEST_PATH_IMAGE024
; Then individual electrical layer resistivity will suitably reduce a bit, i.e. ;
Figure 198960DEST_PATH_IMAGE026
is the thickness correction factor; Usually also get the less number of an absolute value as 0.001, divide 4 kinds of situation to handle respectively according to the mutual relationship between
Figure 173869DEST_PATH_IMAGE027
,
Figure 983562DEST_PATH_IMAGE028
,
Figure 334909DEST_PATH_IMAGE029
:
Figure 31601DEST_PATH_IMAGE044
is when
Figure 442990DEST_PATH_IMAGE034
and
Figure 474400DEST_PATH_IMAGE035
; if
Figure 363859DEST_PATH_IMAGE022
;
Figure 679172DEST_PATH_IMAGE036
; Promptly the suitable attenuate of individual electrical layer thickness thickens and levels thickness is corresponding; if
Figure 514589DEST_PATH_IMAGE024
;
Figure 207739DEST_PATH_IMAGE037
; Promptly
Figure 879023DEST_PATH_IMAGE017
individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness.
Figure DEST_PATH_IMAGE045
is during as
Figure 225690DEST_PATH_IMAGE038
; if
Figure 107059DEST_PATH_IMAGE022
;
Figure 729800DEST_PATH_IMAGE039
; Then this layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down; if
Figure 645803DEST_PATH_IMAGE024
;
Figure 694530DEST_PATH_IMAGE039
; This layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer.
is when
Figure 973513DEST_PATH_IMAGE040
and
Figure 744023DEST_PATH_IMAGE041
;
Figure 839018DEST_PATH_IMAGE022
;
Figure 819612DEST_PATH_IMAGE037
; Then individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness; if
Figure 657173DEST_PATH_IMAGE024
;
Figure 923069DEST_PATH_IMAGE036
; The suitable attenuate of
Figure 266326DEST_PATH_IMAGE017
individual electrical layer thickness thickens and levels thickness is corresponding.
Figure DEST_PATH_IMAGE047
is during as
Figure 502135DEST_PATH_IMAGE043
; if ;
Figure 559401DEST_PATH_IMAGE039
; Then this layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer; if
Figure 389954DEST_PATH_IMAGE024
; ; This layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down.
C) according to a) and b) to each electrically layer circulate, promptly
Figure 29063DEST_PATH_IMAGE017
circulates from 1 to N;
D) according to a), b) and c)
Figure 276242DEST_PATH_IMAGE030
circulated; When the difference of the measured value of apparent resistivity and the theoretical fitting value of just drilling calculating on all periodic points overall relative error
Figure 328512DEST_PATH_IMAGE031
less than a preassigned very little number
Figure 171703DEST_PATH_IMAGE032
, as 0.00001; The time, promptly
Figure 625818DEST_PATH_IMAGE033
To stop circulation.
The resistivity correction factor
Figure 279784DEST_PATH_IMAGE020
, the thickness correction factor
Figure 819350DEST_PATH_IMAGE026
, a pre-specified number
Figure 607177DEST_PATH_IMAGE032
value of the relative difference between the measured value more than three orders of magnitude of the number.
Table 1 is one group of sample data; Can find out; Change by entering resistive formation position, low resistivity zone position along with
Figure 40433DEST_PATH_IMAGE008
; The inverting degree of depth
Figure 989934DEST_PATH_IMAGE048
shows as a sudden change; Meet the transient field signal and in high resistant, propagate the fast characteristics that get;
Figure DEST_PATH_IMAGE049
classic method compute depth then changes gently, is difficult to reflect this resistive formation position.
Table 1
Figure 310112DEST_PATH_IMAGE004
Figure 203430DEST_PATH_IMAGE008
Figure 717588DEST_PATH_IMAGE048
Figure 706272DEST_PATH_IMAGE049
0.84 4.33 362.711608887 84.5 120.539627075
0.9425 4.858 380.533203125 100.8125 139.2837677
1.0575 5.451 392.518920898 348.125 157.410522461
(
Figure 989486DEST_PATH_IMAGE048
is the inverting degree of depth of the present invention in the table, and
Figure 660551DEST_PATH_IMAGE049
is the classic method compute depth)
Be example explanation effect of the present invention with theoretical model and measured data respectively below:
Fig. 2,3 provides a TEM high resistant model (Fig. 2) and a TEM low-resistance model (Fig. 3) respectively, can obviously find out; Main body in the model is unusual; Inversion result all is superior to conventional result of calculation, and unusual interphase is more accurate, and abnormal thickness also coincide better with model.
Fig. 4 and Fig. 5 are conventional result of a TEM survey line in colliery, the Inner Mongol and inversion result comparison diagram; Several black lines are stratigraphic horizon among the figure; Can find out that zone thickness changes greatly, bigger to the rule influence of original apparent resistivity isoline, coal seam absolute altitude drop reaches three, 400 meters; Vertically going up apparent resistivity value and uprised gradually by low, is consistent with the electrical regularity of distribution of actual geologic horizon.But the whole east of apparent resistivity is higher than western in the horizontal, can not well reflect the developmental state on stratum, and analyzing mainly is by the influence of monoclinal stratum to data.Through finding out that the rule of stratum law of development and apparent resistivity isoline is close on the sectional view after the inversion procedure, the sectional view after the inverting is the situation of illustrative actual formation more.
Fig. 6 and Fig. 7 are conventional result of a TEM survey line in colliery, Shaanxi and inversion result comparison diagram; Can see and compare the conventional processing result; Not only obvious through the section after the inversion procedure to the electrical regularity of distribution reflection on stratum; And given prominence to the unusual of deep formation, and make the space distribution of exceptions area more outstanding, intuitively, can better reflect the low-resistance abnormal morphology characteristic in water guide cracks such as doubtful tomography.
Fig. 8 and Fig. 9 are conventional result of a TEM survey line in colliery, Xinjiang and inversion result comparison diagram, and the position that pentagram indicates among the figure is the actual EXIT POINT in down-hole, on processed conventionally figure as a result; This position, coal seam is the high resistant reaction; And on inversion result figure, then be a unusual border of low-resistance of coming by top board, not only corresponding good with EXIT POINT; And having reflected unusual source, this shows that inversion result tallies with the actual situation more.

Claims (10)

1. transient electromagnetic data degree of depth inversion method fast and accurately,, it is characterized in that: realize through following steps:
(1), according to the time delay in the cycle in the magnetotelluric method and the transient electromagnetic method of the simple relation of equity mutually
Can be with the approximate actual measurement apparent resistivity curve that converts magnetotelluric method into of the apparent resistivity curve that transient electromagnetic method calculates;
Figure 2012100797683100001DEST_PATH_IMAGE004
is the cycle in the magnetotelluric method in the formula, and
Figure 2012100797683100001DEST_PATH_IMAGE006
is the time delay in the transient electromagnetic method;
(2), adopt Bostick inversion method body plan initial model in the magnetotelluric method, get the Bostic degree of depth
Figure 2012100797683100001DEST_PATH_IMAGE008
Be the resistivity
Figure 2012100797683100001DEST_PATH_IMAGE012
on the initial inverting degree of depth on the corresponding periodic point
Figure 2012100797683100001DEST_PATH_IMAGE010
, the respective depth point, can convert by the differential of apparent resistivity (
Figure 2012100797683100001DEST_PATH_IMAGE014
) curve
Figure 2012100797683100001DEST_PATH_IMAGE016
；
(3), adjust each electrically thickness and resistivity of layer successively according to the difference of the apparent resistivity theoretical fitting value of just drilling calculating and measured value.
2. transient electromagnetic data according to claim 1 are degree of depth inversion method fast and accurately,, it is characterized in that: step in the described inversion method (3) is adjusted each electrically thickness and resistivity of layer through following steps:
A) just drilling the theoretical fitting value
Figure 2012100797683100001DEST_PATH_IMAGE020
of the apparent resistivity that calculates
Figure 2012100797683100001DEST_PATH_IMAGE018
inferior model; Thought according to your process of iteration of Gauss's one Saden; Here with the 1st layer to
Figure 2012100797683100001DEST_PATH_IMAGE022
layer already correction result join the process of finding the solution
Figure 110180DEST_PATH_IMAGE020
, promptly
Figure 2012100797683100001DEST_PATH_IMAGE026
Figure 2012100797683100001DEST_PATH_IMAGE028
functional in the formula for just drilling;
B) refutation process is pressed the resistivity and the thickness of two formula corrections
Figure 2012100797683100001DEST_PATH_IMAGE030
layer
Figure 2012100797683100001DEST_PATH_IMAGE032
Figure 2012100797683100001DEST_PATH_IMAGE034
In two formulas,
Figure 2012100797683100001DEST_PATH_IMAGE036
is the resistivity correction factor;
Figure 2012100797683100001DEST_PATH_IMAGE038
is the poor of the measured value of apparent resistivity and the theoretical value of just drilling calculating; if
Figure 2012100797683100001DEST_PATH_IMAGE040
; Then
Figure 89637DEST_PATH_IMAGE030
individual electrical layer resistivity will suitably increase a bit; I.e.
Figure 2012100797683100001DEST_PATH_IMAGE042
; if
Figure 2012100797683100001DEST_PATH_IMAGE044
; Then
Figure 400533DEST_PATH_IMAGE030
individual electrical layer resistivity will suitably reduce a bit, i.e.
Figure 2012100797683100001DEST_PATH_IMAGE046
;
Figure 2012100797683100001DEST_PATH_IMAGE048
is the thickness correction factor, handles respectively according to the mutual relationship between
Figure 2012100797683100001DEST_PATH_IMAGE050
,
Figure 2012100797683100001DEST_PATH_IMAGE052
,
Figure 2012100797683100001DEST_PATH_IMAGE054
:
C) according to a) and b) to each electrically layer circulate, promptly circulates from 1 to N;
D) according to a), b) and c)
Figure 2012100797683100001DEST_PATH_IMAGE056
circulated; When the difference of the measured value of apparent resistivity and the theoretical fitting value of just drilling calculating when overall relative error
Figure 2012100797683100001DEST_PATH_IMAGE058
is less than preassigned number
Figure 2012100797683100001DEST_PATH_IMAGE060
on all periodic points, promptly
Figure 2012100797683100001DEST_PATH_IMAGE062
To stop circulation.
3. transient electromagnetic data according to claim 2 are degree of depth inversion method fast and accurately; It is characterized in that: the mutual relationship between described
Figure 498863DEST_PATH_IMAGE050
,
Figure 665533DEST_PATH_IMAGE052
,
Figure 401408DEST_PATH_IMAGE054
is for when
Figure 2012100797683100001DEST_PATH_IMAGE064
and ; if
Figure 397046DEST_PATH_IMAGE040
; ; Promptly the suitable attenuate of
Figure 383368DEST_PATH_IMAGE030
individual electrical layer thickness thickens and levels thickness is corresponding; if
Figure 970207DEST_PATH_IMAGE044
;
Figure 2012100797683100001DEST_PATH_IMAGE070
; Promptly
Figure 537586DEST_PATH_IMAGE030
individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness.
4. transient electromagnetic data according to claim 2 are degree of depth inversion method fast and accurately; It is characterized in that: the mutual relationship between described
Figure 602494DEST_PATH_IMAGE050
,
Figure 594721DEST_PATH_IMAGE052
,
Figure 867308DEST_PATH_IMAGE054
is for as
Figure 2012100797683100001DEST_PATH_IMAGE072
time; if
Figure 312196DEST_PATH_IMAGE040
; ; Then this layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down; if
Figure 259423DEST_PATH_IMAGE044
;
Figure 371736DEST_PATH_IMAGE074
; This layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer.
5. transient electromagnetic data according to claim 2 are degree of depth inversion method fast and accurately; It is characterized in that: the mutual relationship between described
Figure 300377DEST_PATH_IMAGE050
,
Figure 498141DEST_PATH_IMAGE052
, is for when
Figure 2012100797683100001DEST_PATH_IMAGE076
and
Figure 2012100797683100001DEST_PATH_IMAGE078
;
Figure 448834DEST_PATH_IMAGE040
;
Figure 564689DEST_PATH_IMAGE070
; Then
Figure 577644DEST_PATH_IMAGE030
individual electrical layer thickness suitably thickens, and the corresponding attenuate of levels thickness; if
Figure 928991DEST_PATH_IMAGE044
;
Figure 124218DEST_PATH_IMAGE068
; The suitable attenuate of
Figure 535607DEST_PATH_IMAGE030
individual electrical layer thickness thickens and levels thickness is corresponding.
6. transient electromagnetic data according to claim 2 are degree of depth inversion method fast and accurately; It is characterized in that: the mutual relationship between described
Figure 567017DEST_PATH_IMAGE050
,
Figure 456476DEST_PATH_IMAGE052
,
Figure 273253DEST_PATH_IMAGE054
is for as
Figure 2012100797683100001DEST_PATH_IMAGE080
time; if
Figure 449020DEST_PATH_IMAGE040
; ; Then this layer thickness is constant; But last layer bottom boundary buried depth all shoals a bit by identical amplitude respectively with interface, following one deck top buried depth, promptly moves on the layer; if
Figure 927721DEST_PATH_IMAGE044
;
Figure 458059DEST_PATH_IMAGE074
; This layer thickness is constant; But last layer bottom boundary buried depth all deepens a bit by identical amplitude respectively with interface, following one deck top buried depth, and promptly layer moves down.
7. according to the described transient electromagnetic data of above-mentioned any claim degree of depth inversion method fast and accurately;, it is characterized in that: the value of described resistivity correction factor
Figure 335886DEST_PATH_IMAGE036
, thickness correction factor , preassigned number
Figure 589461DEST_PATH_IMAGE060
is for survey the number of 3 above orders of magnitude on the value difference relatively.
8. transient electromagnetic data according to claim 7 are degree of depth inversion method fast and accurately;, it is characterized in that: described resistivity correction factor
Figure 239885DEST_PATH_IMAGE036
is 0.001.
9. transient electromagnetic data according to claim 7 are degree of depth inversion method fast and accurately;, it is characterized in that: described thickness correction factor
Figure 757454DEST_PATH_IMAGE048
is 0.001.
10. transient electromagnetic data according to claim 7 are degree of depth inversion method fast and accurately;, it is characterized in that: described preassigned number
Figure 126118DEST_PATH_IMAGE060
is 0.001.
CN2012100797683A 2012-03-23 2012-03-23 Fast and accurate depth inversion method of transient electromagnetic data Pending CN102608666A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100797683A CN102608666A (en) 2012-03-23 2012-03-23 Fast and accurate depth inversion method of transient electromagnetic data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100797683A CN102608666A (en) 2012-03-23 2012-03-23 Fast and accurate depth inversion method of transient electromagnetic data

Publications (1)

Publication Number Publication Date
CN102608666A true CN102608666A (en) 2012-07-25

Family

ID=46526163

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100797683A Pending CN102608666A (en) 2012-03-23 2012-03-23 Fast and accurate depth inversion method of transient electromagnetic data

Country Status (1)

Country Link
CN (1) CN102608666A (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103614494A (en) * 2013-11-12 2014-03-05 广东省农业科学院动物卫生研究所 Two-color fluorogenic quantitative PCR (Polymerase Chain Reaction) detection kit for CDV (canine distemper viruses) and CPV (canine parvo viruses)
CN103645512A (en) * 2013-12-10 2014-03-19 江苏大学 Multilayer ponding goaf detection method
CN104102814A (en) * 2014-06-11 2014-10-15 中国科学院地质与地球物理研究所 Magnetotelluric data based resistivity and magnetic susceptibility inversion method and system
CN104793254A (en) * 2015-04-28 2015-07-22 中国科学院地质与地球物理研究所 Underground electrical thin layer identification method and device
CN105089663A (en) * 2015-07-31 2015-11-25 中国石油天然气集团公司 High-angle well/horizontal well formation resistivity anisotropy correcting method
CN105277994A (en) * 2015-07-06 2016-01-27 中煤科工集团西安研究院有限公司 Determination method and device of differential conductance data synthetic aperture
CN104280782B (en) * 2013-07-12 2017-02-15 中国石油天然气集团公司 One-dimensional joint inversion method for time-frequency electromagnetic data and magnetotelluric data
CN112253090A (en) * 2020-10-14 2021-01-22 中海油田服务股份有限公司 Data parameter inversion method and device for multi-frequency electrical imaging
CN112346135A (en) * 2020-11-17 2021-02-09 黄河勘测规划设计研究院有限公司 Intelligent inspection detection method for leakage channel of canal embankment project based on Internet of things
CN113177330A (en) * 2021-05-27 2021-07-27 吉林大学 Transient electromagnetic rapid statistical inversion method
CN114859421A (en) * 2021-02-03 2022-08-05 中国科学院声学研究所 Underwater buried target identification method based on multi-parameter simultaneous inversion
CN115201922A (en) * 2022-06-28 2022-10-18 中国科学院地质与地球物理研究所 Semi-aeroelectromagnetic apparent resistivity measurement system and method
CN114859421B (en) * 2021-02-03 2024-05-31 中国科学院声学研究所 Underwater buried target identification method based on multi-parameter simultaneous inversion

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382599A (en) * 2007-09-03 2009-03-11 中国石油天然气集团公司 Transient electromagnetical method for reservoir pore space anisotropy
CN102053281A (en) * 2009-11-10 2011-05-11 中国石油化工集团公司 Method for realizing oil-gas detection by applying long-offset distance transient electromagnetic array method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382599A (en) * 2007-09-03 2009-03-11 中国石油天然气集团公司 Transient electromagnetical method for reservoir pore space anisotropy
CN102053281A (en) * 2009-11-10 2011-05-11 中国石油化工集团公司 Method for realizing oil-gas detection by applying long-offset distance transient electromagnetic array method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
《石油地球物理勘探》 20020415 苏朱刘等 大地电磁测深"正演修正法"一维反演 , 第02期 *
严良俊等: "长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法", 《石油地球物理勘探》 *
张跃彬: "磁法勘探在圈定煤田自燃区及在煤田地质勘察中的应用研究", 《中国新技术新产品》 *
苏朱刘等: "大地电磁测深"正演修正法"一维反演", 《石油地球物理勘探》 *
陈乐寿等: "几种简易有效的大地电磁一维反演方法", 《石油地球物理勘探》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104280782B (en) * 2013-07-12 2017-02-15 中国石油天然气集团公司 One-dimensional joint inversion method for time-frequency electromagnetic data and magnetotelluric data
CN103614494B (en) * 2013-11-12 2015-04-08 广东省农业科学院动物卫生研究所 Two-color fluorogenic quantitative PCR (Polymerase Chain Reaction) detection kit for CDV (canine distemper viruses) and CPV (canine parvo viruses)
CN103614494A (en) * 2013-11-12 2014-03-05 广东省农业科学院动物卫生研究所 Two-color fluorogenic quantitative PCR (Polymerase Chain Reaction) detection kit for CDV (canine distemper viruses) and CPV (canine parvo viruses)
CN103645512B (en) * 2013-12-10 2016-08-17 江苏大学 A kind of multilamellar hydrops Mined-Out Areas method
CN103645512A (en) * 2013-12-10 2014-03-19 江苏大学 Multilayer ponding goaf detection method
CN104102814A (en) * 2014-06-11 2014-10-15 中国科学院地质与地球物理研究所 Magnetotelluric data based resistivity and magnetic susceptibility inversion method and system
CN104102814B (en) * 2014-06-11 2017-07-11 中国科学院地质与地球物理研究所 A kind of method and system based on magnetotelluric data inverting resistivity and magnetic susceptibility
CN104793254A (en) * 2015-04-28 2015-07-22 中国科学院地质与地球物理研究所 Underground electrical thin layer identification method and device
CN105277994A (en) * 2015-07-06 2016-01-27 中煤科工集团西安研究院有限公司 Determination method and device of differential conductance data synthetic aperture
CN105089663A (en) * 2015-07-31 2015-11-25 中国石油天然气集团公司 High-angle well/horizontal well formation resistivity anisotropy correcting method
CN105089663B (en) * 2015-07-31 2017-11-07 中国石油天然气集团公司 A kind of High angle/horizontal well formation resistivity anisotropy bearing calibration
CN112253090A (en) * 2020-10-14 2021-01-22 中海油田服务股份有限公司 Data parameter inversion method and device for multi-frequency electrical imaging
CN112346135A (en) * 2020-11-17 2021-02-09 黄河勘测规划设计研究院有限公司 Intelligent inspection detection method for leakage channel of canal embankment project based on Internet of things
CN112346135B (en) * 2020-11-17 2023-04-14 黄河勘测规划设计研究院有限公司 Intelligent inspection detection method for leakage channel of canal embankment project based on Internet of things
CN114859421A (en) * 2021-02-03 2022-08-05 中国科学院声学研究所 Underwater buried target identification method based on multi-parameter simultaneous inversion
CN114859421B (en) * 2021-02-03 2024-05-31 中国科学院声学研究所 Underwater buried target identification method based on multi-parameter simultaneous inversion
CN113177330A (en) * 2021-05-27 2021-07-27 吉林大学 Transient electromagnetic rapid statistical inversion method
CN113177330B (en) * 2021-05-27 2022-07-22 吉林大学 Transient electromagnetic rapid statistical inversion method
CN115201922A (en) * 2022-06-28 2022-10-18 中国科学院地质与地球物理研究所 Semi-aeroelectromagnetic apparent resistivity measurement system and method

Similar Documents

Publication Publication Date Title
CN102608666A (en) Fast and accurate depth inversion method of transient electromagnetic data
Yin et al. A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling
Salem et al. Interpretation of magnetic data using an enhanced local wavenumber (ELW) method
Zhang et al. An improved and stable downward continuation of potential field data: The truncated Taylor series iterative downward continuation method
Cherevatova et al. Magnetotelluric array data analysis from north-west Fennoscandia
Feng et al. 3D gravity inversion of basement relief for a rift basin based on combined multinorm and normalized vertical derivative of the total horizontal derivative techniques
Sasaki et al. Frequency and time domain three-dimensional inversion of electromagnetic data for a grounded-wire source
Nguyen et al. Comparing large-scale 3D Gauss–Newton and BFGS CSEM inversions
Cevallos Automatic generation of 3D geophysical models using curvatures derived from airborne gravity gradient data
Bansal et al. Modelling of magnetic data for scaling geology
CN105259571A (en) Stratum inclination angle detection method
Song et al. Adaptive phase k-means algorithm for waveform classification
Datta et al. Full-waveform inversion of salt models using shape optimization and simulated annealing
Mulder et al. Time-domain modeling of electromagnetic diffusion with a frequency-domain code
Yu Using cylindrical surface-based curvature change rate to detect faults and fractures
Wang et al. Multispectral phase-based geosteering coherence attributes for deep stratigraphic feature characterization
Sajid et al. Logarithm of short‐time Fourier transform for extending the seismic bandwidth
Zhou et al. Improvement of ground-penetrating radar full-waveform inversion images using cone penetration test data
Rui‐Xue et al. The 3D transient electromagnetic forward modeling of volcanogenic massive sulfide ore deposits
Pardo et al. Fast and automatic inversion of LWD resistivity measurements for petrophysical interpretation
Zhang et al. Sand Bodies Delineation by Fusing Multifrequency Attributes via t-SNE
CN104991281B (en) A kind of detection method and device of coal seam buried depth
Bensdorp et al. An approximate 3D computational method for real-time computation of induction logging responses
Yu et al. An improved cylindrical surface fitting-related method for fault characterization
Zhang et al. Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120725