CN111983701B - 横向大地电磁波测深方法及其测深装置 - Google Patents
横向大地电磁波测深方法及其测深装置 Download PDFInfo
- Publication number
- CN111983701B CN111983701B CN202010905975.4A CN202010905975A CN111983701B CN 111983701 B CN111983701 B CN 111983701B CN 202010905975 A CN202010905975 A CN 202010905975A CN 111983701 B CN111983701 B CN 111983701B
- Authority
- CN
- China
- Prior art keywords
- measuring point
- sounding
- magnetotelluric
- transverse
- power spectrum
- 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
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/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种横向大地电磁波测深方法,包括选定测线、主测点和辅助测点;在主测点及辅助测点布设传感器、获取监测时间序列并估算功率谱;计算各辅助测点与主测点之间的相位差频率曲线或空间自相关系数频率曲线并计算大地电阻率,完成横向大地电磁波测深。本发明还公开了实现所述横向大地电磁波测深方法的测深装置。本发明可以通过相速度计算电阻率、也可以结合传统大地电磁法计算卡尼亚电阻率,实现一个装置完成两种方法的测深工作,提高测量精度。最后,本发明是基于横向大地电磁波,适用于人工源或天然源,方便快捷、简单高效。
Description
技术领域
本发明属于地球物理勘探领域,具体涉及一种横向大地电磁波测深方法及其测深装置。
背景技术
随着经济技术的发展和人们生活水平的提高,人们对于基础物理领域的研究也越来越深入。
在地球物理勘探领域,1950年杰洪诺夫(Tichonov)和1953年卡尼亚(Cagnaird)提出了基于平面电磁波理论的大地电磁测深法(MT),并逐步发展出音频大地电磁测深法(AMT)和可控源音频大地电磁测深法(CSAMT),统称为平面波大地电磁测深法。这些方法都是基于平面电磁波垂直入射至大地的理论基础上。实际上不管是人工源还是天然源,大地电磁平面波垂直入射都是简化的理想状态,由于复杂地形、各向异性和三维不均匀体的影响,实际中的大地电磁波可能沿各个方向扩散传播,但都可以分解为沿垂直方向传播的波(定义为纵向大地电磁波)和沿水平方向传播的波(定义为横向大地电磁波)的矢量和。
目前,传统的大地电磁测深法在应用时,由于需要观测的数据较多,使得观测装置极为沉重,而且深部数据还容易受到浅表不均匀地电体或者地形起伏引起的静态效应影响。
发明内容
本发明的目的之一在于提供一种基于横向大地电磁波,而且方便快捷、简单高效,不受静态效应影响的横向大地电磁波测深方法。
本发明的目的之二在于提供一种实现所述横向大地电磁波测深方法的测深装置。
本发明提供的这种横向大地电磁波测深方法,包括如下步骤:
S1.选定测线,并在测线上选定主测点;
S2.在步骤S1选定的主测点周围,选定N个辅助测点;N为不小于2的正整数;所述辅助测点均匀分布在以主测点为中心、半径为L的圆周上;
S3.在步骤S2选定的主测点及各个辅助测点处沿平行于测线方向和垂直于测线方向布设传感器,同步监测并获取各测点上的监测时间序列;
S4.对步骤S3得到的各监测时间序列进行功率谱估算;
S5.根据步骤S4得到的功率谱计算各辅助测点与主测点之间的相位差频率曲线或空间自相关系数频率曲线;
S6.根据步骤S5得到的相位差频率曲线或空间自相关系数频率曲线,计算大地电阻率,完成横向大地电磁波测深。
步骤S3所述的传感器,具体为电场传感器或磁场传感器。
步骤S4所述的功率谱估算,具体为估算各监测时间序列的自功率谱以及主测点与辅助测点监测信号之间的互功率谱;
步骤S5中各辅助测点与主测点的相位差频率曲线,具体为采用如下算式计算相位差频率曲线:
式中:为第i辅助测点和主测点的相位差频率曲线;angle()为取复数辐角主值操作,/>表示第i辅助测点和主测点间监测信号互功率谱,ω为角频率。
步骤S5中各辅助测点与主测点的空间自相关系数频率曲线,具体为采用如下算式计算空间自相关系数频率曲线:
式中表示第i辅助测点和主测点间监测信号互功率谱的实部,表示第i辅助测点监测信号的自功率谱,/>表示主测点监测信号自功率谱;表示主测点与第i辅助测点连线与测线的夹角,L为辅助测点和主测点间距离。
步骤S6所述的计算大地电阻率,具体为采用如下算式计算大地电阻率:
式中:μ为磁导率,ω为角频率,k与ω相关且计算公式如下:
若N=2,则
若N>2,则针对频率ω0,将空间自相关系数频率曲线与0阶贝塞尔函数拟合,采用如下算式计算拟合误差其中J0为0阶贝塞尔函数;并取拟合误差δ(ω0)最小时所对应的k的取值为最终的k值。
本发明还公开了一种实现所述横向大地电磁波测深方法的测深装置,包括测点传感器模块和上位机模块;测点传感器模块和上位机模块连接;测点传感器模块安装在测点上,用于监测并获取各测点上的监测时间序列并上传上位机模块;上位机模块用于根据接收的监测时间序列进行功率谱统计估算和计算相位差频率曲线或空间自相关系数频率曲线,并计算得到最终的大地电阻率。
所述的测点传感器模块为电场传感器模块或磁场传感器模块。
本发明提供的这种横向大地电磁波测深方法及其测深装置,可以计算电场与电场、磁场与磁场的相速度电阻率,也可以计算磁场与电场的相速度电阻率。本发明方案观测装置多样化,可适用不同地形条件,而且可以单独采用本发明方案实现电磁测深,也可以结合传统大地电磁法,实现一个装置完成两种方法的测深工作。最后,本发明基于横向大地电磁波,适用于人工源或天然源,方便快捷、简单高效。
附图说明
图1为本发明方法的方法流程示意图。
图2为本发明装置的功能模块图。
具体实施方式
如图1所示为本发明方法的方法流程示意图:本发明提供的这种横向大地电磁波测深方法,包括如下步骤:
S1.选定测线,并在测线上选定主测点;
S2.在步骤S1选定的主测点周围,选定N个辅助测点;N为不小于2的正整数;所述辅助测点均匀分布在以主测点为中心、半径为L的圆周上;
S3.在步骤S2选定的主测点及各个辅助测点处沿平行于测线方向和垂直于测线方向布设传感器,同步监测并获取各测点上的监测时间序列;
在具体实施时,布设的传感器可以为电场传感器或磁场传感器;
S4.对步骤S3得到的各监测时间序列进行功率谱估算;具体为估算各监测时间序列的自功率谱以及主测点与辅助测点监测信号之间的互功率谱;
S5.根据步骤S4得到的功率谱计算各辅助测点与主测点之间的相位差频率曲线或空间自相关系数频率曲线;
在具体实施时,采用如下算式计算相位差频率曲线:
式中:为第i辅助测点和主测点的相位差频率曲线;angle()为取复数辐角主值操作,/>表示第i辅助测点和主测点间监测信号互功率谱,ω为角频率;
采用如下算式计算空间自相关系数频率曲线:
式中表示第i辅助测点和主测点间监测信号互功率谱的实部,表示第i辅助测点监测信号的自功率谱,/>表示主测点监测信号自功率谱;表示主测点与第i辅助测点连线与测线的夹角,L为辅助测点和主测点间距离;
S6.根据步骤S5得到的相位差频率曲线或空间自相关系数频率曲线,计算大地电阻率,完成横向大地电磁波测深;具体为采用如下算式计算大地电阻率:
式中:μ为磁导率,ω为角频率,k与ω相关且计算公式如下:
若N=2,则
若N>2,则针对频率ω0,将空间自相关系数频率曲线与0阶贝塞尔函数拟合,采用如下算式计算拟合误差其中J0为0阶贝塞尔函数;并取拟合误差δ(ω0)最小时所对应的k的取值为最终的k值。
如图2所述为本发明装置的功能模块图:本发明还公开了一种实现所述横向大地电磁波测深方法的测深装置,包括测点传感器模块和上位机模块;测点传感器模块和上位机模块连接;测点传感器模块安装在测点上,用于监测并获取各测点上的监测时间序列并上传上位机模块;上位机模块用于根据接收的监测时间序列进行功率谱统计估算和计算相位差频率曲线或空间自相关系数频率曲线,并计算得到最终的大地电阻率。
在具体实施时,测点传感器模块为电场传感器模块或磁场传感器模块。
横向大地电磁波测深方法基于相速度与大地电阻的关系实现测深,由于电磁波在相同的媒质中相速度是相同的,对于地形起伏,传播媒质不变,相速度不变,所以根据相速度计算得到的电阻率也不受地形影响,因而不会出现传统大地电磁法中由于地形起伏引起的静态效应现象。而浅表不均匀地电体仅影响其深度附近一定体积范围内的相速度,并不影响距离其较远的深度范围的相速度,因而也不会出现传统大地电磁测深中的电场振幅引起的静态效应现象。
本发明方法的理论基础为:对于纵向大地电磁波,其等相位面为水平面,也就是说当同步观测的两测点处于相同水平面时,其纵向大地电磁波相位差为0。
对于横向大地电磁波,其等相位面为垂直面,也即是说水平面上相隔一定距离的同步观测的两个测点存在相位差,该相位差与横向大地电磁波在大地传播的波速相关,而该波速与大地电导率相关,因此可以通过研究两测点间相位差来评估大地电导率。
根据电磁波传播理论,角频率ω的电磁波在导电媒质中相速度v表达式为:
式中ω=2πf为角频率,k为相位常数,单位为rad/m;ε为介电常数;μ为磁导率;σ为电导率;真空中,μ=μ0=4π×10-7H/m;
对于大地电磁波常用频段及常见大地岩矿电导率分布范围,一般有所以相位常数/>所以相速度为
此外,根据传播距离差L带来的相位差也可以计算得到相速度
联立以上公式可得电阻率
电阻率即为该频率趋肤深度内的视电阻率。
由于天然大地电磁波来自各个方向且不固定,为提高观测质量,可采取多个辅助测点辅助观测,辅助测点在主测点为圆心的外围圆周上均匀分布,通过观测信号统计的空间自相关系数频率函数来统计求得相位常数,进一步求得电阻率。
此外,以上本发明所提供的方法和装置,是针对剖面测量,实际应用中也可以用于面积观测。面积性观测时,各横向剖面和各纵向剖面分别按上述方法分析即可。
Claims (7)
1.一种横向大地电磁波测深方法,包括如下步骤:
S1.选定测线,并在测线上选定主测点;
S2.在步骤S1选定的主测点周围,选定N个辅助测点;N为不小于2的正整数;所述辅助测点均匀分布在以主测点为中心、半径为L的圆周上;
S3.在步骤S2选定的主测点及各个辅助测点处沿平行于测线方向和垂直于测线方向布设传感器,同步监测并获取各测点上的监测时间序列;
S4.对步骤S3得到的各监测时间序列进行功率谱估算;
S5.根据步骤S4得到的功率谱计算各辅助测点与主测点之间的相位差频率曲线或空间自相关系数频率曲线;
其中,所述的各辅助测点与主测点的空间自相关系数频率曲线,具体为采用如下算式计算:
式中表示第i辅助测点和主测点间监测信号互功率谱的实部,/>表示第i辅助测点监测信号的自功率谱,/>表示主测点监测信号自功率谱;/>表示主测点与第i辅助测点连线与测线的夹角,L为辅助测点和主测点间距离;
S6.根据步骤S5得到的相位差频率曲线或空间自相关系数频率曲线,计算大地电阻率,完成横向大地电磁波测深。
2.根据权利要求1所述的横向大地电磁波测深方法,其特征在于:步骤S3所述的传感器,具体为电场传感器或磁场传感器。
3.根据权利要求1或2所述的横向大地电磁波测深方法,其特征在于:步骤S4所述的功率谱估算,具体为估算各监测时间序列的自功率谱以及主测点与辅助测点监测信号之间的互功率谱。
4.根据权利要求3所述的横向大地电磁波测深方法,其特征在于:步骤S5中各辅助测点与主测点的相位差频率曲线,具体为采用如下算式计算:
式中:为第i辅助测点和主测点的相位差频率曲线;angle()为取复数辐角主值操作,/>表示第i辅助测点和主测点间监测信号互功率谱,ω为角频率。
5.根据权利要求4所述的横向大地电磁波测深方法,其特征在于:步骤S6所述的计算大地电阻率,具体为采用如下算式计算:
式中:μ为磁导率,ω为角频率,k与ω相关且计算公式如下:
若N=2,则
若N>2,则针对频率ω0,将空间自相关系数频率曲线与0阶贝塞尔函数拟合,采用如下算式计算拟合误差其中J0为0阶贝塞尔函数;并取拟合误差δ(ω0)最小时所对应的k的取值为最终的k值。
6.一种实现权利要求1~5之一所述的横向大地电磁波测深方法的测深装置,其特征在于:包括测点传感器模块和上位机模块;测点传感器模块和上位机模块连接;测点传感器模块安装在测点上,用于监测并获取各测点上的监测时间序列并上传上位机模块;上位机模块用于根据接收的监测时间序列进行功率谱统计估算和计算相位差频率曲线或空间自相关系数频率曲线,并计算得到最终的大地电阻率。
7.根据权利要求6所述的测深装置,其特征在于:所述的测点传感器模块为电场传感器模块或磁场传感器模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010905975.4A CN111983701B (zh) | 2020-09-01 | 2020-09-01 | 横向大地电磁波测深方法及其测深装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010905975.4A CN111983701B (zh) | 2020-09-01 | 2020-09-01 | 横向大地电磁波测深方法及其测深装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111983701A CN111983701A (zh) | 2020-11-24 |
CN111983701B true CN111983701B (zh) | 2023-10-10 |
Family
ID=73448369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010905975.4A Active CN111983701B (zh) | 2020-09-01 | 2020-09-01 | 横向大地电磁波测深方法及其测深装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111983701B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102419456A (zh) * | 2011-06-30 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 瞬变电磁测深数据的直接时间域处理方法 |
CN103389514A (zh) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | 一种大地电磁去噪装置及大地电磁去噪方法 |
CN104102814A (zh) * | 2014-06-11 | 2014-10-15 | 中国科学院地质与地球物理研究所 | 一种基于大地电磁数据反演电阻率和磁化率的方法及系统 |
CN105911603A (zh) * | 2016-05-04 | 2016-08-31 | 湖南科技大学 | 基于天然电场的四维物探方法 |
CN109765627A (zh) * | 2019-03-08 | 2019-05-17 | 中南大学 | 适用于山地的平面波电磁测深方法 |
CN110488365A (zh) * | 2019-09-11 | 2019-11-22 | 湖南五维地质科技有限公司 | 一种多极化大地电磁测深方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9784882B2 (en) * | 2012-07-04 | 2017-10-10 | Genie Ip B.V. | Mapping hydrocarbon liquid properties of a kerogencontaining source rock |
-
2020
- 2020-09-01 CN CN202010905975.4A patent/CN111983701B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102419456A (zh) * | 2011-06-30 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 瞬变电磁测深数据的直接时间域处理方法 |
CN103389514A (zh) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | 一种大地电磁去噪装置及大地电磁去噪方法 |
CN104102814A (zh) * | 2014-06-11 | 2014-10-15 | 中国科学院地质与地球物理研究所 | 一种基于大地电磁数据反演电阻率和磁化率的方法及系统 |
CN105911603A (zh) * | 2016-05-04 | 2016-08-31 | 湖南科技大学 | 基于天然电场的四维物探方法 |
CN109765627A (zh) * | 2019-03-08 | 2019-05-17 | 中南大学 | 适用于山地的平面波电磁测深方法 |
CN110488365A (zh) * | 2019-09-11 | 2019-11-22 | 湖南五维地质科技有限公司 | 一种多极化大地电磁测深方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111983701A (zh) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2007295028B2 (en) | Rapid inversion of electromagnetic reconnaissance survey data | |
AU2007277410B2 (en) | Method for determining physical properties of structures | |
Nelson | Calculation of the magnetic gradient tensor from total field gradient measurements and its application to geophysical interpretation | |
WO2009146041A1 (en) | Constructing a reduced order model of an electromagnetic response in a subterranean structure | |
CN110488365B (zh) | 一种多极化大地电磁测深方法 | |
CN110007357B (zh) | 一种航空tem和航空mt联合反演方法 | |
CN111239827B (zh) | 基于局部相似系数的三维地震数据多次波压制方法 | |
CN104656156A (zh) | 音频大地电磁测深三维采集资料的磁参考处理方法 | |
Wang et al. | Forward modeling and inversion of tensor CSAMT in 3D anisotropic media | |
CN114036745A (zh) | 各向异性介质三维大地电磁正演方法及系统 | |
Qiu et al. | 3-D full-wave inversion of helicopter transient electromagnetic data in frequency domain | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
Yin et al. | Calculation of the magnetic gradient tensor from total magnetic anomaly field based on regularized method in frequency domain | |
CN111983701B (zh) | 横向大地电磁波测深方法及其测深装置 | |
CN108376204B (zh) | 一种基于地下粗糙介质模型的电磁广义趋肤深度计算方法 | |
Yan et al. | An analytical method to estimate shadow and source overprint effects in CSAMT sounding | |
US11816401B2 (en) | Providing for uncertainty in non-linear inversions of geophysical data | |
Utada et al. | ACTIVE system for monitoring volcanic activity: A case study of the Izu-Oshima Volcano, Central Japan | |
Rozier et al. | Inversion of a cylindrical vibrating body in shallow water from aspect-limited data using filtered SVD and the L-curve | |
KR101290332B1 (ko) | 장파장 속도 모델링에 의한 지하 구조의 영상화 장치 및 방법 | |
CN113064211B (zh) | 一种基于海洋磁力异常小波变换的居里面深度计算方法 | |
Chang et al. | Lateral variations of Rayleigh-wave dispersions in the Philippine Sea region | |
Plotkin et al. | Recording the vertical component of the electric field in magnetotelluric sounding | |
Vasilyev et al. | The method of the calibration of the geoelectric systems of the geodynamic control | |
Voronovich et al. | Horizontal refraction modal tomography of the ocean with mode interactions |
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 |