CN103630941A - 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 - Google Patents
长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 Download PDFInfo
- Publication number
- CN103630941A CN103630941A CN201310035114.5A CN201310035114A CN103630941A CN 103630941 A CN103630941 A CN 103630941A CN 201310035114 A CN201310035114 A CN 201310035114A CN 103630941 A CN103630941 A CN 103630941A
- Authority
- CN
- China
- Prior art keywords
- pseudorandomcode
- observation
- signal
- work area
- ground
- 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.)
- Granted
Links
Images
Abstract
本发明公开了一种长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法,该系统包括发射机、与发射机连接的两个接地电极A和B以及分布在整个工区中的多条测线。所述方法包括:采用伪随机编码信号作为激励场波形;将整个工区沿着测线连续分区,每个分区选择一个最佳伪随机编码信号的时钟频率;在所有观测点上阵列式同步观测大地同向电场,得到所有观测点的相干性数据;基于相干性数据和同步记录的激励场波形得到大地脉冲响应;基于大地脉冲响应得到整个工区的电阻率填图。本发明对工区进行非均匀分区,不同分区采用不同的伪随机编码时钟频率,使得所有观测点都有良好的信噪比和分辨率。本发明能够快速高效准确的进行大规模深部勘探和油气资源勘探。
Description
技术领域
本发明涉及长线源瞬变电磁探测技术、阵列式接收瞬变电磁法技术领域,特别涉及一种长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法。
背景技术
改革开放以来,随着国家政策的良好导向和科学技术的发展促进,我国经济得到全面快速的发展。与此同时,能源消耗更是和经济增长成比例的增加。我国作为全球最大的发展中国家,亟待解决对经济发展掣肘的能源问题。然而国际大环境给我国的能源问题的解决带来极大挑战。数次的中东战争使得国际油价被美国掌控,给我国的油气进口增加了高昂的成本;阿富汗、伊拉克和伊朗等国作为我国的主要油气进口国,要么被美国武力控制要么遭到美国的经济制裁和封锁,减少了我国能源进口的渠道。另一方面,西方等资源勘探强国盯住了我国亟需解决能源问题的迫切性,西方各国地球物理勘探公司纷纷不约而同的提高地球物理勘探设备和勘探服务的价格,提高价格达到几倍到几十倍不等,少数先进勘探技术更是对我国实行禁运封锁。
基于上述情况,在我国十二五战略规划重视下,以“地壳探测”项目为代表的资源探测项目陆续深入开展。其中,油气勘探可谓是重中之重,是关系到我国国民经济的战略问题。地震方法是有效勘探油气资源的传统方法,但对于地表及地质条件复杂的区域,却无能为力。而我国西北部等油气储量丰富的地区,多数是地质条件复杂且人迹罕至的区域,对于这些区域,电性源瞬变电磁法特别是长线源瞬变电磁法能从不同于地震方法的角度反应油气层的物理特性。目前电性源瞬变电磁法已经成为俄罗斯等国家油气勘探的主要手段。
长偏移距瞬变电磁测深法(LOTEM)是电性源瞬变电磁法的典型代表,LOTEM在设计之初就是为了油气勘探和深部地壳研究,由于采用大偏移距,导致信号较弱,信号的信噪比不够高。英国爱丁堡大学的MTEM方法在LOTEM方法的基础上,对接收方式和发射信号做了改进,可有效的进行油气勘查。然而MTEM方法在工作中,一次发射只有一条测线接收信号,在进行大规模油气探勘时,效率不够高。同时,由于受到测线几乎均匀分段的限制,MTEM方法不能兼顾众多观测点的信噪比和分辨率;另一方面,把电性源作为电偶极子近似处理,又对方法的有效性增加了限制。
发明内容
本发明的目的是为了克服上述现有技术中的不足,提出一种长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法。本发明对工区进行偏移距非均匀分区,不同分区采用不同的伪随机编码时钟频率,同时工区多条测线的所有观测点都有较好的信噪比和分辨率。本发明能够快速高效准确的进行大规模深部勘探和油气资源勘探。
本发明同步记录伪随机编码发射信号和同向电场分量电压,对二者进行反卷积处理得到大地脉冲响应值。采用长伪随机编码信号并周期重复发射,信噪比提高40dB以上;同时扩宽频带,加大了探测深度。对工区进行偏移距非均匀分区,不同分区采用不同的伪随机编码时钟频率,保证了所有观测点的信号都有较好的信噪比和分辨率,更有效的提高了工作效率,能实现整个工区的快速电阻率填图。
为了实现上述目的,本发明采用下述技术方案:
根据本发明的一方面,提出一种长线源伪随机编码发射-阵列式接收瞬变电磁系统,该系统包括发射模块和接收模块,其中,
所述发射模块包括发射机和与所述发射机连接的两个接地电极A和B;
所述发射机用于产生伪随机编码电流信号,并通过所述两个接地电极A、B将所述伪随机编码电流注入大地;
所述接收模块包括分布在整个工区中的多条测线,其中每条测线上设有若干个接收机,每个接收机作为一个观测点,用于通过与其连接的两个接地电极C、D来观测大地的AB向电场。
根据本发明的另一方面,提出一种长线源伪随机编码发射-阵列式接收瞬变电磁方法,该方法包括以下步骤:
步骤1,周期性重复发射伪随机编码电流信号,并同步记录激励场电流波形和工区中所有观测点的接收数据;
步骤2,根据观测点的最低可分辨电平,选择一个兼顾信号的信噪比和分辨率的最佳的伪随机编码电流信号的时钟频率,具有同一时钟频率的区域作为一个分区;
步骤3,对于每个分区内的所有观测点,阵列式同步观测平行于测线方向的感应电场,得到所有观测点上与发射的所述伪随机编码电流信号具有相干性的观测数据;
步骤4,采用反卷积信号处理方法对相干性观测数据和同步记录的激励场电流波形进行数据预处理,得到具有地电信息的大地脉冲响应;
步骤5,使用所述大地脉冲响应来进行反演,得到大地的地电信息,从而得到整个工区的电阻率填图。
本发明的有益效果是:1.采用长线源伪随机编码发射,有效的提高了观测数据的信噪比,使得本发明方法在深部勘探和油气勘探方面的能力进一步提高;2.扩宽了发射信号的频带,加大勘探深度;3.整个工区的非均匀分区并采用不同的伪随机编码时钟频率,兼顾所有测线上所有观测点的信噪比和分辨率;4.阵列式同步观测,使得工区所有测点能同步进行观测,大大的提高了施工效率,节约了时间和成本。
附图说明
图1是本发明长线源伪随机编码发射-阵列式接收瞬变电磁系统布置图;
图2是本发明长线源伪随机编码发射-阵列式接收瞬变电磁方法流程图;
图3是根据本发明一实施例的伪随机编码电流信号示意图;
图4是根据本发明一实施例的伪随机编码信号频谱示意图;
图5是根据本发明一实施例的伪随机编码信号自相关函数示意图;
图6是根据本发明一实施例得到的大地脉冲响应示意图;
图7是根据本发明一实施例的阵列式接收工区非均匀分区及对应时钟频率布置示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明中的阵列式,与以往在其它电磁法、地震法和雷达天线等出现的阵列式方法有本质的区别。以往的阵列式方法只是一些在空间有一定规则分布的点阵,各个数据点或者观测点相互独立的进行数据采集;而本发明中的阵列式是相干性阵列式接收,数据在阵列内的各观测点具有相干性,可以在同一时间完成工区所有观测点的同向电场和激励场电流波形的数据采集。相邻观测点的数据在频率域都有相干性,这给数据处理带来很大益处,可显著提高信号质量,这是本发明采用的阵列式接收系统的数据处理方法的特点。
图1是本发明长线源伪随机编码发射-阵列式接收瞬变电磁系统布置图,如图1所示,该系统包括:发射模块和接收模块,其中,
所述发射模块包括发射机和与所述发射机连接的两个接地电极A和B;
所述发射机用于产生伪随机编码电流信号,并通过所述两个接地电极A、B将所述伪随机编码电流注入大地;
所述接收模块包括分布在整个工区中的多条(比如九条)测线,其中每条测线上设有若干个接收机,每个接收机作为一个观测点,用于通过与其连接的两个接地电极C、D来观测大地的AB向电场;
其中,相邻的接收机共用一个接地电极。另外,根据整个工区的实际地质条件,多条测线是大致平行且等间距分布的,并且以两个接地电极A、B形成的发射线为对称轴。
图2是本发明长线源伪随机编码发射-阵列式接收瞬变电磁方法流程图,如图2所示,该方法包括以下步骤:
步骤1,周期性重复发射伪随机编码电流信号,并同步记录激励场电流波形和工区中所有观测点的接收数据,时变的伪随机编码电流信号通过接地良好的电极注入大地,从而激励大地产生电磁感应;
所述伪随机编码电流信号的示意图如图3所示,其频谱示意图如图4所示,图4中,所述伪随机编码电流信号的时钟频率为100Hz。
步骤2,根据观测点的最低可分辨电平,选择一个兼顾信号的信噪比和分辨率的最佳的伪随机编码电流信号的时钟频率,具有同一时钟频率的区域作为一个分区,其中,该时钟频率与发射机与观测点之间的距离,即偏移距的二次方成反比;
由于不同分区采用不同的最佳时钟频率,从而能够进一步提高观测数据的信号质量和分辨率。
另外,根据伪随机编码的特性,在时钟频率固定时,增加码长,数据处理后的信噪比会增大,因此可根据实际需要选择伪随机编码的码长以增大信噪比,具体地,可通过调节编码的阶数来选择码长:若伪随机编码的阶数为N,则对应的码长为2N-1,比如若伪随机编码阶数的调节范围是8阶到18阶,则对应的码长为256到262144。同时,在周期性重复发射伪随机编码电流信号时,还可根据实际需要选择其结尾的码元游程,所述码元游程的调节范围为1到N。
步骤3,对于每个分区内的所有观测点,阵列式同步观测平行于测线方向的感应电场,得到所有观测点上与发射的所述伪随机编码电流信号具有相干性的观测数据;
所述伪随机编码电流信号的自相关函数示意图如图5所示,其中,所述伪随机编码电流信号的码长为512。
步骤4,采用反卷积等信号处理方法对相干性观测数据和同步记录的激励场电流波形进行数据预处理,得到具有地电信息的大地脉冲响应,所述大地脉冲响应如图6所示,图中,tp表示峰值时刻,4tp、15tp和30tp分别是不同的观测结束时刻;
该步骤中,将同步记录的激励场电流(即发射的所述伪随机编码电流)信号记为i(t),将一定偏移距处的观测电压信号记为V(t),则通过对i(t)和V(t)数据的进行反卷积就可以获得具有地电信息的大地脉冲响应g(t):
V(t)=i(t)*g(t)+n(t) (1)
其中:i(t)为发射电流信号,g(t)为大地脉冲响应,V(t)是观测电压信号,n(t)为没有相关性的噪声;星号“*”表示卷积操作。
上述反卷积的步骤可描述为:
1)将(1)式转换到频域:
V=I·G+n (2)
其中,V、I、G和n分别为V(t)、I(t)、G(t)和n(t)对应的频域信号;
2)对上式两边同时乘以所述激励场电流波形的复共轭,得到:
I*V=I*I·G+I*n (3)
3)理论上,将I*I转换到时域上得到的是类冲激信号δ(t),但由于接收数据的长度有限及噪声的干扰,所述类冲激信号的长度不是无限窄,而是和发射电流信号的记录长度一致,因此将类冲激信号记为 将(3)式转换到时间域,得到:
其中,F-1(I*V)表示对于I*V进行反傅氏变换;
然后去除上式中的不相关噪声n1(t)(I*n对应的时域噪声),则得到的即是包含地电信息的大地脉冲响应,也就是是本发明要恢复的目标,其中,表明:时间越短,恢复得到的大地脉冲响应越接近真实值,而这是由伪随机编码的时钟频率决定的,同时也是工区分区选择伪随机编码时钟频率的下限的依据,这使得在选择伪随机编码时钟频率时,为了提高观测数据的信号质量,会降低时钟频率,但为了使大地脉冲响应接近真实值,又需要为时钟频率设置一个选择下限。
步骤5,使用所述大地脉冲响应来进行反演,得到大地的地电信息,从而得到整个工区的电阻率填图。
图7是根据本发明一实施例的阵列式接收工区非均匀分区及对应时钟频率布置示意图,在该实施例中,设置伪随机编码信号最大发射电流为70A;伪随机编码码长为4096;结尾码元游程为8;发射电极距为500m;接收电极距为200m;观测点间距为200m;预期探测深度d为2km;最大偏移距rmax为8km;最小偏移距rmin为3km;接收机最小分辨电平为0.2μV;接收阵列具有9条测线,与发射线呈对称分布,测线间距为500m;测线5与发射线同在x轴上,测线1在距离x轴2km处,测线9在另一侧距离x轴2km处;测线1-4和测线9-6的参数设置相同;测线1-4的非均匀分区如表1所示,整个工区的非均匀分区如图7所示;由于测线4和测线5在偏移距8km处的测点的接收信号低于接收机最小分辨电平,因此去除这两个测点。
表1工区分区和相应的时钟频率fs,其中,Vmin是各分区接收信号最小值,分区的单位为km
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种长线源伪随机编码发射-阵列式接收瞬变电磁系统,其特征在于,该系统包括发射模块和接收模块,其中,
所述发射模块包括发射机和与所述发射机连接的两个接地电极A和B;
所述发射机用于产生伪随机编码电流信号,并通过所述两个接地电极A、B将所述伪随机编码电流注入大地;
所述接收模块包括分布在整个工区中的多条测线,其中每条测线上设有若干个接收机,每个接收机作为一个观测点,用于通过与其连接的两个接地电极C、D来观测大地的AB向电场。
2.根据权利要求1所述的系统,其特征在于,相邻的接收机共用一个接地电极。
3.根据权利要求1所述的系统,其特征在于,所述多条测线大致平行且等间距分布。
4.根据权利要求1所述的系统,其特征在于,所述多条测线以两个接地电极A、B形成的发射线为对称轴。
5.一种长线源伪随机编码发射-阵列式接收瞬变电磁方法,其特征在于,该方法包括以下步骤:
步骤1,周期性重复发射伪随机编码电流信号,并同步记录激励场电流波形和工区中所有观测点的接收数据;
步骤2,根据观测点的最低可分辨电平,选择一个兼顾信号的信噪比和分辨率的最佳的伪随机编码电流信号的时钟频率,具有同一时钟频率的区域作为一个分区;
步骤3,对于每个分区内的所有观测点,阵列式同步观测平行于测线方向的感应电场,得到所有观测点上与发射的所述伪随机编码电流信号具有相干性的观测数据;
步骤4,采用反卷积信号处理方法对相干性观测数据和同步记录的激励场电流波形进行数据预处理,得到具有地电信息的大地脉冲响应;
步骤5,使用所述大地脉冲响应来进行反演,得到大地的地电信息,从而得到整个工区的电阻率填图。
6.根据权利要求5所述的方法,其特征在于,所述时钟频率与发射机与观测点之间的距离,即偏移距的二次方成反比。
7.根据权利要求5所述的方法,其特征在于,为了进一步提高观测数据的信噪比,可通过调节编码的阶数来选择伪随机编码的码长:若伪随机编码的阶数为N,则对应的码长为2N-1。
8.根据权利要求5所述的方法,其特征在于,所述伪随机编码的码元游程可调,其调节范围为1到N。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310035114.5A CN103630941B (zh) | 2013-01-30 | 2013-01-30 | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310035114.5A CN103630941B (zh) | 2013-01-30 | 2013-01-30 | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630941A true CN103630941A (zh) | 2014-03-12 |
CN103630941B CN103630941B (zh) | 2016-08-24 |
Family
ID=50212188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310035114.5A Active CN103630941B (zh) | 2013-01-30 | 2013-01-30 | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630941B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929455A (zh) * | 2016-04-18 | 2016-09-07 | 中国科学院地质与地球物理研究所 | 一种多通道瞬变电磁法三维探测的装置和方法 |
CN104375194B (zh) * | 2014-11-10 | 2017-01-11 | 山东能源集团有限公司 | 深部开采矿井富水区电性源瞬变电磁探查方法 |
CN106842334A (zh) * | 2016-12-12 | 2017-06-13 | 中国石油天然气集团公司 | 一种电磁勘探方法和装置 |
CN107589461A (zh) * | 2017-09-04 | 2018-01-16 | 长江大学 | 一种基于双向编码的尖脉冲时域电磁深层探测方法 |
CN111580169A (zh) * | 2020-05-29 | 2020-08-25 | 河海大学 | 一种用于堰塞坝渗漏通道的抗干扰探测系统及方法 |
CN111694060A (zh) * | 2020-05-13 | 2020-09-22 | 东华理工大学 | 基于footprint技术的多道瞬变电磁反演方法、瞬变电磁勘测装置 |
CN113050175A (zh) * | 2021-03-08 | 2021-06-29 | 吉林大学 | 直升机航空电磁收录装置及发射源参数识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050237063A1 (en) * | 2001-09-07 | 2005-10-27 | Wright David A | Detection of subsurface resistivity contrasts with application to location of fluids |
CN101176017A (zh) * | 2005-03-14 | 2008-05-07 | 睦泰姆有限责任公司 | 真实幅度瞬时电磁系统响应测量 |
US20100201367A1 (en) * | 2009-02-11 | 2010-08-12 | Antoni Marjan Ziolkowski | Short-offset transient electromagnetic geophysical surveying |
CN102419452A (zh) * | 2011-06-30 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 电性源瞬变电磁法全场区探测方法 |
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
-
2013
- 2013-01-30 CN CN201310035114.5A patent/CN103630941B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050237063A1 (en) * | 2001-09-07 | 2005-10-27 | Wright David A | Detection of subsurface resistivity contrasts with application to location of fluids |
CN101176017A (zh) * | 2005-03-14 | 2008-05-07 | 睦泰姆有限责任公司 | 真实幅度瞬时电磁系统响应测量 |
US20100201367A1 (en) * | 2009-02-11 | 2010-08-12 | Antoni Marjan Ziolkowski | Short-offset transient electromagnetic geophysical surveying |
CN102419452A (zh) * | 2011-06-30 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 电性源瞬变电磁法全场区探测方法 |
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
Non-Patent Citations (2)
Title |
---|
向葵: "《长偏移距瞬变电磁法》", 《内蒙古石油化工》, vol. 33, no. 3, 31 March 2007 (2007-03-31), pages 205 * |
唐新功等: "《地堑地形对长偏移距瞬变电磁测深的影响研究》", 《工程地球物理学报》, vol. 1, no. 4, 31 August 2004 (2004-08-31), pages 313 - 316 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375194B (zh) * | 2014-11-10 | 2017-01-11 | 山东能源集团有限公司 | 深部开采矿井富水区电性源瞬变电磁探查方法 |
CN105929455A (zh) * | 2016-04-18 | 2016-09-07 | 中国科学院地质与地球物理研究所 | 一种多通道瞬变电磁法三维探测的装置和方法 |
CN105929455B (zh) * | 2016-04-18 | 2016-12-14 | 中国科学院地质与地球物理研究所 | 一种多通道瞬变电磁法三维探测的装置和方法 |
CN106842334A (zh) * | 2016-12-12 | 2017-06-13 | 中国石油天然气集团公司 | 一种电磁勘探方法和装置 |
CN107589461A (zh) * | 2017-09-04 | 2018-01-16 | 长江大学 | 一种基于双向编码的尖脉冲时域电磁深层探测方法 |
CN107589461B (zh) * | 2017-09-04 | 2019-11-05 | 长江大学 | 一种基于双向编码的尖脉冲时域电磁深层探测方法 |
CN111694060A (zh) * | 2020-05-13 | 2020-09-22 | 东华理工大学 | 基于footprint技术的多道瞬变电磁反演方法、瞬变电磁勘测装置 |
CN111580169A (zh) * | 2020-05-29 | 2020-08-25 | 河海大学 | 一种用于堰塞坝渗漏通道的抗干扰探测系统及方法 |
CN113050175A (zh) * | 2021-03-08 | 2021-06-29 | 吉林大学 | 直升机航空电磁收录装置及发射源参数识别方法 |
CN113050175B (zh) * | 2021-03-08 | 2022-05-13 | 吉林大学 | 直升机航空电磁收录装置及发射源参数识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103630941B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103630941A (zh) | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 | |
CN102721982B (zh) | 基于扩频编码技术的地面电磁勘探方法及其探测系统 | |
CN102426393B (zh) | 电法勘探方法 | |
CN101258422B (zh) | 用多个发射器执行受控源电磁探测的方法 | |
CA2671434C (en) | Systems and methods for measuring sea-bed resistivity | |
CN103207413B (zh) | 电法勘探装置及系统 | |
US7667464B2 (en) | Time segmentation of frequencies in controlled source electromagnetic (CSEM) applications | |
CN202330736U (zh) | 一种电法勘探装置 | |
CN105158808A (zh) | 一种浅海瞬变电磁海空探测及其解释方法 | |
CN103837897B (zh) | 煤矿井下槽波精细探测方法 | |
CN204556851U (zh) | 一种海底电缆全方位地震采集观测系统 | |
CN1239922C (zh) | 人工源时间频率电磁测深方法 | |
Li et al. | Application of grounded electrical source airborne transient electromagnetic (GREATEM) system in goaf water detection | |
CN107589461A (zh) | 一种基于双向编码的尖脉冲时域电磁深层探测方法 | |
Yang et al. | Distributed wide field electromagnetic method based on high-order 2n sequence pseudo random signal | |
CN110376651B (zh) | 基于水平双极电流源的时频电磁装置及地球物理勘探方法 | |
CN105891892A (zh) | 全区测量电流源频率域磁场垂直分量的装置及方法 | |
CN102636807A (zh) | 电磁式可控震源地震信号检测方法 | |
CN103605158A (zh) | 一种最大炮检距确定方法及装置 | |
WO2015026650A1 (en) | Systems and methods for remote electromagnetic exploration for mineral and energy resources using stationary long-range transmitters | |
Duncan et al. | SALTMAP—Airborne EM for the environment | |
Engelmark et al. | Towed streamer EM comes of age | |
Ingerov et al. | Non-grounded surface electroprospecting technique | |
CN103336302B (zh) | 基于高次余弦幅度加权的地震波束形成方法 | |
CN104502986B (zh) | 物探激电测深数据层析法处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |