CN109188542B - 一种波区相关性检测的远参考大地电磁阻抗计算方法 - Google Patents

一种波区相关性检测的远参考大地电磁阻抗计算方法 Download PDF

Info

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
Application number
CN201811338060.9A
Other languages
English (en)
Other versions
CN109188542A (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.)
National Science Chongqing Instrument Co ltd
Original Assignee
National Science Chongqing Instrument Co ltd
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 National Science Chongqing Instrument Co ltd filed Critical National Science Chongqing Instrument Co ltd
Priority to CN201811338060.9A priority Critical patent/CN109188542B/zh
Publication of CN109188542A publication Critical patent/CN109188542A/zh
Application granted granted Critical
Publication of CN109188542B publication Critical patent/CN109188542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

本发明公开了一种波区相关性检测的远参考大地电磁阻抗计算方法，包括以下步骤：S1：分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号并通过傅里叶变换转换为频率域数据并进行滤波和放大转换，分别得到基站、第一参考站、第二参考站的电场、磁场数字信号，S2‑4：根据基站的电场、磁场数字信号形成基站的电场和磁场的互功率谱，并根据互功率谱计算基站的电磁场相干度，同理得到第一、第二参考站的电磁场相干度，S5：分别计算基站、第一参考站、第二参考站彼此间的磁场分量相干度
Figure DDA0001861746620000011
Figure DDA0001861746620000012
S6：判断
Figure DDA0001861746620000013
与
Figure DDA0001861746620000014
值的大小。该方法能够明确对基站与参考站噪声的非相关性，从而根据参考站数据测到测深曲线。

Description

一种波区相关性检测的远参考大地电磁阻抗计算方法
技术领域
本发明涉及地球物理勘探领域，特别是涉及一种波区相关性检测的远参考大地电磁阻抗计算方法。
背景技术
大地电磁测探法(MT)是在地面上一点或多点同时观测天然变化的、互为垂直的电磁场分量，用以探测地球内部的电性构造，广泛用于矿场勘查、油气勘探等方面。然而在测量大地电磁数据时，常会受到噪声的干扰，严重影响后续张量阻抗估算的准确性。远参考大地电磁方法是消除噪声中非相干噪声的有效方法，但是这种处理方式利用不受非相关噪声干扰的参考站数据，并基于基站与参考站噪声的非相关性以修正本地测站的功率谱，从而得到休整之后的测深曲线，以得到张量阻抗的真实值。201610204315.7的《一种基于相关性检测的远参考大地电磁阻抗计算方法》(以下简称“前案”)中有提到能够提高张量阻抗估算准确性的方法，但是这种方法依然不能很明确基站与参考站之间的非波区信号干扰。大地电磁阻抗的估算，需要利用波区的电磁信号。常规远参考判断不了相关度高的信号是否是波区信号。
发明内容
针对上述现有技术的不足，本发明提供了一种波区相关性检测的远参考大地电磁阻抗计算方法，该方法能够明确对基站与参考站噪声的非相关性，从而根据参考站数据测到测深曲线。
为了解决上述技术问题，本发明采用了如下的技术方案：
一种波区相关性检测的远参考大地电磁阻抗计算方法，包括以下步骤：
S1：分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号，再对所述时间序列信号通过傅里叶变换转换为频率域数据，并对电场和磁场的频率域数据进行滤波和放大转换，分别得到基站、第一参考站、第二参考站的电场数字信号和磁场数字信号，在采集时需保证基站、第一参考站、第二参考站的时间段和采样率的一致性，其中，第二参考站位于第一参考站和基站之间；
S2：根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱，并根据互功率谱计算基站的电磁场相干度
Figure BDA0001861746600000021
当
Figure BDA0001861746600000022
Figure BDA0001861746600000023
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到基站电磁场不受非相关噪声干扰的数据段Segmentsst，其中，t为数据段编号，CohSSthr为预设基站电磁场相干度的阈值，其取值范围为0≤CohSSthr＜1，
Figure BDA0001861746600000024
定义如下：
Figure BDA0001861746600000025
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000026
为Ei与Hj的互功率谱，
Figure BDA0001861746600000027
为Ei的自功率谱，
Figure BDA0001861746600000028
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA0001861746600000029
脚标中的S表示基站；
S3：根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱，并根据互功率谱计算第一参考站的电磁场相干度
Figure BDA00018617466000000210
当
Figure BDA00018617466000000211
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m，其中，m为数据段编号，CohRR1thr为预设第一参考值电磁场相干度的阈值，其取值范围为0≤CohRR1thr＜1，
Figure BDA00018617466000000212
定义如下：
Figure BDA00018617466000000213
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000031
为Ei与Hj的互功率谱，
Figure BDA0001861746600000032
为Ei的自功率谱，
Figure BDA0001861746600000033
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA0001861746600000034
脚标中的R1表示第一参考站；
S4：根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱，并根据互功率谱计算第二参考站的电磁场相干度
Figure BDA0001861746600000035
当
Figure BDA0001861746600000036
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h，其中，h为数据段编号，CohRR2thr为预设第二参考值电磁场相干度的阈值，其取值范围为0≤CohRR2thr＜1，
Figure BDA0001861746600000037
定义如下：
Figure BDA0001861746600000038
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000039
为Ei与Hj的互功率谱，
Figure BDA00018617466000000310
为Ei的自功率谱，
Figure BDA00018617466000000311
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA00018617466000000312
脚标中的R2表示第二参考站；
S5：分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度
Figure BDA00018617466000000313
基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure BDA00018617466000000314
第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure BDA00018617466000000315
当
Figure BDA00018617466000000316
Figure BDA00018617466000000317
且
Figure BDA00018617466000000318
且
Figure BDA00018617466000000319
时，直接进入步骤S6，否则将
Figure BDA00018617466000000320
对应的数据段置零，其中，CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值，其取值范围为0≤CohHRSthr＜1，其中，
Figure BDA0001861746600000041
的定义如下：
Figure BDA0001861746600000042
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000043
为HSi与HR1i的互功率谱，也就是基站与第一参考站之间的磁场分量的相干度，
Figure BDA0001861746600000044
为HSi的自功率谱，
Figure BDA0001861746600000045
为HR1i的自功率谱，HSi为基站i方向的磁场强度，HR1i为第一参考站i方向的磁场强度，
Figure BDA0001861746600000046
脚标中的S表示基站，
Figure BDA0001861746600000047
脚标中的R1表示第一参考站；
Figure BDA0001861746600000048
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000049
为HSi与HR2i的互功率谱，也就是基站与第二参考站之间的磁场分量的相干度，
Figure BDA00018617466000000410
为HSi的自功率谱，
Figure BDA00018617466000000411
为HR2i的自功率谱，HSi为基站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure BDA00018617466000000412
脚标中的S表示基站，
Figure BDA00018617466000000413
脚标中的R2表示第二参考站；
Figure BDA00018617466000000414
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA00018617466000000415
为HR1i与HR2i的互功率谱，也就是第一参考站与第二参考站之间的磁场分量的相干度，
Figure BDA00018617466000000416
为HR1i的自功率谱，
Figure BDA00018617466000000417
为HR2i的自功率谱，HR1i为第一参考站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure BDA00018617466000000418
脚标中的R1表示第一参考，
Figure BDA00018617466000000419
脚标中的R2表示第二参考站；
S6：判断
Figure BDA00018617466000000420
与
Figure BDA00018617466000000421
值的大小，保留满足
Figure BDA00018617466000000422
对应的磁场和电场分量的数据段并进入S7，舍弃
Figure BDA00018617466000000423
对应的磁场和电场分量的数据段，将这些数据段置零，得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段，并统计这些数据段的段数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：计算基站的电磁场相干度
Figure BDA0001861746600000061
具体是根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱，并根据互功率谱计算基站的电磁场相干度
Figure BDA0001861746600000062
当
Figure BDA0001861746600000063
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对电场和磁场的数据段进行标记，得到基站电磁场不受非相关噪声干扰的数据段Segmentsst，其中，t为数据段编号，CohSSthr为预设基站电磁场相干度的阈值，其取值范围为0≤CohSSthr＜1，
Figure BDA0001861746600000064
定义如下：
Figure BDA0001861746600000065
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000066
为Ei与Hj的互功率谱，
Figure BDA0001861746600000067
为Ei的自功率谱，
Figure BDA0001861746600000068
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA0001861746600000069
脚标中的S表示基站；
S3：计算第一参考站的电磁场相干度
Figure BDA00018617466000000610
具体是根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱，并根据互功率谱计算第一参考站的电磁场相干度
Figure BDA00018617466000000611
当
Figure BDA00018617466000000612
Figure BDA00018617466000000613
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对电场和磁场的数据段进行标记，得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m，其中，m为数据段编号，CohRR1thr为预设第一参考值电磁场相干度的阈值，其取值范围为0≤CohRR1thr＜1，
Figure BDA0001861746600000071
定义如下：
Figure BDA0001861746600000072
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000073
为Ei与Hj的互功率谱，
Figure BDA0001861746600000074
为Ei的自功率谱，
Figure BDA0001861746600000075
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA0001861746600000076
脚标中的R1表示第一参考站；
S4：计算第二参考站的电磁场相干度
Figure BDA0001861746600000077
具体是根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱，并根据互功率谱计算第二参考站的电磁场相干度
Figure BDA0001861746600000078
当
Figure BDA0001861746600000079
Figure BDA00018617466000000710
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对电场和磁场的数据段进行标记，得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h，其中，h为数据段编号，CohRR2thr为预设第二参考值电磁场相干度的阈值，其取值范围为0≤CohRR2thr＜1，
Figure BDA00018617466000000711
定义如下：
Figure BDA00018617466000000712
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA00018617466000000713
为Ei与Hj的互功率谱，
Figure BDA00018617466000000714
为Ei的自功率谱，
Figure BDA00018617466000000715
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure BDA00018617466000000716
脚标中的R2表示第二参考站；
S5：分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度
Figure BDA00018617466000000717
基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure BDA0001861746600000081
第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure BDA0001861746600000082
当
Figure BDA0001861746600000083
Figure BDA0001861746600000084
且
Figure BDA0001861746600000085
且
Figure BDA0001861746600000086
时，直接进入步骤S6，否则将
Figure BDA0001861746600000087
对应的数据段置零，其中，CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值，其取值范围为0≤CohHRSthr＜1。只有当数据段同时满足
Figure BDA0001861746600000088
且
Figure BDA0001861746600000089
Figure BDA00018617466000000825
且
Figure BDA00018617466000000810
时，才会采用此数据段，否则就将此数据段置零。其中，
Figure BDA00018617466000000811
的定义如下：
Figure BDA00018617466000000812
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA00018617466000000813
为HSi与HR1i的互功率谱，也就是基站与第一参考站之间的磁场分量的相干度，
Figure BDA00018617466000000814
为HSi的自功率谱，
Figure BDA00018617466000000815
为HR1i的自功率谱，HSi为基站i方向的磁场强度，HR1i为第一参考站i方向的磁场强度，
Figure BDA00018617466000000816
脚标中的S表示基站，
Figure BDA00018617466000000817
脚标中的R1表示第一参考站；
Figure BDA00018617466000000818
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA00018617466000000819
为HSi与HR2i的互功率谱，也就是基站与第二参考站之间的磁场分量的相干度，
Figure BDA00018617466000000820
为HSi的自功率谱，
Figure BDA00018617466000000821
为HR2i的自功率谱，HSi为基站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure BDA00018617466000000822
脚标中的S表示基站，
Figure BDA00018617466000000823
脚标中的R2表示第二参考站；
Figure BDA00018617466000000824
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure BDA0001861746600000091
为HR1i与HR2i的互功率谱，也就是第一参考站与第二参考站之间的磁场分量的相干度，
Figure BDA0001861746600000092
为HR1i的自功率谱，
Figure BDA0001861746600000093
为HR2i的自功率谱，HR1i为第一参考站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure BDA0001861746600000094
脚标中的R1表示第一参考，
Figure BDA0001861746600000095
脚标中的R2表示第二参考站；
S6、判断
Figure BDA0001861746600000096
与
Figure BDA0001861746600000097
值的大小，当
Figure BDA0001861746600000098
时，说明采集的数据段在第一参考站和第二参考站之间，也就是说，采集点远离基站，保留满足
Figure BDA0001861746600000099
对应的磁场和电场分量的数据段并进入S7。当
Figure BDA00018617466000000910
说明采集的数据段在基站和第二参考站之间，也就是说，采集点远离第一参考点，要舍弃此数据段，将这些数据段置零，最后得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段，并统计这些数据段的段数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)
本发明的工作原理：本发明的方法是在大地电磁测探法中，在进行张量阻抗估算之前，利用基站中的磁场分量与第一参考站中磁场分量相干度
Figure BDA0001861746600000101
基站中的磁场分量与第二参考站中的磁场分量相干度
Figure BDA0001861746600000102
第一参考站中的磁场分量与第二参考站中的磁场分量相干度
Figure BDA0001861746600000103
来筛选位于第一参考站和第二参考站的数据段，避免基站与第二参考站非相关噪声的干扰，同时避免磁场噪声对张量阻抗计算结果的影响，从而得到能反映地下真实电性结构的探测曲线，提高了大地电磁的勘探效果。
最后应说明的是：本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等统计数的范围之内，则本发明也意图包含这些改动和变型。

Claims (3)

1.一种波区相关性检测的远参考大地电磁阻抗计算方法，其特征在于，包括以下步骤：
S1：分别采集并测量基站、第一参考站、第二参考站的电场和磁场分量的时间序列信号，再对所述时间序列信号通过傅里叶变换转换为频率域数据，并对电场和磁场的频率域数据进行滤波和放大转换，分别得到基站、第一参考站、第二参考站的电场数字信号和磁场数字信号，在采集时需保证基站、第一参考站、第二参考站的时间段和采样率的一致性，其中，第二参考站位于第一参考站和基站之间；
S2：根据基站的电场数字信号和磁场数字信号形成基站的电场和磁场的互功率谱，并根据互功率谱计算基站的电磁场相干度
Figure FDA0002360551360000011
当
Figure FDA0002360551360000012
Figure FDA0002360551360000013
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到基站电磁场不受非相关噪声干扰的数据段Segmentsst，其中，Segmentsst的脚标中的t为数据段编号，Segmentsst的脚标中的S为基站；CohSSthr为预设基站电磁场相干度的阈值，其取值范围为0≤CohSSthr＜1，
Figure FDA0002360551360000014
定义如下：
Figure FDA0002360551360000015
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA0002360551360000016
为Ei与Hj的互功率谱，
Figure FDA0002360551360000017
为Ei的自功率谱，
Figure FDA0002360551360000018
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure FDA0002360551360000019
脚标中的S表示基站；
S3：根据第一参考站的电场数字信号和磁场数字信号形成第一参考站的电场和磁场的互功率谱，并根据互功率谱计算第一参考站的电磁场相干度
Figure FDA00023605513600000110
当
Figure FDA00023605513600000111
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到第一参考站的电磁场不受非相关噪声干扰的数据段SegmentRR1m，其中，m为数据段编号，CohRR1thr为预设第一参考值电磁场相干度的阈值，其取值范围为0≤CohRR1thr＜1，
Figure FDA0002360551360000021
定义如下：
Figure FDA0002360551360000022
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA0002360551360000023
为Ei与Hj的互功率谱，
Figure FDA0002360551360000024
为Ei的自功率谱，
Figure FDA0002360551360000025
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure FDA0002360551360000026
脚标中的R1表示第一参考站；
S4：根据第二参考站的电场数字信号和磁场数字信号形成第二参考站的电场和磁场的互功率谱，并根据互功率谱计算第二参考站的电磁场相干度
Figure FDA0002360551360000027
当
Figure FDA0002360551360000028
时，保留互功率谱中电磁场相干度对应的电场和磁场的数据段，并对所述电场和磁场的数据段进行标记，得到第二参考站的电磁场不受非相关噪声干扰的数据段SegmentRR2h，其中，SegmentRR2h的脚标中的h为数据段编号，SegmentRR2h的脚标中的R为参考站；CohRR2thr为预设第二参考值电磁场相干度的阈值，其取值范围为0≤CohRR2thr＜1，
Figure FDA0002360551360000029
定义如下：
Figure FDA00023605513600000210
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA00023605513600000211
为Ei与Hj的互功率谱，
Figure FDA00023605513600000212
为Ei的自功率谱，
Figure FDA00023605513600000213
为Hj的自功率谱，Ei为i方向的电场强度，Hj为j方向的磁场强度，
Figure FDA00023605513600000214
脚标中的R2表示第二参考站；
S5：分别计算基站Segmentsst中的磁场分量与第一参考站SegmentRR1m中磁场分量相干度
Figure FDA0002360551360000031
基站Segmentsst中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure FDA0002360551360000032
第一参考站SegmentRR1m中的磁场分量与第二参考站SegmentRR2h中的磁场分量相干度
Figure FDA0002360551360000033
当
Figure FDA0002360551360000034
Figure FDA00023605513600000323
且
Figure FDA0002360551360000035
且
Figure FDA0002360551360000036
时，直接进入步骤S6，否则将
Figure FDA0002360551360000037
对应的数据段置零，其中，CohHRSthr为预设基站与第一参考站、第二参考站的磁场相干度的共同阈值，其取值范围为0≤CohHRSthr<1，其中，
Figure FDA0002360551360000038
的定义如下：
Figure FDA0002360551360000039
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA00023605513600000310
为HSi与HR1i的互功率谱，也就是基站与第一参考站之间的磁场分量的相干度，
Figure FDA00023605513600000311
为HSi的自功率谱，
Figure FDA00023605513600000312
为HR1i的自功率谱，HSi为基站i方向的磁场强度，HR1i为第一参考站i方向的磁场强度，
Figure FDA00023605513600000313
脚标中的S表示基站，
Figure FDA00023605513600000314
脚标中的R1表示第一参考站；
Figure FDA00023605513600000315
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA00023605513600000316
为HSi与HR2i的互功率谱，也就是基站与第二参考站之间的磁场分量的相干度，
Figure FDA00023605513600000317
为HSi的自功率谱，
Figure FDA00023605513600000318
为HR2i的自功率谱，HSi为基站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure FDA00023605513600000319
脚标中的S表示基站，
Figure FDA00023605513600000320
脚标中的R2表示第二参考站；
Figure FDA00023605513600000321
其中，i、j表示x、y方向，n为独立观测的数据段段数，k为独立观测的数据段的序号，
Figure FDA00023605513600000322
为HR1i与HR2i的互功率谱，也就是第一参考站与第二参考站之间的磁场分量的相干度，
Figure FDA0002360551360000041
为HR1i的自功率谱，
Figure FDA0002360551360000042
为HR2i的自功率谱，HR1i为第一参考站i方向的磁场强度，HR2i为第二参考站i方向的磁场强度，
Figure FDA0002360551360000043
脚标中的R1表示第一参考，
Figure FDA0002360551360000044
脚标中的R2表示第二参考站；
S6、判断
Figure FDA0002360551360000045
与
Figure FDA0002360551360000046
值的大小，保留满足
Figure FDA0002360551360000047
对应的磁场和电场分量的数据段并进入S7，舍弃
Figure FDA0002360551360000048
对应的磁场和电场分量的数据段，将这些数据段置零，得到磁场不受噪声干扰的第一参考站与第二参考站之间的数据段，并统计这些数据段的段数SegmentHSR；
S7、将数据段SegmentHSR运用基于Robust的远参考大地电磁阻抗计算方法进行张量阻抗估算得到张量估算值，并利用张量估算值计算电阻率和阻抗相位值。
2.根据权利要求1所述的一种波区相关性检测的远参考大地电磁阻抗计算方法，其特征在于，对所述时间序列信号做傅里叶变换转换的具体方法为：
1.1)对每个单位时间段做傅里叶变换；
1.2)取每个单位时间段的中点为新端点；
1.3)将新端点形成的新的单位时间段做傅里叶变换，并与原端点做的傅里叶变换进行叠加。
3.根据权利要求1所述的一种波区相关性检测的远参考大地电磁阻抗计算方法，其特征在于，所述CohSSthr、CohRR1thr、CohRR2thr、CohHRSthr的区间为[0.6，0.9)。
CN201811338060.9A 2018-11-12 2018-11-12 一种波区相关性检测的远参考大地电磁阻抗计算方法 Active CN109188542B (zh)

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)

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

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

* Cited by examiner, † Cited by third party
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 西南科技大学 一种基于相关性检测的远参考大地电磁阻抗计算方法

Patent Citations (5)

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

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