CN109188542B - 一种波区相关性检测的远参考大地电磁阻抗计算方法 - Google Patents
一种波区相关性检测的远参考大地电磁阻抗计算方法 Download PDFInfo
- Publication number
- CN109188542B CN109188542B CN201811338060.9A CN201811338060A CN109188542B CN 109188542 B CN109188542 B CN 109188542B CN 201811338060 A CN201811338060 A CN 201811338060A CN 109188542 B CN109188542 B CN 109188542B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- reference station
- power spectrum
- base station
- coherence
- 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.)
- Active
Links
Images
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)
- Mobile Radio Communication Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
技术领域
本发明涉及地球物理勘探领域,特别是涉及一种波区相关性检测的远参考大地电磁阻抗计算方法。
背景技术
大地电磁测探法(MT)是在地面上一点或多点同时观测天然变化的、互为垂直的电磁场分量,用以探测地球内部的电性构造,广泛用于矿场勘查、油气勘探等方面。然而在测量大地电磁数据时,常会受到噪声的干扰,严重影响后续张量阻抗估算的准确性。远参考大地电磁方法是消除噪声中非相干噪声的有效方法,但是这种处理方式利用不受非相关噪声干扰的参考站数据,并基于基站与参考站噪声的非相关性以修正本地测站的功率谱,从而得到休整之后的测深曲线,以得到张量阻抗的真实值。201610204315.7的《一种基于相关性检测的远参考大地电磁阻抗计算方法》(以下简称“前案”)中有提到能够提高张量阻抗估算准确性的方法,但是这种方法依然不能很明确基站与参考站之间的非波区信号干扰。大地电磁阻抗的估算,需要利用波区的电磁信号。常规远参考判断不了相关度高的信号是否是波区信号。
发明内容
针对上述现有技术的不足,本发明提供了一种波区相关性检测的远参考大地电磁阻抗计算方法,该方法能够明确对基站与参考站噪声的非相关性,从而根据参考站数据测到测深曲线。
为了解决上述技术问题,本发明采用了如下的技术方案:
一种波区相关性检测的远参考大地电磁阻抗计算方法,包括以下步骤:
S1:分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号,再对所述时间序列信号通过傅里叶变换转换为频率域数据,并对电场和磁场的频率域数据进行滤波和放大转换,分别得到基站、第一参考站、第二参考站的电场数字信号和磁场数字信号,在采集时需保证基站、第一参考站、第二参考站的时间段和采样率的一致性,其中,第二参考站位于第一参考站和基站之间;
S2:根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱,并根据互功率谱计算基站的电磁场相干度当 时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到基站电磁场不受非相关噪声干扰的数据段Segmentsst,其中,t为数据段编号,CohSSthr为预设基站电磁场相干度的阈值,其取值范围为0≤CohSSthr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的S表示基站;
S3:根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱,并根据互功率谱计算第一参考站的电磁场相干度当时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m,其中,m为数据段编号,CohRR1thr为预设第一参考值电磁场相干度的阈值,其取值范围为0≤CohRR1thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R1表示第一参考站;
S4:根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱,并根据互功率谱计算第二参考站的电磁场相干度当时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h,其中,h为数据段编号,CohRR2thr为预设第二参考值电磁场相干度的阈值,其取值范围为0≤CohRR2thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R2表示第二参考站;
S5:分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度当 且且时,直接进入步骤S6,否则将对应的数据段置零,其中,CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值,其取值范围为0≤CohHRSthr<1,其中,的定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR1i的互功率谱,也就是基站与第一参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR1i的自功率谱,HSi为基站i方向的磁场强度,HR1i为第一参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R1表示第一参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR2i的互功率谱,也就是基站与第二参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR2i的自功率谱,HSi为基站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R2表示第二参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HR1i与HR2i的互功率谱,也就是第一参考站与第二参考站之间的磁场分量的相干度,为HR1i的自功率谱,为HR2i的自功率谱,HR1i为第一参考站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的R1表示第一参考,脚标中的R2表示第二参考站;
S6:判断与值的大小,保留满足对应的磁场和电场分量的数据段并进入S7,舍弃对应的磁场和电场分量的数据段,将这些数据段置零,得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段,并统计这些数据段的段数SegmentHSR;
S7:将数据段SegmentHSR运用基于Robust的远参考大地电磁阻抗计算方法进行张量阻抗估算得到张量估算值,并利用张量估算值计算电阻率和阻抗相位值。
作为优化,对所述时间序列信号做傅里叶变换转换的具体方法为:
1.1)对每个单位时间段做傅里叶变换;
1.2)取每个单位时间段的中点为新端点;
1.3)将新端点形成的新的单位时间段做傅里叶变换,并与原端点做的傅里叶变换进行叠加。
作为优化,所述CohSSthr、CohRR1thr、CohRR2thr、CohHRSthr的优选区间为[0.6,0.9)。
本发明的有益效果是:
1.通过对多层时间序列信号做傅里叶叠加,提高了采集数据的准确性,同时。
2.在第一参考站和基站之间设立第二参考站,提升了采集数据基于基站与第一参考站噪声的非相关性,从而保证了第一参考站所受的干扰与基站噪声无关。
附图说明
图1为本发明所述的一种基于波区相关性检测的远参考大地电磁阻抗计算方法的流程示意图;
图2为本发明所述的一种基于波区相关性检测的远参考大地电磁阻抗计算方法的的傅里叶变换的示意图。
具体实施方式
下面结合附图对本发明作进一步的详细说明。
如图1所示,一种基于波区相关性检测的远参考大地电磁阻抗计算方法,包括以下步骤:
S1:分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号,再对时间序列信号通过傅里叶变换转换为频率域数据,并对电场和磁场的频率域数据进行滤波和放大转换,分别得到基站、第一参考站、第二参考站的电场数字信号和磁场数字信号,在采集时需保证基站、第一参考站、第二参考站的时间段和采样率的一致性,其中,第二参考站位于第一参考站和基站之间,本实施例中,采样率fs=1。
S2:计算基站的电磁场相干度具体是根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱,并根据互功率谱计算基站的电磁场相干度当时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对电场和磁场的数据段进行标记,得到基站电磁场不受非相关噪声干扰的数据段Segmentsst,其中,t为数据段编号,CohSSthr为预设基站电磁场相干度的阈值,其取值范围为0≤CohSSthr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的S表示基站;
S3:计算第一参考站的电磁场相干度具体是根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱,并根据互功率谱计算第一参考站的电磁场相干度当 时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对电场和磁场的数据段进行标记,得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m,其中,m为数据段编号,CohRR1thr为预设第一参考值电磁场相干度的阈值,其取值范围为0≤CohRR1thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R1表示第一参考站;
S4:计算第二参考站的电磁场相干度具体是根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱,并根据互功率谱计算第二参考站的电磁场相干度当 时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对电场和磁场的数据段进行标记,得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h,其中,h为数据段编号,CohRR2thr为预设第二参考值电磁场相干度的阈值,其取值范围为0≤CohRR2thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R2表示第二参考站;
S5:分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度当 且且时,直接进入步骤S6,否则将对应的数据段置零,其中,CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值,其取值范围为0≤CohHRSthr<1。只有当数据段同时满足且 且时,才会采用此数据段,否则就将此数据段置零。其中,的定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR1i的互功率谱,也就是基站与第一参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR1i的自功率谱,HSi为基站i方向的磁场强度,HR1i为第一参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R1表示第一参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR2i的互功率谱,也就是基站与第二参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR2i的自功率谱,HSi为基站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R2表示第二参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HR1i与HR2i的互功率谱,也就是第一参考站与第二参考站之间的磁场分量的相干度,为HR1i的自功率谱,为HR2i的自功率谱,HR1i为第一参考站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的R1表示第一参考,脚标中的R2表示第二参考站;
S6、判断与值的大小,当时,说明采集的数据段在第一参考站和第二参考站之间,也就是说,采集点远离基站,保留满足对应的磁场和电场分量的数据段并进入S7。当说明采集的数据段在基站和第二参考站之间,也就是说,采集点远离第一参考点,要舍弃此数据段,将这些数据段置零,最后得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段,并统计这些数据段的段数SegmentHSR。
S7、将数据段SegmentHSR运用基于Robust的远参考大地电磁阻抗计算方法进行张量阻抗估算得到张量估算值,并利用张量估算值计算电阻率和阻抗相位值。
本实施例中,对时间序列信号做傅里叶变换转换的具体方法为:
1.1)对每个单位时间段做傅里叶变换。本实施例中,单位时间段为1秒。如图2所示,a1a2、a2a3、a3a4、a4a5、a5a6分别为一个单位时间段,分别对此单位时间段做傅里叶变换得到FFT1、FFT2、FFT3、FFT4、FFT5。
1.2)取每个单位时间段的中点为新端点。如图2所示,b1、b2、b3、b4、b5为新端点。
1.3)将新端点形成的新的单位时间段做傅里叶变换,并与原端点做的傅里叶变换进行叠加。新的单位时间段为b1b2、b2b3、b3b4、b4b5,对应的傅里叶变换为FFT6、FFT7、FFT8、FFT9,最终的傅里叶变换为FFT=FFT1+FFT2+FFT3+FFT4+FFT5+FFT6+FFT7+FFT8+FFT9。这样可以对更多时间序列信号做傅里叶变换,使得到的数据更加精确。
本实施例中,CohSSthr、CohRR1thr、CohRR2thr、CohHRSthr的优选区间为[0.6,0.9)
本发明的工作原理:本发明的方法是在大地电磁测探法中,在进行张量阻抗估算之前,利用基站中的磁场分量与第一参考站中磁场分量相干度基站中的磁场分量与第二参考站中的磁场分量相干度第一参考站中的磁场分量与第二参考站中的磁场分量相干度来筛选位于第一参考站和第二参考站的数据段,避免基站与第二参考站非相关噪声的干扰,同时避免磁场噪声对张量阻抗计算结果的影响,从而得到能反映地下真实电性结构的探测曲线,提高了大地电磁的勘探效果。
最后应说明的是:本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等统计数的范围之内,则本发明也意图包含这些改动和变型。
Claims (3)
1.一种波区相关性检测的远参考大地电磁阻抗计算方法,其特征在于,包括以下步骤:
S1:分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号,再对所述时间序列信号通过傅里叶变换转换为频率域数据,并对电场和磁场的频率域数据进行滤波和放大转换,分别得到基站、第一参考站、第二参考站的电场数字信号和磁场数字信号,在采集时需保证基站、第一参考站、第二参考站的时间段和采样率的一致性,其中,第二参考站位于第一参考站和基站之间;
S2:根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱,并根据互功率谱计算基站的电磁场相干度当 时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到基站电磁场不受非相关噪声干扰的数据段Segmentsst,其中,Segmentsst的脚标中的t为数据段编号,Segmentsst的脚标中的S为基站;CohSSthr为预设基站电磁场相干度的阈值,其取值范围为0≤CohSSthr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的S表示基站;
S3:根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱,并根据互功率谱计算第一参考站的电磁场相干度当时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m,其中,m为数据段编号,CohRR1thr为预设第一参考值电磁场相干度的阈值,其取值范围为0≤CohRR1thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R1表示第一参考站;
S4:根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱,并根据互功率谱计算第二参考站的电磁场相干度当时,保留互功率谱中电磁场相干度对应的电场和磁场的数据段,并对所述电场和磁场的数据段进行标记,得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h,其中,SegmentRR2h的脚标中的h为数据段编号,SegmentRR2h的脚标中的R为参考站;CohRR2thr为预设第二参考值电磁场相干度的阈值,其取值范围为0≤CohRR2thr<1,定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为Ei与Hj的互功率谱,为Ei的自功率谱,为Hj的自功率谱,Ei为i方向的电场强度,Hj为j方向的磁场强度,脚标中的R2表示第二参考站;
S5:分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度当 且且时,直接进入步骤S6,否则将对应的数据段置零,其中,CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值,其取值范围为0≤CohHRSthr<1,其中,的定义如下:
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR1i的互功率谱,也就是基站与第一参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR1i的自功率谱,HSi为基站i方向的磁场强度,HR1i为第一参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R1表示第一参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HSi与HR2i的互功率谱,也就是基站与第二参考站之间的磁场分量的相干度,为HSi的自功率谱,为HR2i的自功率谱,HSi为基站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的S表示基站,脚标中的R2表示第二参考站;
其中,i、j表示x、y方向,n为独立观测的数据段段数,k为独立观测的数据段的序号,为HR1i与HR2i的互功率谱,也就是第一参考站与第二参考站之间的磁场分量的相干度,为HR1i的自功率谱,为HR2i的自功率谱,HR1i为第一参考站i方向的磁场强度,HR2i为第二参考站i方向的磁场强度,脚标中的R1表示第一参考,脚标中的R2表示第二参考站;
S6、判断与值的大小,保留满足对应的磁场和电场分量的数据段并进入S7,舍弃对应的磁场和电场分量的数据段,将这些数据段置零,得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段,并统计这些数据段的段数SegmentHSR;
S7、将数据段SegmentHSR运用基于Robust的远参考大地电磁阻抗计算方法进行张量阻抗估算得到张量估算值,并利用张量估算值计算电阻率和阻抗相位值。
2.根据权利要求1所述的一种波区相关性检测的远参考大地电磁阻抗计算方法,其特征在于,对所述时间序列信号做傅里叶变换转换的具体方法为:
1.1)对每个单位时间段做傅里叶变换;
1.2)取每个单位时间段的中点为新端点;
1.3)将新端点形成的新的单位时间段做傅里叶变换,并与原端点做的傅里叶变换进行叠加。
3.根据权利要求1所述的一种波区相关性检测的远参考大地电磁阻抗计算方法,其特征在于,所述CohSSthr、CohRR1thr、CohRR2thr、CohHRSthr的区间为[0.6,0.9)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811338060.9A CN109188542B (zh) | 2018-11-12 | 2018-11-12 | 一种波区相关性检测的远参考大地电磁阻抗计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811338060.9A CN109188542B (zh) | 2018-11-12 | 2018-11-12 | 一种波区相关性检测的远参考大地电磁阻抗计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109188542A CN109188542A (zh) | 2019-01-11 |
CN109188542B true CN109188542B (zh) | 2020-04-14 |
Family
ID=64938763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811338060.9A Active CN109188542B (zh) | 2018-11-12 | 2018-11-12 | 一种波区相关性检测的远参考大地电磁阻抗计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109188542B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110058319B (zh) * | 2019-01-16 | 2020-11-27 | 南方科技大学 | 一种大地电磁数据采集方法、装置及终端设备 |
CN111190234B (zh) * | 2020-01-13 | 2021-03-19 | 山东大学 | 一种针对人工电性源频率域电磁法的噪声观测方法与装置 |
CN111736226B (zh) * | 2020-06-30 | 2021-06-11 | 吉林大学 | 一种确定大地电磁时间域阻抗的方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686476A (en) * | 1981-11-23 | 1987-08-11 | Mobil Oil Corporation | Plotting of magnetotelluric impedance to determine structure and resistivity variation of anomalies |
CN101142498A (zh) * | 2005-01-19 | 2008-03-12 | Ksn能源有限责任公司 | 利用电磁阻抗层析术的用于石油开采中的温度测量和流体流量的地下成像术 |
CN102116872A (zh) * | 2009-12-31 | 2011-07-06 | 核工业北京地质研究院 | 一种stratagem大地电磁测量系统阻抗张量的稳健估算方法 |
CN102944901A (zh) * | 2012-11-30 | 2013-02-27 | 中国船舶重工集团公司第七二二研究所 | 一种大地电磁阻抗估计方法 |
CN105629317A (zh) * | 2016-04-08 | 2016-06-01 | 中国矿业大学(北京) | 一种基于站间传递函数的大地电磁噪声压制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104656156A (zh) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | 音频大地电磁测深三维采集资料的磁参考处理方法 |
CN105204078B (zh) * | 2014-06-09 | 2016-11-23 | 中国石油化工股份有限公司 | 一种计算三维大地电磁测深主阻抗的方法 |
CN105866852B (zh) * | 2016-04-01 | 2017-10-13 | 西南科技大学 | 一种基于相关性检测的远参考大地电磁阻抗计算方法 |
-
2018
- 2018-11-12 CN CN201811338060.9A patent/CN109188542B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686476A (en) * | 1981-11-23 | 1987-08-11 | Mobil Oil Corporation | Plotting of magnetotelluric impedance to determine structure and resistivity variation of anomalies |
CN101142498A (zh) * | 2005-01-19 | 2008-03-12 | Ksn能源有限责任公司 | 利用电磁阻抗层析术的用于石油开采中的温度测量和流体流量的地下成像术 |
CN102116872A (zh) * | 2009-12-31 | 2011-07-06 | 核工业北京地质研究院 | 一种stratagem大地电磁测量系统阻抗张量的稳健估算方法 |
CN102944901A (zh) * | 2012-11-30 | 2013-02-27 | 中国船舶重工集团公司第七二二研究所 | 一种大地电磁阻抗估计方法 |
CN105629317A (zh) * | 2016-04-08 | 2016-06-01 | 中国矿业大学(北京) | 一种基于站间传递函数的大地电磁噪声压制方法 |
Non-Patent Citations (1)
Title |
---|
重复中位数估计Robust算法在长周期大地电磁中的应用;张刚等;《物探与化探》;20111231;第35卷(第6期);第813-816、828页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109188542A (zh) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109188542B (zh) | 一种波区相关性检测的远参考大地电磁阻抗计算方法 | |
Jones et al. | A comparison of techniques for magnetotelluric response function estimation | |
CN105549097A (zh) | 一种瞬变电磁信号工频及其谐波干扰消除方法及装置 | |
Li et al. | Wavelet-based higher order correlative stacking for seismic data denoising in the curvelet domain | |
CN109581516B (zh) | 曲波域统计量自适应阈值探地雷达数据去噪方法及系统 | |
CN105204073B (zh) | 一种张量视电导率测量方法 | |
CN108345039A (zh) | 一种消除地面核磁共振数据中邻频谐波干扰的方法 | |
Zhou et al. | Migration velocity analysis and prestack migration of common-transmitter GPR data | |
JI et al. | A study on solution of transient electromagnetic response during transmitting current turn‐off in the ATTEM system | |
CN104849590A (zh) | 一种混合噪声干扰下微弱脉冲信号检测方法 | |
Liu et al. | A modified empirical mode decomposition method for multiperiod time-series detrending and the application in full-waveform induced polarization data | |
CN110554428A (zh) | 一种基于变分模态分解的地震波低频能量变化率提取方法 | |
CN109117816A (zh) | 基于六阶样条插值小波的信号奇异点检测方法 | |
CN114265124A (zh) | 基于时间域瞬变电磁概率反演的不良地质体定位方法 | |
CN109581481B (zh) | 一种便携式高频可控震源地震信号谐波干扰消除方法 | |
CN112766224B (zh) | 畸变信号中提取真实信号的方法、装置、设备和存储介质 | |
Sørensen et al. | TDEM in groundwater mapping—a continuous approach | |
CN108776357A (zh) | 沉积地层瞬变电磁法电磁干扰的校正方法及装置 | |
LIN et al. | Segmented time-frequency peak filtering for random noise reduction of MRS oscillating signal | |
CN107829729B (zh) | 过套管微分电阻率测井的频率域信号处理方法 | |
CN106125132B (zh) | 含单频干扰地震道的迭代识别和压制方法 | |
Sungkono et al. | Application of Noise-Assisted Multivariate Empirical Mode Decomposition in VLF-EM Data to Identify Underground River | |
CN110109184B (zh) | 一种基于多日变点的被动场源类三维电场勘探方法 | |
CN110135022B (zh) | 一种基于极化介质模型的广义趋肤深度计算方法 | |
CN115407412A (zh) | 一种井旁异常体时间域电磁波探测方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |