CN108897051B - 一种高分辨率的地面核磁共振成像方法 - Google Patents
一种高分辨率的地面核磁共振成像方法 Download PDFInfo
- Publication number
- CN108897051B CN108897051B CN201810459074.XA CN201810459074A CN108897051B CN 108897051 B CN108897051 B CN 108897051B CN 201810459074 A CN201810459074 A CN 201810459074A CN 108897051 B CN108897051 B CN 108897051B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- wave field
- nuclear magnetic
- field
- seismic wave
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 title claims abstract description 11
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 44
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 22
- 238000003384 imaging method Methods 0.000 claims abstract description 21
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 230000005012 migration Effects 0.000 claims abstract description 11
- 238000013508 migration Methods 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000002595 magnetic resonance imaging Methods 0.000 claims abstract description 5
- 230000004044 response Effects 0.000 claims description 19
- 230000010354 integration Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 7
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims 1
- 230000001427 coherent effect Effects 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
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/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种高分辨率的地面核磁共振成像方法，解决了传统地面核磁共振成像方法在进一步提高地下水解释精度方面遇到的瓶颈问题。应用阻尼法最小二乘法，求解核磁共振电磁场与地震波动场方程，实现波场变换，得到测线上各接收线圈拟地震波场离散数据；通过对测线上各接收线圈拟地震波场离散数据求取反褶积，消除波场变换的波形展宽效应；基于相关叠加原理求取合成孔径范围内各点互相关系数，实现相关点叠加，提高探测信噪比，获取测线上的合成孔径虚拟地震波场合成值；采用克希霍夫偏移成像理论，求解拟地震波波动方程，实现地下含水结构偏移成像。本发明规避了传统核磁共振数据解释方法含水层边缘分辨率差等缺点。
Description
技术领域
本发明涉及一种高分辨率的地面核磁共振成像方法，属于地面核磁共振地球物理勘探方法技术领域。
背景技术
地面核磁共振(Surface Nuclear Magnetic Resonance，简称SNMR)是近年来新兴的地球物理方法，主要应用于浅层地下水勘探及水源性地质灾害预测。探测时，通过地面铺设线圈向地下发射氢质子自旋频率的交流场，激发地下水产生核磁共振信号，并通过反演、成像等过程判定含水层含水量及孔隙度等信息，具有直接非入侵探测、定量定性及信息量丰富等优点。
然而，现有SNMR技术信号十分微弱，仅为纳伏级别，虽然能够实现含水位置及含量信息的获取，但探测信噪比及含水层边缘分辨率较低，难以实现精确的含水范围圈定。为此，探索高分辨率的地面核磁共振地下水成像方法具有重要现实意义。
发明内容
本发明的目的在于针对现有核磁共振成像技术的不足，提供一种高分辨率的地面核磁共振成像方法。包括应用虚拟波场变换技术，将地面核磁共振探测中的电磁接收信号转换为拟地震波，有效突出传统地面核磁共振解释方法难以分辨的含水层边缘；消除波场变换后虚拟子波存在的波形展宽现象，进一步提高核磁共振成像空间分辨能力；利用合成孔径方法对不同接收线圈波场转换后的拟地震信号叠加，从而提高探测信噪比。求解拟地震波波动方程，实现高精度地下含水结构成像。
本发明是通过以下技术方案方式实现的：
一种高分辨率的地面核磁共振成像方法，包括：应用阻尼法最小二乘法，求解核磁共振电磁场与地震波动场方程，实现波场变换，得到测线上各接收线圈拟地震波场离散数据；通过对测线上各接收线圈拟地震波场离散数据求取反褶积，消除波场变换的波形展宽效应；基于相关叠加原理求取合成孔径范围内各点互相关系数，实现相关点叠加，提高探测信噪比，获取测线上的合成孔径虚拟地震波场合成值；采用克希霍夫偏移成像理论，求解拟地震波波动方程，实现地下含水结构偏移成像。
实现地下含水结构偏移成像。
一种高分辨率的地面核磁共振成像方法，包括以下步骤：
a、建立核磁共振电磁响应场波场变换关系；
通过SNMR仪器控制发射线圈产生激发电磁场，激发地下水中氢质子吸收能量并发生能级跃迁；激发结束后，氢质子释放能量并产生核磁共振响应电磁场；测线上的多个接收线圈在接收感应场作用下，产生核磁共振信号：
其中，t为时间，r为对应的空间位置，HNMR(r,t)为核磁共振响应电磁场，μ为磁导率；
核磁共振响应电磁场与虚拟波场满足：
其中，U(r,τ)为以波速传播的虚拟波场，自变量τ是时间平方根的量纲；
将公式(2)改写为数值积分形式为：
其中m为采样个数，n为数值积分时积分区间剖分个数，pj为积分步长；
b、对公式(3)离散化得到对应公式(3)的线性方程组：
H＝AU, (4)
其中，U＝(u1,…,uj,…,un)T为离散化的虚拟子波，H＝(h1,…,hi,…,hm)T为离散化的核磁共振响应场，A＝(Aij)为m×n阶矩阵，且m≥n，
采用阻尼法最小二乘法求解公式(4)，即可得到对应核磁共振响应场H的拟地震波场U。
c、对步骤b得到的测线上各接收线圈拟地震波场离散数据U求取反褶积，消除波场变换的波形展宽效应。
U(r,τ)为步骤b求解得到的虚拟地震子波，通过反褶积滤波因子获取宽度，得到压缩后的新的子波：
其中，h(t)为反褶积滤波因子，可通过最小平方反褶积求得。
d、在测线上选定合成孔径范围-N,…,i,…,N，第i点为中心点，U(ri,τ)为该点的虚拟波场信号，ri为i点到i-N,…,i+N点的距离，τ为相对时移量。求各点到中心点的归一化互相关系数，并采用相关叠加方法对测线上各点的拟地震波场进行合成，得到中心点虚拟地震波场合成值；
e、在测线上依次移动，得到中心点外其他点处的虚拟地震波场合成值；
f、获取测线上的合成孔径虚拟地震波场值后，基于克希霍夫偏移成像理论，求解核磁共振虚拟波场波动方程，完成地下水三维偏移成像。
本发明的有益效果在于：本专利提出了一种高分辨率的地面核磁共振成像方法，解决了传统地面核磁共振成像方法在进一步提高地下水解释精度方面遇到的瓶颈问题。与现有技术相比，本发明基于核磁共振响应扩散场与拟地震波动场之间的数学积分变换实现含水层的高精度成像，规避了传统核磁共振数据解释方法含水层边缘分辨率差等缺点，对地面核磁共振技术进一步的应用及推广具有一定意义。
附图说明
图1地面核磁共振合成孔径示意图；
图2克希霍夫成像闭区域与边界示意图。
1SNMR仪器，2发射线圈，3第1接收线圈，4第2接收线圈，5第3接收线圈，6第4接收线圈，7第5接收线圈，8第6接收线圈，9第7接收线圈，10含水层，13测线。
具体实施方式
一种高分辨率的地面核磁共振成像方法，主要适用于赋水结构复杂的地下水探测成像。其利用核磁共振信号满足的扩散方程与地震信号满足的波动方程间的数学关系，将核磁共振响应电磁信号HNMR转化为拟地震波U，从而实现不同探测方向数据的叠加，提高信噪比。由于地震波场对几何截面更为敏感，所以合成孔径成像能够有效增强含水结构边缘的分辨率，对复杂含水结构的精确成像和准确解释水文地质参数，实现地下赋水结构复杂多变地区水资源的勘查具有重要意义。
下面结合附图和实施例对本发明作进一步详细说明：
一种高分辨率的地面核磁共振成像方法，包括以下顺序和步骤：
a、建立核磁共振电磁响应场波场变换关系；
如图1所示，通过SNMR仪器1控制发射线圈2产生激发电磁场，激发含水层10地下水中氢质子吸收能量并发生能级跃迁；激发结束后，氢质子释放能量并产生核磁共振响应电磁场；测线13上的第1接收线圈3，第2接收线圈4，第3接收线圈5，第4接收线圈6，第5接收线圈7，第6接收线圈8，第7接收线圈9在接收感应场作用下，产生核磁共振信号：
其中，t为时间，r为对应的空间位置，HNMR(r,t)为核磁共振响应电磁场，μ为磁导率；
核磁共振响应电磁场与虚拟波场满足：
其中，U(r,τ)为以波速传播的虚拟波场，自变量τ是时间平方根的量纲；
将公式(2)改写为数值积分形式为：
其中m为采样个数，n为数值积分时积分区间剖分个数，pj为积分步长；
b、对公式(3)离散化得到对应公式(3)的线性方程组：
H＝AU, (4)
其中，U＝(u1,…,uj,…,un)T为离散化的虚拟子波，H＝(h1,…,hi,…,hm)T为离散化的核磁共振响应场，A＝(Aij)为m×n阶矩阵，且m≥n，
采用阻尼法最小二乘法求解方程组(4)，即可得到对应核磁共振响应场H的拟地震波场U。
c、对步骤b得到的测线13上各接收线圈拟地震波场离散数据U求取反褶积，消除波场变换的波形展宽效应。
U(r,τ)为步骤b求解得到的虚拟地震子波，通过反褶积滤波因子获取宽度，得到压缩后的新的子波：
其中，h(t)为反褶积滤波因子，可通过最小平方反褶积求得。
d、在测线13上选定合成孔径范围：第1接收线圈3为第i-3点，第2接收线圈4为第i-2点，第3接收线圈5为第i-1点，第5接收线圈7为第i+1点，第6接收线圈8为第i+2点，第7接收线圈9为第第1接收线圈3为第i-3点，第2接收线圈4为第i-2点，第3接收线圈5为第i-1点，第5接收线圈7为第i+1点，第6接收线圈8为第i+2点，第7接收线圈9为第i+3点，第4接收线圈6为中心点i，U(ri,τ)为该点的虚拟波场信号，ri为i点到i-3,…,i+3内某点的距离，τ为相对时移量。求各接收线圈到第4接收线圈6的归一化互相关系数，并采用相关叠加方法对测线13上的各接收线圈拟地震波场进行合成，得到第4接收线圈6处的虚拟地震波场合成值；
e、在测线13上依次移动，得到第1接收线圈3，第2接收线圈4，第3接收线圈5，第5接收线圈7，第6接收线圈8，第7接收线圈9处的虚拟地震波场场合成值；
f、获取测线13上的合成孔径虚拟地震波场值后，基于克希霍夫偏移成像理论，求解核磁共振虚拟波场波动方程，完成地下水三维偏移成像。
步骤b中所述求解方法：
建立波场变换求解目标函数：
其中，ε2为阻尼因子，令公式(7)的梯度为0，可得
(ATA-ε2I)U＝ATH, (7)
进而有
U＝(ATA-ε2I)+ATH, (8)
其中，(ATA-ε2I)+为(ATA-ε2I)的Moore-Penrose广义逆。
步骤d中所述方法：
如图1，中心点i外的点分别与中心点做相关，其归一化互相关系数：
其中，m为每一测点的时间道数。
改变时移量τ找出最大相关系数ρmax(ri,τm)，此时的时移量称为最佳延时τm。
通过计算各个接收线圈对应点最大相关系数及最佳延时并应用相关计算得到最大相关系数作为权系数，各点权系数与波场值乘积即为叠加到中心点的拟地震波场合成值：
步骤f中所述方法：
如图2，设空间中闭区域具有闭合面Q，Q＝Q0+Q1，根据希尔霍夫理论，围绕震源所在闭合面Q上已知波动的位移位U'(x,y,z,t)及其导数，且这些值是连续的，则Q面外任意观测点M(x1,y1,z1)上由源引起的位移位U'的解：
其中，[U']为延迟位，其表达式为
虚拟波长偏移是获取记录的逆过程，已知地表观测点数据，需确定反射界面上作为二次虚拟波源的空间位置。令U'(x,y,z,t)＝u(x,y,z,-t)，则u(x,y,z,-t)仍然满足方程(11)，且u(x,y,z,-t)可看做时间的倒退，即把反射界面各点等效为上行波源。将接收点信号拟时间方向还原到二次波源，以寻找反射界面的波场函数，确定反射界面位置。令G(x,y,z0,t)为自激自收的虚拟波动信号，是地下反射界面产生的二次源激发的波长g(x,y,z,t)在地表z0上的值。由公式(11)可知
即为波场的向下延拓。由于实测数据为离散数据，所以采用边界元技术进行克希霍夫积分离散。对曲面边界Q0进行剖分，则公式(12)变化为
将上式变换为矩阵形式，带入已知地表波场函数值、波场法相导数及波场速度，即可求解出地下任意点处波场函数值，实现地下含水结构偏移成像。
Claims (5)
1.一种高分辨率的地面核磁共振成像方法，其特征在于，包括如下的步骤:
应用阻尼法最小二乘法，求解核磁共振电磁场与地震波动场方程，实现波场变换，得到测线上各接收线圈拟地震波场离散数据；
通过对测线上各接收线圈拟地震波场离散数据求取反褶积，消除波场变换的波形展宽效应；
基于相关叠加原理求取合成孔径范围内各点互相关系数，实现相关点叠加，提高探测信噪比，获取测线上的合成孔径虚拟地震波场合成值；
采用克希霍夫偏移成像理论，求解拟地震波波动方程，实现地下含水结构偏移成像。
2.按照权利要求1所述的一种高分辨率的地面核磁共振成像方法，其特征在于，用阻尼法最小二乘法，求解核磁共振电磁场与地震波动场方程，实现波场变换包括步骤a：
建立核磁共振电磁响应场波场变换关系；
通过SNMR仪器控制发射线圈产生激发电磁场，激发地下水中氢质子吸收能量并发生能级跃迁；激发结束后，氢质子释放能量并产生核磁共振响应电磁场；测线上的多个接收线圈在接收感应场作用下，产生核磁共振信号：
其中，t为时间，r为对应的空间位置，HNMR(r,t)为核磁共振响应电磁场，μ为磁导率；
核磁共振响应电磁场与虚拟波场满足：
其中，U(r,τ)为以波速传播的虚拟波场，自变量τ是时间平方根的量纲；
将公式(2)改写为数值积分形式为：
其中m为采样个数，n为数值积分时积分区间剖分个数，pj为积分步长。
3.按照权利要求2所述的一种高分辨率的地面核磁共振成像方法，其特征在于，还包括步骤b：对公式(3)离散化得到对应公式(3)的线性方程组：
H＝AU, (4)
其中，U＝(u1,…,uj,…,un)T为离散化的虚拟子波，H＝(h1,…,hi,…,hm)T为离散化的核磁共振响应场，A＝(Aij)为m×n阶矩阵，且m≥n，
采用阻尼法最小二乘法求解公式(4)，得到对应核磁共振响应场H的拟地震波场U。
4.按照权利要求3所述的一种高分辨率的地面核磁共振成像方法，其特征在于，对步骤b得到的测线上各接收线圈拟地震波场U求取反褶积，消除波场变换的波形展宽效应；
U(r,τ)为步骤b求解得到的虚拟地震子波，通过反褶积滤波因子获取宽度，得到压缩后的新的子波：
其中，h(t)为反褶积滤波因子，通过最小平方反褶积求得。
5.按照权利要求1所述的一种高分辨率的地面核磁共振成像方法，其特征在于，基于相关叠加原理求取合成孔径范围内各点互相关系数，实现相关点叠加，提高探测信噪比，包括步骤d：在测线上选定合成孔径范围-N,…,i,…,N，第i点为中心点，U(ri,τ)为该点的虚拟波场信号，ri为i点到i-N,…,i+N点的距离，τ为相对时移量；求各点到中心点的归一化互相关系数，并采用相关叠加方法对测线上各点的拟地震波场进行合成，得到中心点虚拟地震波场合成值；
e、在测线上依次移动，得到中心点外其他点处的虚拟地震波场合成值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810459074.XA CN108897051B (zh) | 2018-05-15 | 2018-05-15 | 一种高分辨率的地面核磁共振成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810459074.XA CN108897051B (zh) | 2018-05-15 | 2018-05-15 | 一种高分辨率的地面核磁共振成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108897051A CN108897051A (zh) | 2018-11-27 |
| CN108897051B true CN108897051B (zh) | 2019-10-01 |
Family
ID=64342811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810459074.XA Active CN108897051B (zh) | 2018-05-15 | 2018-05-15 | 一种高分辨率的地面核磁共振成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108897051B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109633761B (zh) * | 2018-12-13 | 2021-05-28 | 吉林大学 | 基于小波变换模极大值法的磁共振信号工频噪声消减方法 |
| CN113504574B (zh) * | 2021-04-30 | 2022-07-01 | 青海省第三地质勘查院 | 基于归一化函数的电磁数据拟地震剖面显示方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104330828A (zh) * | 2014-10-27 | 2015-02-04 | 中国石油天然气股份有限公司 | 甜点储层的预测方法及预测装置 |
| WO2016094170A1 (en) * | 2014-12-10 | 2016-06-16 | Schlumberger Canada Limited | Formation property characteristic determination methods |
| CN106842318A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | 微观溶蚀孔隙地球物理二维表征确定方法和装置 |
| CN107167846A (zh) * | 2017-05-19 | 2017-09-15 | 吉林大学 | 地空结合多功能高精度快速地磁信息测量装置及测量方法 |
-
2018
- 2018-05-15 CN CN201810459074.XA patent/CN108897051B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104330828A (zh) * | 2014-10-27 | 2015-02-04 | 中国石油天然气股份有限公司 | 甜点储层的预测方法及预测装置 |
| WO2016094170A1 (en) * | 2014-12-10 | 2016-06-16 | Schlumberger Canada Limited | Formation property characteristic determination methods |
| CN106842318A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | 微观溶蚀孔隙地球物理二维表征确定方法和装置 |
| CN107167846A (zh) * | 2017-05-19 | 2017-09-15 | 吉林大学 | 地空结合多功能高精度快速地磁信息测量装置及测量方法 |
Non-Patent Citations (1)
| Title |
|---|
| "基于自适应遗传算法的MRS－TEM联合反演方法研究";万玲 等;《地球物理学报》;20131130;第56卷(第11期);第3728-3739页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108897051A (zh) | 2018-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Klotzsche et al. | Review of crosshole ground-penetrating radar full-waveform inversion of experimental data: Recent developments, challenges, and pitfalls | |
| RU2503035C2 (ru) | Оценивание свойств почвы с использованием волновых сигналов сейсмических поверхностных волн | |
| CN108037526B (zh) | 基于全波波场vsp/rvsp地震资料的逆时偏移方法 | |
| Nemeth et al. | Least-squares migration of incomplete reflection data | |
| CN102692645B (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
| Zhu et al. | Imaging diffraction points using the local image matrices generated in prestack migration | |
| AU2005294550A1 (en) | Amplitude preserving prestack migration method | |
| Zhou et al. | Migration velocity analysis and prestack migration of common-transmitter GPR data | |
| CN105158808A (zh) | 一种浅海瞬变电磁海空探测及其解释方法 | |
| US20100110830A1 (en) | Method of determining properties of the earth | |
| CN104280777A (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
| CN106443770A (zh) | 一种页岩气地质甜点的预测方法 | |
| CN108897051B (zh) | 一种高分辨率的地面核磁共振成像方法 | |
| Maurer et al. | Optimized experimental design in the context of seismic full waveform inversion and seismic waveform imaging | |
| Liu et al. | Linear prediction-based DOA estimation for directional borehole radar 3-D imaging | |
| CN114779355A (zh) | 基于发射电流全波形的地面瞬变电磁法反演方法及装置 | |
| Song et al. | Utilization of multimode surface wave dispersion for characterizing roadbed structure | |
| Wang et al. | Road sinkhole detection with 2D ambient noise tomography | |
| CN102385066B (zh) | 一种叠前地震定量成像方法 | |
| van der Kruk et al. | GPR full-waveform inversion, recent developments, and future opportunities | |
| CN104280703A (zh) | 一种基于脉冲梯度磁场的新型d-t2脉冲施加方法 | |
| Kumar et al. | Resolving phase wrapping by using sliding transform for generation of dispersion curves | |
| Liu et al. | Weighted diffraction-based migration velocity analysis of common-offset GPR reflection data | |
| CN105259573A (zh) | 一种地震沉积体信息的获取方法及装置 | |
| Su et al. | High-resolution seismic processing technique with broadband, wide-azimuth, and high-density seismic data—A case study of thin-sand reservoirs in eastern China |
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 |