CN109541703A - A kind of quantitative preferred method of magnetotelluric sounding curve - Google Patents

A kind of quantitative preferred method of magnetotelluric sounding curve Download PDF

Info

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
magnetic field
curve
distance
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
CN201811449764.3A
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.)
Southwest University of Science and Technology
Original Assignee
Southwest University of Science and Technology
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 Southwest University of Science and Technology filed Critical Southwest University of Science and Technology
Priority to CN201811449764.3A priority Critical patent/CN109541703A/en
Publication of CN109541703A publication Critical patent/CN109541703A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing 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

大地电磁测深曲线的自动优选方法,属于地球物理勘探领域,以克服人工筛选测深曲线,效率低下且主观性强的缺点。包括以下步骤:采用N台大地电磁测深仪进行分量数据采集,每个测深点利用不同的处理方式进行处理,得到相应测深曲线,计算测深曲线中视电阻率相互之间的离散弗雷歇距离以及抗相位相互之间的离散弗雷歇距离,设两条测深曲线P和Q,其长度分别为m和n,则P和Q的离散弗雷歇距离为:将离散弗雷歇距离计算结果加权,得到本地站和其他采集站测深曲线的离散弗雷歇距离;筛选上述结果最小者,得到该测深曲线即为最光滑合理的测深曲线。本发明适用于进行大地电磁勘探。

The automatic optimization method of the magnetotelluric sounding curve belongs to the field of geophysical exploration to overcome the shortcomings of manual screening of the sounding curve, low efficiency and strong subjectivity. It includes the following steps: using N sets of magnetotelluric sounders to collect component data, each sounding point is processed by different processing methods to obtain a corresponding sounding curve, and calculating the discrete Frey difference between apparent resistivity in the sounding curve. The discrete Frescher distance and the discrete Frescher distance between the anti-phases, suppose two sounding curves P and Q, whose lengths are m and n respectively, then the discrete Frescher distance of P and Q is: The discrete Frechet distance calculation results are weighted to obtain the discrete Frechet distances of the sounding curves of the local station and other acquisition stations; the one with the smallest results is selected, and the sounding curve is the most smooth and reasonable sounding curve. The invention is suitable for conducting magnetotelluric exploration.

Description

A kind of quantitative preferred method of magnetotelluric sounding curve
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 α123It is Mon ([1:m+ N], [0:m]) in element, β12It is the element in Mon ([1:m+n], [0:n]), then ψm,n={ (α11),(α12), (α21),(α22),(α31),(α32)}.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.一种大地电磁测深曲线的定量优选方法,其特征在于,采用N台大地电磁测深仪分别对N个测深点进行分量数据采集,包括以下步骤:1. a quantitative optimization method of magnetotelluric sounding curve, is characterized in that, adopts N magnetotelluric sounders to carry out component data acquisition to N sounding points respectively, comprises the following steps: A.每个测深点利用常规单点张量阻抗估算处理方式、远参考处理方式、基于磁场相关远参考处理方式、磁场共用处理方式中的至少两种进行分布式数据处理,得到每种处理方式后的测深曲线,并记录每一组测深曲线的频点个数;A. Each sounding point uses at least two of the conventional single-point tensor impedance estimation processing method, the remote reference processing method, the magnetic field correlation-based remote reference processing method, and the magnetic field sharing processing method to perform distributed data processing, and obtain each processing method. The sounding curve after the method, and record the number of frequency points of each group of sounding curves; B.计算任意两组测深曲线中视电阻率相互之间的离散弗雷歇距离以及抗相位相互之间的离散弗雷歇距离,B. Calculate the discrete Frecher distance between the apparent resistivity and the anti-phase between any two sets of sounding curves, 离散弗雷歇距离按照下式计算:设两条测深曲线P和Q,其频点个数分别为m和n,则P和Q的离散弗雷歇距离为:The discrete Frechet distance is calculated according to the following formula: Suppose two sounding curves P and Q, whose frequency points are m and n respectively, then the discrete Frechet distance of P and Q is: 其中P:[0:m]→V,P是从集合[0:m]={0,1,...,m}到空间V的一个映射,α:[1:m+n]→[0,m],α是从集合[1:m+n]={1,2,...,m+n}到集合[0:m]={0,1,...,m}的一个映射,α(k),k=1,2,...m+n表示映射α在自变量取k时的映射值;Q:[0:n]→V,Q是从集合[0:n]={0,1,...,n}到空间V的一个映射,β:[1:m+n]→[0,n],β是从集合[1:m+n]={1,2,...,m+n}到集合[0:n]={0,1,...,n}的一个映射,β(k),k=1,2,...m+n表示映射β在自变量取k时的映射值,表示算子的复合运算,||·||2表示欧几里得距离;where P:[0:m]→V, P is a mapping from the set [0:m]={0,1,...,m} to the space V, α:[1:m+n]→[ 0,m], α is from set [1:m+n]={1,2,...,m+n} to set[0:m]={0,1,...,m} A mapping, α(k), k=1,2,...m+n represents the mapping value of mapping α when the independent variable takes k; Q:[0:n]→V, Q is from the set [0: n]={0,1,...,n} to a mapping of space V, β:[1:m+n]→[0,n], β is from the set [1:m+n]={ A mapping from 1,2,...,m+n} to the set [0:n]={0,1,...,n}, β(k),k=1,2,...m +n represents the mapping value of the mapping β when the independent variable takes k, Represents the composite operation of the operator, ||·|| 2 represents the Euclidean distance; ψm,n=Mon([1:m+n],[0:m])×Mon([1:m+n],[0:n]) (2)ψ m,n =Mon([1:m+n],[0:m])×Mon([1:m+n],[0:n]) (2) ψm,n是由两个连续、不减且满射算子构成的所有可能组合,式(2)中,Mon(X,Y)是由集合X到集合Y中连续、不减且满射的算子所构成,[1:m+n]={1,2,…,m+n},[0:m]={0,1,…,m},[0:n]={0,1,…,n};ψ m,n is all possible combinations composed of two continuous, non-decreasing and surjective operators. In formula (2), Mon(X, Y) is a continuous, non-decreasing and surjective from set X to set Y composed of operators, [1:m+n]={1,2,…,m+n}, [0:m]={0,1,…,m}, [0:n]={0 ,1,…,n}; 按照上述计算过程计算本地站B和其他采集站Rj的视电阻率弗雷歇距离和阻抗相位的弗雷歇距离其中,i表示xy或者yx方向,j=1,…,N,N为测深曲线组数;Calculate the apparent resistivity Frechet distance of the local station B and other acquisition stations R j according to the above calculation process Frecher distance and impedance phase Among them, i represents the xy or yx direction, j=1,...,N, N is the number of sounding curve groups; C.将离散弗雷歇距离计算结果加权,得到本地站和其他采集站测深曲线的离散弗雷歇距离:C. Weighting the discrete Frechet distance calculation results to obtain the discrete Frechet distances of the local station and other acquisition stations' bathymetric curves: 其中,i表示xy或者yx方向,j=1,…,N,N为测深曲线组数,l和d为加权因子;Among them, i represents the xy or yx direction, j=1,...,N, N is the number of sounding curve groups, and l and d are weighting factors; D.筛选出δF(B,Rj)最小者,即得到最优的测深曲线并输出。D. Screen out the one with the smallest δ F (B, R j ), that is, get the optimal sounding curve and output it. 2.如权利要求1所述的方法,其特征在于,所述分量数据采集包括5分量采集与4分量采集,其中,5分量包括电场2分量,磁场3分量;4分量包括电场2分量,磁场2分量。2 . The method according to claim 1 , wherein the component data acquisition includes 5-component acquisition and 4-component acquisition, wherein the 5-component includes 2 components of the electric field and 3 components of the magnetic field; the 4-component includes 2 components of the electric field and the magnetic field 2 servings. 3.如权利要求1所述的方法,其特征在于,步骤A中常规单点张量阻抗估算处理方式利用本地站的电场数据和磁场数据进行处理得到1组处理结果,远参考处理、基于磁场相关远参考处理利用其它N-1个采集站的磁场数据和本地站的电磁场数据进行处理分别得到N-1组处理结果;磁场共用处理方式利用本地电场数据和其它N-1个采集站的磁场数据来进行处理,得到N-1组处理结果。3. The method of claim 1, wherein in step A, the conventional single-point tensor impedance estimation processing method utilizes the electric field data and magnetic field data of the local station to process to obtain 1 group of processing results, remote reference processing, based on the magnetic field. The relative remote reference processing uses the magnetic field data of other N-1 acquisition stations and the electromagnetic field data of the local station to process N-1 sets of processing results respectively; the magnetic field sharing processing method uses the local electric field data and the magnetic field of other N-1 acquisition stations. The data are processed to obtain N-1 sets of processing results. 4.如权利要求1所述的方法,其特征在于,l和d的范围分别为[0,1]。4. The method of claim 1, wherein the ranges of l and d are [0, 1], respectively.
CN201811449764.3A 2018-11-30 2018-11-30 A kind of quantitative preferred method of magnetotelluric sounding curve Pending CN109541703A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
US7519474B2 (en) Method and apparatus for measuring the resistivity of electromagnetic waves of the earth
CA2921340C (en) Automated horizon auto-picking on multiple volumes
CN112083508B (en) A denoising method and system for artificial source electromagnetic exploration signal with noise reference track
Kaufmann et al. Simultaneous multichannel multi‐offset ground‐penetrating radar measurements for soil characterization
CN103809204B (en) A kind of collecting method of field audio magnetotelluric method
Küperkoch et al. Automated determination of S-phase arrival times using autoregressive prediction: application to local and regional distances
EP2135114A2 (en) Csem surveying method using an arrangement of plural signal sources
US8706462B2 (en) System and method for providing a physical property model
CN103210323B (en) Process the system and method for geological data
CN110531422A (en) A kind of tensor artificial source acquired electromagnetic data acquiring and processing method and device
Beggan et al. Geolectric field measurement, modelling and validation during geomagnetic storms in the UK
CN111708094A (en) A method and system for identifying sandstone-type uranium ore sand bodies based on wide-area electromagnetic method
CN101825722B (en) A Robust Method for Estimating the Instantaneous Frequency of Seismic Signals
Pipan et al. Multifold ground-penetrating radar and resistivity to study the stratigraphy of shallow unconsolidated sediments
CN109541703A (en) A kind of quantitative preferred method of magnetotelluric sounding curve
Li et al. Eliminating above-surface diffractions from ground-penetrating radar data using iterative Stolt migration
Sungkono et al. Application of noise-assisted multivariate empirical mode decomposition in VLF-EM data to identify underground river
Tronicke The influence of high frequency uncorrelated noise on first-break arrival times and crosshole traveltime tomography
WO2019209114A1 (en) Providing for uncertainty in non-linear inversions of geophysical data
GB2583906A (en) Method of estimating a mineral content of a geological structure
CN110109184B (en) A passive field source-like three-dimensional electric field exploration method based on multiple diurnal variation points
CN111538080B (en) Method of seismic imaging
Kang et al. Deep-learning-based airborne transient electromagnetic inversion providing the depth of investigation
Wang et al. Noise suppression of airborne transient electromagnetic data with minimum curvature

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

Application publication date: 20190329

WD01 Invention patent application deemed withdrawn after publication