CN109541703A - A kind of quantitative preferred method of magnetotelluric sounding curve - Google Patents
A kind of quantitative preferred method of magnetotelluric sounding curve Download PDFInfo
- Publication number
- CN109541703A CN109541703A CN201811449764.3A CN201811449764A CN109541703A CN 109541703 A CN109541703 A CN 109541703A CN 201811449764 A CN201811449764 A CN 201811449764A CN 109541703 A CN109541703 A CN 109541703A
- Authority
- CN
- China
- Prior art keywords
- sounding
- discrete
- curve
- distance
- magnetic field
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The automatic preferred method of magnetotelluric sounding curve, belongs to field of geophysical exploration, to overcome artificial screening sounding curve, inefficiency and the strong disadvantage of subjectivity.The following steps are included: carrying out component data acquisition using N platform telluric electromagnetic sounding instrument, each depth measurement point is handled using different processing modes, obtain corresponding sounding curve, calculate the discrete Fu Leixie distance of apparent resistivity in sounding curve mutual discrete Fu Leixie distance and anti-phase between each other, if two sounding curves P and Q, its length is respectively the discrete Fu Leixie distance of m and n, then P and Q are as follows:Discrete Frechet distance result is weighted, the discrete Fu Leixie distance of local station He other acquisition station sounding curves is obtained；The above results reckling is screened, obtaining the sounding curve is most smooth reasonable sounding curve.The present invention is suitable for carrying out magneto-electrotelluric exploration.
Description
Technical field
The invention belongs to field of geophysical exploration, are related to magnetotelluric data processing technique, divide more particularly, to being directed to
The quantitative preferred method of magnetotelluric sounding curve after cloth data processing.
Background technique
Magnetotelluric sounding method (MT) is one-point or multi-point while to observe natural variation, electricity vertical each other on the ground
Magnetic-field component, to detect the electrical construction of earth interior.Since it has many advantages, such as that field construction is easy, investigation depth is big,
Therefore this method study Deep Geological Structures, mineral exploration, in terms of be widely applied.
Currently, telluric electromagnetic sounding is gradually developed to distribution observation from traditional Single Point Surveying, to improve observation effect
Rate, meanwhile, the magnetotelluric data after distribution observation is convenient for carrying out distributed treatment using a variety of tensor evaluation methods.It compares
Single Point Surveying is only capable of one group of sounding curve of processing, and distribution observation can be carried out from the multiple groups sounding curve after distributed treatment
It is preferred that.
There are many tensor evaluation methods at present to carry out magnetotelluric sounding curve estimation, comprising: standalone processes (Sims,
Etal., 1971), using local station related and noise non-correlation to reference station signal remote reference process (GAMBLE, 1979,
Referring specifically to: GAMBLE T D.magnetotellurics with a remote reference [J] .Geophysics,
1979,81 (1): 87-94) (hereinafter referred to as remote reference process), utilize the magnetic field phase of local station and reference station magnetic field correlation
Close remote reference process (Varentsov, 2003, referring specifically to VARENTSOV I M.Chapter 10Arrays
ofSimultaneous Electromagnetic Soundings:Design,Data Processing and Analysis
[J].Methods in Geochemistry&Geophysics,2003,40:259-273；Zhang Gang, etc., 2017, referring specifically to:
Zhang Gang, TuoXian state, Wang Xuben wait the magnetic field correlation far with reference to application [J] petroleum earth in magnetotelluric data processing
Physical prospecting, 2017 (6): 1333-1343) (hereinafter referred to as based on the related remote reference process in magnetic field), shared track processing (
Wen Xiu, etc., 2012, referring specifically to: Wenxiu ZHANG, Zhou Fengdao, Lin Jun wait distributed electrical Magnetic Detection System to provide in Deep Groundwater
Source reconnoitre in application [J] Jilin University journal (geoscience version), 2012,42 (4): 1207-1213.).For above-mentioned four
Kind data tensor estimation, if having N platform telluric electromagnetic sounding instrument in same survey area while carrying out distributed observation, for a certain sight
For measuring point, standalone processes mode has 1 group of processing result using electric field data and the magnetic field data processing of local station；Remote reference portion
Reason, based on the related remote reference process in magnetic field using other N-1 acquisition stations magnetic field data and local station electromagnetic field data into
Row processing has N-1 group processing result respectively；Shared track processing utilizes the magnetic field data of this earth electric field and other N-1 acquisition station
It is handled, also there is a N-1 group processing result；So all processing result sums of the acquisition station are 3N-2 group, it is assumed that distribution
Formula observation shares 10 acquisition stations, then up to 28 groups of processing results.The depth measurement that above-mentioned four kinds of data processing methods are calculated is bent
Knot fruit is not fully consistent.In general, the sounding curve of remote reference process will be got well compared to standalone processes effect, premise item
Part be it is more demanding to the quality of data of reference station, need to be laid in the environment of weak/noiseless interference, but generally exist
Carry out distributed capture when, the quality of data of each acquisition station cannot be estimated in advance, it is possible that will appear remote reference process it
The case where sounding curve afterwards is not as good as standalone processes, same situation may occur based on the related remote reference process in magnetic field
In shared track processing method.So how quantitative assessment particular acquisition station four kinds of results calculated which is better and which is worse, sieve
Best in quality, most reasonable sounding curve is selected, is the important guarantee for obtaining local reliable electrical conductivity structure model.
It is to carry out preferably, having in the way of artificial screening in the prior art to filter out reliable sounding curve
Certain subjectivity and low efficiency, and how quantitative screening goes out to characterize the sounding curve needs of underground electrical structure studies
Problem.
Summary of the invention
The technical problem to be solved by the present invention is to the disadvantages mentioned above in order to overcome the prior art, provide a kind of magnetotelluric
The quantitative preferred method of sounding curve is not necessarily to artificial screening, can quantitative screening go out to characterize underground electrical structure depth measurement it is bent
Line improves treatment effeciency, enhances the accuracy of preferred objectivity and result.
The technical solution adopted by the present invention to solve the technical problems is: a kind of magnetotelluric sounding curve quantifies preferably
Method carries out component data acquisition to N number of depth measurement point respectively using N platform telluric electromagnetic sounding instrument, comprising the following steps:
A. each depth measurement point by conventional single point tensor impedance estimate processing mode, in the way of remote reference process, be based on magnetic field
Related remote reference process mode, magnetic field share at least two progress distributed data processings in processing mode, obtain every kind of place
Sounding curve after reason mode, and record the frequency points of each group of sounding curve；
B. apparent resistivity in any two groups of sounding curves mutual discrete Fu Leixie distance and anti-phase phase are calculated
Discrete Fu Leixie distance between mutually, discrete Fu Leixie distance calculate according to the following formula: setting two sounding curves P and Q, frequency point
Number is respectively m and n, then the discrete Fu Leixie distance of P and Q are as follows:
Wherein P:[0:m] → V, P is a mapping from set [0:m]={ 0,1 ..., m } to SPACE V, α: [1:m+n]
→ [0, m], α are a mapping from set [1:m+n]={ 1,2 ..., m+n } to set [0:m]={ 0,1 ..., m }, α
(k), k=1,2 ... m+n indicates mapping value of the mapping α when independent variable takes k；Q:[0:n] → V, Q be from set [0:n]=
{ 0,1 ..., n } to a mapping of SPACE V, β: [1:m+n] → [0, n], β is from set [1:m+n]={ 1,2 ..., m+n }
It is mapped to set [0:n]={ 0,1 ..., n } one, β (k), k=1,2 ... m+n indicates mapping β when independent variable takes k
Mapping value,Indicate the compound operation of operator, | | | |2Indicate Euclidean distance；ψm,nBe by two it is continuous, do not subtract and surjection
The all possible combinations that operator is constituted,
ψm,n=Mon ([1:m+n], [0:m]) × Mon ([1:m+n], [0:n]) (2)
In formula (2), Mon (X, Y) be by it is continuous in set X to set Y, do not subtract and the operator of surjection is constituted, [1:m+n]
={ 1,2 ..., m+n }, [0:m]={ 0,1 ..., m }, [0:n]={ 0,1 ..., n }；
Local station B and other acquisition stations R are calculated according to above-mentioned calculating processjApparent resistivity Fu Leixie distance(B,Rj)
With the Fu Leixie distance of impedance phase(B,Rj), wherein i indicates the direction xy yx, j=1 ..., N, and N is sounding curve group
Number；
C. discrete Frechet distance result is weighted, obtain local station and other acquisition station sounding curves it is discrete not
Thunder is had a rest distance:
Wherein, i indicates the direction xy yx, j=1 ..., N, and N is sounding curve group number, and l and d are weighted factor；
D. δ is filtered outF(B,Rj) reckling, obtaining the sounding curve is most smooth reasonable sounding curve.
Specifically, the component data acquisition includes that the acquisition of 5 components is acquired with 4 components, wherein 5 components include that electric field 2 divides
Amount, 3 component of magnetic field；4 components include 2 component of electric field, 2 component of magnetic field.
Specifically, electric field data and magnetic field of the conventional single point tensor impedance estimation processing mode using local station in step A
Data are handled to obtain 1 group of processing result, and remote reference process is adopted based on the related remote reference process in magnetic field using other N-1
The electromagnetic field data of the magnetic field data and local station that collect station carries out processing and respectively obtains N-1 group processing result；Magnetic field shares processing
Mode is handled using the magnetic field data of local electric field data and other N-1 acquisition station, obtains N-1 group processing result.
Preferably, the range of l and d is respectively [0,1].
The beneficial effects of the present invention are: quantified using discrete Fu Leixie distance, automatic screening sounding curve, it is existing to overcome
When having artificial screening sounding curve in technology subjectivity it is strong, it is time-consuming long, inefficiency, and the shortcomings that without quantitative target, it is specific and
Yan Shi, using discrete Fu Leixie distance come the multiple groups sounding curve after quantitative assessment distributed data processing, discrete Fu Leixie away from
It is formed by two curves with a distance from all discrete points from the path using two targets, measures the similarity of two curves,
The optimal result of sounding curve after then filtering out distributed treatment；This method saves human cost, and filtering out can characterize ground
The sounding curve of lower electrical structure improves result objectivity and accuracy, improves treatment effeciency, and use scope is wider.The present invention
Suitable for screening sounding curve.
Detailed description of the invention
Fig. 1 is flow diagram of the invention；
Fig. 2 is that forward modeling sounding curve adds Gauss, Rayleigh, F distribution, logarithm normal distribution, Poisson 5 respectively
Result after kind noise；
Fig. 3 is 7 groups of sounding curves after the distributed data processing of L39 measuring point.
Specific embodiment
With reference to the accompanying drawings and embodiments, technical scheme in the embodiment of the invention is clearly and completely described, shows
So, described embodiment is only a part of the embodiments of the present invention, instead of all the embodiments.Below at least one
The description only actually of exemplary embodiment be it is illustrative, never as to the present invention and its application or any limit used
System.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts
Every other embodiment, shall fall within the protection scope of the present invention.Non-elaborated part of the present invention belongs to techniques well known.
Technology, method and apparatus known to person of ordinary skill in the relevant may be not discussed in detail, but suitable
In the case of, the technology, method and apparatus should be considered as authorizing part of specification.
As shown in Figure 1, the quantitative preferred method of magnetotelluric sounding curve of the invention, comprising the following steps:
(1) component data acquisition is carried out to N number of depth measurement point respectively using N platform telluric electromagnetic sounding instrument, it is commonly used in the art
Be the acquisition of 5 components or the acquisition of 4 components, i.e. component acquisition includes that the acquisition of 5 components is acquired with 4 components, wherein 5 component packets
Include 2 component of electric field, 3 component of magnetic field；4 components include 2 component of electric field, 2 component of magnetic field；
(2) each depth measurement point by conventional single point tensor impedance estimate processing mode, in the way of remote reference process, be based on magnetic field
Related remote reference process mode, magnetic field share at least two progress distributed data processings in processing mode, obtain every kind of place
Sounding curve after reason mode, and record the frequency points of each group of sounding curve；
Conventional single point tensor impedance estimates that processing mode is handled using the electric field data and magnetic field data of local station
Obtain 1 group of processing result, remote reference process, the magnetic field number that other N-1 acquisition stations are utilized based on the related remote reference process in magnetic field
Processing, which is carried out, according to the electromagnetic field data with local station respectively obtains N-1 group processing result；Magnetic field shares processing mode using locally
The magnetic field data of electric field data and other N-1 acquisition station is handled, and obtains N-1 group processing result.
In the examples below, in order to select optimal sounding curve, be all handled with four kinds of processing modes, and
After carry out screening comparison.In actual operation, the several processing mode greater than 1 can arbitrarily be selected to be screened.Then same
When selecting above-mentioned four kinds of processing modes, 3N-2 group sounding curve can be obtained altogether.
(3) the discrete Fu Leixie distance that apparent resistivity is mutual in above-mentioned any two groups of sounding curves is calculated, and each
The mutual discrete Fu Leixie distance of group impedance phase.When using four kinds of processing modes, that is, calculate 3N-2 group sounding curve
The discrete Fu Leixie distance of middle relevant parameter.
Discrete Fu Leixie distance calculates according to the following formula: setting two sounding curves P and Q, frequency points are respectively m and n, then
The discrete Fu Leixie distance (Sriraghavendra et al, 2007) of P and Q are as follows:
Wherein P:[0:m] → V, P is a mapping from set [0:m]={ 0,1 ..., m } to SPACE V, α: [1:m+n]
→ [0, m], α are a mappings from set [1:m+n]={ 1,2 ..., m+n } to set [0:m]={ 0,1 ..., m }, α
(k), k=1,2 ... m+n indicates mapping value of the mapping α when independent variable takes k；Q:[0:n] → V, Q be from set [0:n]=
{ 0,1 ..., n } to a mapping of SPACE V, β: [1:m+n] → [0, n], β is from set [1:m+n]={ 1,2 ..., m+n }
It is mapped to set [0:n]={ 0,1 ..., n } one, β (k), k=1,2 ... m+n indicates mapping β when independent variable takes k
Mapping value,Indicate the compound operation of operator, | | | |2Indicate Euclidean distance；
ψm,n=Mon ([1:m+n], [0:m]) × Mon ([1:m+n], [0:n]) (2)
ψm,nBe by two it is continuous, do not subtract and all possible combinations that Surjectivity Operators are constituted, such as α1,α2,α3It is Mon ([1:m+
N], [0:m]) in element, β1,β2It is the element in Mon ([1:m+n], [0:n]), then ψm,n={ (α1,β1),(α1,β2),
(α2,β1),(α2,β2),(α3,β1),(α3,β2)}.In formula (2), Mon (X, Y) be by it is continuous in set X to set Y, do not subtract and full
The operator penetrated is constituted, [1:m+n]={ 1,2 ..., m+n }, indicates the set of integers from 1 to m+n, and 1≤m+n；[0:m]=0,
1 ..., m }, indicate the set of integers from 0 to m, and 0≤m；[0:n]={ 0,1 ..., n }, set of integers of the expression from 0 to n, and 0≤
N, m and n are respectively the frequency points of sounding curve P and Q, according to different magnetotelluric methods, general value 40-100, also and frequently
The density degree of point is related.
Local station B and other acquisition stations R are calculated according to above-mentioned calculating processjApparent resistivity Fu Leixie distance(B,Rj)
With the Fu Leixie distance of impedance phase(B,Rj), wherein i indicates the direction xy yx, and j=1 ..., N, N is sounding curve
Group number.
(4) discrete Frechet distance result is weighted, obtains the discrete of local station and other acquisition station sounding curves
Fu Leixie distance:
Wherein, i indicates the direction xy yx, j=1 ..., N, and N is sounding curve group number, and l and d are weighted factor, l and d
Range be [0,1], based on optimize consider, take l=0.6, d=0.4
(5) δ is filtered outF(B,Rj) reckling, obtaining this group of sounding curve is most smooth reasonable sounding curve, is obtained
Optimal sounding curve and output.
This method is based on discrete Fu Leixie distance, and the multiple groups sounding curve after quantitative screening distributed data processing avoids
Utilize the subjectivity and time-consuming and laborious drawback of artificial screening.
Embodiment 1
The correctness of method is mainly described in theoretical modeling, and optimal curve can be filtered out using the method.Generally,
Field of geophysical exploration has a new method to propose, it is necessary first to carry out theoretical model verifying, then be generalized to again in practice
It goes.
This example is theoretical modeling calculating, is all handled with four kinds of processing modes, bent using above-mentioned telluric electromagnetic sounding
The quantitative preferred method of line includes the following steps:
(1) establish stratified model: in the case where 0-100km depth resistivity is the background of 1000 Ω m, 5-20km depth is 10
The low-resistance of Ω m is abnormal, and the sounding curve after obtaining forward modeling is as shown in Fig. 2 solid black lines.
(2) in order to simulate the exception of sounding curve caused by various interference, 5 kinds of noises are added in original sounding curve, are divided
It is not: Gaussian Profile (Gauss) noise, rayleigh distributed (Rayleigh) noise, F partition noise, logarithm normal distribution noise, pool
Loose (Poisson) noise, obtains multiple groups sounding curve as shown in discrete point in Fig. 2；
(3) the discrete Fu Leixie distance that apparent resistivity is mutual in above-mentioned 6 groups of sounding curves is calculated(B,Rj), and
The mutual discrete Fu Leixie distance of 6 groups of sounding curve middle impedance phases(B,Rj), calculate discrete Fu Leixie distance such as public affairs
Shown in formula (1)；
(4) discrete Fu Leixie distance is weighted, weight computation process such as formula (3), obtained weighted results such as table 1
It is shown；
1 forward modeling sounding curve of table with plus after making an uproar every group of sounding curve Fu Leixie distance weights
(5) the discrete Fu Leixie distance δ of weighting is filtered outF(B,Rj) in reckling, obtain this group of sounding curve be most light
Sliding reasonable sounding curve.It is computed, sounding curve and the discrete Fu Leixie of other groups distance are 26.54 after original forward modeling, are
Minimum, thus after selecting forward modeling in 6 groups of sounding curves sounding curve as final result.
The calculated results demonstrate the correctness using discrete Fu Leixie apart from automatic screening sounding curve.
Embodiment 2
This example is practical field data, and a kind of quantitative preferred method of magnetotelluric sounding curve, L33, L36, L39 are position
Three telluric electromagnetic sounding points of observation are synchronized in across Longmenshan Structural Belt south west section, in order to from L39 measuring point distribution
One group of optimal sounding curve is filtered out in reason result, is included the following steps:
(1) for L39 measuring point, distributed data processing is carried out using four kinds of processing modes, specifically, conventional single point
Amount impedance estimation processing (L39SS) handles out one group of sounding curve, carries out remote reference process to L39 measuring point using L33, L36
(L39vsL33_RR, L39vsL36_RR) handle out 2 groups of sounding curves, the related remote reference process in magnetic field (L39vsL33_RRHC,
L39vsL36_RRHC 2 groups of sounding curves) are handled out, magnetic field shares processing (L39vsL33_UH, L39vsL36_UH) and handles out 2
Group sounding curve, is obtained 7 groups of sounding curves, as shown in Figure 3；
(2) the discrete Fu Leixie distance and 7 groups of impedances that apparent resistivity is mutual in above-mentioned 7 groups of sounding curves are calculated
The mutual discrete Fu Leixie distance of phase.Discrete Fu Leixie distance is calculated according to (1) formula；
(3) discrete Fu Leixie distance is weighted, weight computation process such as formula (3), obtained weighted results such as table 2
It is shown；
The weighted value of every group of Fu Leixie distance after the distributed data processing of table 2L39 and L33, L36
(4) the discrete Fu Leixie distance δ of weighting is filtered outF(B,Rj) in reckling, this group of sounding curve is most smooth conjunction
The sounding curve of reason.It is computed, in 7 groups of sounding curves, remote reference process L39vsL33_RR is discrete with other group of sounding curve
Fu Leixie distance is 122.74, and for minimum, the result after thus selecting remote reference process is final result.
Use scope of the present invention is wider, utilizes the audio-frequency magnetotelluric magnetic method (AMT) of Natural electromagnetic field exploration, wideband the earth electricity
Magnetic method (MT) and long period magnetotelluric method (LMT) all can be used.
Claims (4)
1. a kind of quantitative preferred method of magnetotelluric sounding curve, which is characterized in that distinguished using N platform telluric electromagnetic sounding instrument
Component data acquisition is carried out to N number of depth measurement point, comprising the following steps:
A. each depth measurement point by conventional single point tensor impedance estimate processing mode, in the way of remote reference process, it is related based on magnetic field
Remote reference process mode, magnetic field share at least two progress distributed data processings in processing mode, obtain every kind of processing side
Sounding curve after formula, and record the frequency points of each group of sounding curve；
B. calculate the mutual discrete Fu Leixie distance of apparent resistivity in any two groups of sounding curves and anti-phase mutually it
Between discrete Fu Leixie distance,
Discrete Fu Leixie distance calculates according to the following formula: setting two sounding curves P and Q, frequency points are respectively m and n, then P and Q
Discrete Fu Leixie distance are as follows:
Wherein P:[0:m] → V, P is that one from set [0:m]={ 0,1 ..., m } to SPACE V maps, α: [1:m+n] →
[0, m], α are mapped from set [1:m+n]={ 1,2 ..., m+n } to one of set [0:m]={ 0,1 ..., m }, α (k),
K=1,2 ... m+n indicates mapping value of the mapping α when independent variable takes k；Q:[0:n] → V, Q be from set [0:n]=0,
1 ..., n } to a mapping of SPACE V, β: [1:m+n] → [0, n], β is arrived from set [1:m+n]={ 1,2 ..., m+n }
One for gathering [0:n]={ 0,1 ..., n } maps, β (k), k=1, and 2 ... m+n indicates mapping β reflecting when independent variable takes k
Value is penetrated,Indicate the compound operation of operator, | | | |2Indicate Euclidean distance；
ψm,n=Mon ([1:m+n], [0:m]) × Mon ([1:m+n], [0:n]) (2)
ψm,nBe by two it is continuous, do not subtract and all possible combinations that Surjectivity Operators are constituted, in formula (2), Mon (X, Y) is by gathering
It is continuous in X to set Y, do not subtract and the operator of surjection is constituted, [1:m+n]={ 1,2 ..., m+n }, [0:m]=0,1 ...,
M }, [0:n]={ 0,1 ..., n }；
Local station B and other acquisition stations R are calculated according to above-mentioned calculating processjApparent resistivity Fu Leixie distanceAnd resistance
The Fu Leixie distance of anti-phaseWherein, i indicates the direction xy yx, j=1 ..., N, and N is sounding curve group number；
C. discrete Frechet distance result is weighted, obtains the discrete Fu Leixie of local station He other acquisition station sounding curves
Distance:
Wherein, i indicates the direction xy yx, j=1 ..., N, and N is sounding curve group number, and l and d are weighted factor；
D. δ is filtered outF(B,Rj) reckling to get arrive optimal sounding curve and output.
2. the method as described in claim 1, which is characterized in that the component data acquisition includes that the acquisition of 5 components is adopted with 4 components
Collection, wherein 5 components include 2 component of electric field, 3 component of magnetic field；4 components include 2 component of electric field, 2 component of magnetic field.
3. the method as described in claim 1, which is characterized in that conventional single point tensor impedance estimation processing mode benefit in step A
It is handled to obtain 1 group of processing result with the electric field data of local station and magnetic field data, remote reference process is based on magnetic field correlation far
Reference process carries out processing using the magnetic field data of other N-1 acquisition stations and the electromagnetic field data of local station and respectively obtains N-1
Group processing result；Magnetic field is shared processing mode and is carried out using the magnetic field data of local electric field data and other N-1 acquisition station
Processing, obtains N-1 group processing result.
4. the method as described in claim 1, which is characterized in that the range of l and d is respectively [0,1].
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811449764.3A CN109541703A (en) | 2018-11-30 | 2018-11-30 | A kind of quantitative preferred method of magnetotelluric sounding curve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811449764.3A CN109541703A (en) | 2018-11-30 | 2018-11-30 | A kind of quantitative preferred method of magnetotelluric sounding curve |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109541703A true CN109541703A (en) | 2019-03-29 |
Family
ID=65851518
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811449764.3A Pending CN109541703A (en) | 2018-11-30 | 2018-11-30 | A kind of quantitative preferred method of magnetotelluric sounding curve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109541703A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120261135A1 (en) * | 2011-04-15 | 2012-10-18 | Nowak Ethan J | Shape-Based Metrics In Reservoir Characterization |
| US20140129194A1 (en) * | 2012-11-02 | 2014-05-08 | Technolmaging, Llc. | Methods of three-dimensional potential field modeling and inversion for layered earth models |
| CN104375195A (en) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | Time-frequency electromagnetic multi-source multi-component three-dimensional joint inversion method |
| CN105866852A (en) * | 2016-04-01 | 2016-08-17 | 西南科技大学 | Remote reference magnetotelluric impedance calculation method based on correlation detection |
-
2018
- 2018-11-30 CN CN201811449764.3A patent/CN109541703A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120261135A1 (en) * | 2011-04-15 | 2012-10-18 | Nowak Ethan J | Shape-Based Metrics In Reservoir Characterization |
| US20140129194A1 (en) * | 2012-11-02 | 2014-05-08 | Technolmaging, Llc. | Methods of three-dimensional potential field modeling and inversion for layered earth models |
| CN104375195A (en) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | Time-frequency electromagnetic multi-source multi-component three-dimensional joint inversion method |
| CN105866852A (en) * | 2016-04-01 | 2016-08-17 | 西南科技大学 | Remote reference magnetotelluric impedance calculation method based on correlation detection |
Non-Patent Citations (4)
| Title |
|---|
| R. SRIRAGHAVENDRA，等: "Frechet Distance based Approach for Searching Online Handwritten Documents", 《IEEE》 * |
| T. D. GAMBLE ，等: "Magnetotellurics with a remote magnetic reference", 《GEOPHYSICS》 * |
| 张刚，等: "分布式大地电磁测深曲线的定量优选方法研究", 《中国地球科学联合学术年会2018》 * |
| 张刚，等: "磁场相关性在远参考大地电磁数据处理中的应用", 《石油地球物理勘探》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101556340B (en) | Three-dimensional small surface element magnetotelluric continuous array data acquisition method | |
| Labson et al. | Geophysical exploration with audiofrequency natural magnetic fields | |
| Nippress et al. | Optimized automatic pickers: application to the ANCORP data set | |
| CN105116453B (en) | Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone | |
| CN106022339B (en) | A kind of extracting method of Reclaimed Land shallow embedding underground pipe depth | |
| CN103809204B (en) | A kind of collecting method of field audio magnetotelluric method | |
| CN103210323B (en) | Process the system and method for geological data | |
| Küperkoch et al. | Automated determination of S-phase arrival times using autoregressive prediction: application to local and regional distances | |
| WO2008121770A2 (en) | Csem surveying method using an arrangement of plural signal sources | |
| US8706462B2 (en) | System and method for providing a physical property model | |
| CN109541695A (en) | Artificial field source frequency domain electric-force gradient far-zone apparent resistivity fast imaging method | |
| CN110007357A (en) | A kind of aviation TEM and aviation MT joint inversion method | |
| CN104570076A (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| Beggan et al. | Geolectric field measurement, modelling and validation during geomagnetic storms in the UK | |
| CN109061752A (en) | A kind of resistivity curve bearing calibration on the stratum containing grey matter | |
| CN108169802B (en) | A kind of slow diffusion imaging method of time domain electromagnetic data of harsh media model | |
| CN108594301B (en) | A kind of method and processing terminal of the seismic data fusion with difference characteristic | |
| CN109782355A (en) | The detection method and device of OBS detection point drift | |
| Acworth | The electrical image method compared with resistivity sounding and electromagnetic profiling for investigation in areas of complex geology: A case study from groundwater investigation in a weathered crystalline rock environment | |
| CN109541703A (en) | A kind of quantitative preferred method of magnetotelluric sounding curve | |
| Abraham et al. | Use of a handheld broadband EM induction system for deriving resistivity depth images | |
| CN110109184B (en) | Passive field source type three-dimensional electric field exploration method based on multiple daily change points | |
| CN108363739A (en) | A kind of seismic data low-and high-frequency expanding method based on sparse acquisition | |
| CN107329183A (en) | A kind of controlled-source audiomagnetotellurics sounding collecting method and device | |
| Soliman et al. | Evaluation of Saudi National Seismic Network (SNSN) detectability |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190329 |