CN110531416B - 一种基于三分量反射信号时频域极化参数的断层确定方法 - Google Patents
一种基于三分量反射信号时频域极化参数的断层确定方法 Download PDFInfo
- Publication number
- CN110531416B CN110531416B CN201910774252.2A CN201910774252A CN110531416B CN 110531416 B CN110531416 B CN 110531416B CN 201910774252 A CN201910774252 A CN 201910774252A CN 110531416 B CN110531416 B CN 110531416B
- Authority
- CN
- China
- Prior art keywords
- component
- main polarization
- fault
- reflection
- signals
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000005641 tunneling Effects 0.000 claims abstract description 16
- 229910052704 radon Inorganic materials 0.000 claims abstract description 5
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 3
- 101150072344 argA gene Proteins 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 7
- 239000003245 coal Substances 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于三分量反射信号时频域极化参数的断层确定方法，先建立观测***，在炮点激发地震波后，将地震仪采集到的各个地震波信号整理为三个分量的共炮点道集，在三个分量的共炮点道集中利用拉冬变换提取负速度同相轴，得到每个三分量检波器关于断层界面的反射三分量信号；利用获得的反射三分量信号，对上述各个检波器的反射三分量信号进行时频域主极化方向计算，求取瞬时极化倾角
与方位角φi，利用空间交汇关系确定虚震源点的位置。本发明能精确确定断层的具体倾角参数以及位置信息，从而保证超前探测断层的位置精度，为后续掘进提供数据支撑。
Description
技术领域
本发明涉及一种断层确定方法，具体是一种基于三分量反射信号时频域极化参数的断层确定方法。
背景技术
隧道的建设给交通网络的扩展带来极大的便利，而隧道在开掘过程中难以避免的会遇到安全问题。除去生产过程中不遵守规范导致的安全事故，更多的事故是未提前准确查明隧道前方的地质构造情况导致的，如掘进前方存在地质断层而仍然沿原路线掘进则会发生安全事故。地质断层的发育伴随着地层的失稳和岩性的突变，使原先连续渐变的地层在断层处产生畸变，影响隧道掘进的效率；另外某些断层形成后会与其他含水异常体形成导水通道，从而导致掘进到该处时发生涌水情况，甚至会对施工人员人身安全产生威胁。基于上述原因，目前针对断层的地球物理特点，利用地震超前探测手段进行探查，但是现有的反射地震隧道超前探测手段，采用求取速度曲线并偏移成像的步骤，这个步骤比较复杂，由于隧道空间的特性，导致速度求取难度较大，速度求取不正确导致断层位置确定不准、根据偏移结果推断反射界面导致的断层角度求取不准的问题；因此难以确定断层的具体倾角参数以及位置信息，最终导致超前探测断层的位置精度较低，无法满足对后续掘进的指导作用。因此，提出一种不依赖于速度模型的断层确定方法很重要。
专利号为：201310123575.8的中国发明专利提出了基于单一炮检对观测***的煤巷超前探测断失翼煤层的方法，利用煤层断点绕射波反演出断失翼煤层位置，该方法利用一个X、Z两分量接收点地震记录计算出断点的绕射波主极化方向，通过接收点的二维坐标及主极化方向形成一条射线，以接收点为圆心并以断点绕射波距离为半径形成一个圆弧，射线与圆弧相交即为断失翼煤层位置。与上述专利类似，201610316770.6的中国发明专利提出了一种突出煤层巷道小断层超前定性预报方法，通过顶板或底板的两个接收点的射线相交，确定断点位置。但是，上述方法都存在很大的局限性，主要体现在：上述方法对主极化方向确定的精度要求极高，一旦得出的方向略有误差则会导致一个射线与一个圆弧(或射线)相交得出的断失翼煤层位置偏差很大，然而实际地震记录复杂导致主极化方向的求解误差很大，特别是三维空间中两条射线之间可能不相交，从而不能得到交点。综上所述，现有的方法均不能达到断层的精细探测需求。
发明内容
针对上述现有技术存在的问题，本发明提供一种基于三分量反射信号时频域极化参数的断层确定方法，能精确确定断层的具体倾角参数以及位置信息，从而保证超前探测断层的位置精度，为后续掘进提供数据支撑。
为了实现上述目的，本发明采用的技术方案是：一种基于三分量反射信号时频域极化参数的断层确定方法，该方法的具体步骤为：
A、在隧道两帮根据现场施工条件选取其中一帮，在该帮上沿着隧道掘进方向在隧道腰线位置以等间距2m布设16个三分量检波器，由后向前依次编号为R1～R16，在三分量检波器R1的后方10m位置同一高度布置一个炮点，记为S1，在炮点向前激发地震波，16个三分量检波器将接收到的地震波信号传递给地震仪进行记录；
B、根据已知的炮点S1和三分量检波器R1～R16各个位置，建立观测***，并将地震仪采集到的各个地震波信号整理为三个分量的共炮点道集；
C、在三个分量的共炮点道集中利用拉冬变换提取负速度同相轴，得到每个三分量检波器关于断层界面的反射三分量信号；
D、利用步骤C中获得的反射三分量信号，对上述16个三分量检波器的反射三分量信号进行时频域主极化方向计算，求取主极化倾角
与主极化方位角φi，利用空间交汇关系确定断层反射界面倾角与位置，具体确定过程为：
Ⅰ、选取三分量检波器R1关于断层界面的反射三分量信号得出的主极化倾角
与主极化方位角φi，其步骤为：
①反射三分量信号采用解析信号Ai(t)来表示：
式中，
表示采样点优势频率，argAi(t)表示求取Ai(t)的幅角；
②构建协方差矩阵：
式中
表示t时刻的l分量时窗叠加振幅，第l分量在时窗内的均值μl，Tlh(t)表示t时刻
的自适应时窗长度，h为1、2、3；
③通过上述协方差矩阵求取的特征向量归一化后表示为[x1(t2),y1(t2),z1(t2)]；
④根据上述特征向量求取主极化倾角，其表示为：
其中
表示主极化方向与XOZ面之间的夹角，即主极化倾角；
⑤根据上述特征向量求取主极化方位角，可表示为：
其中
表示主极化方向与XOY面之间的夹角，即主极化方位角；
Ⅱ、确定该断层界面反射波信号传播的方向
其中(m n p)由主极化倾角
与主极化方位角φi确定，(xR yR zR)表示三分量检波器的空间坐标：
Ⅲ、重复步骤Ⅰ、Ⅱ，确定其他15个三分量检波器关于断层界面反射波信号传播的方向
Ⅳ、将各个三分量检波器的采样点优势频率的倒数取平均值，以该平均值的1/4为网格大小，并将三维坐标系网格化，将16个断层界面反射波信号传播的方向
延伸经过各个空间网格，并记录每个网格内经过的射线数目；
Ⅴ、统计各个网格内经过的射线数目，经过射线数目最多的网格即为所有向量的空间交汇点，确定空间交汇点的网格中心为虚震源点，在三维坐标系中做连接炮点与虚震源点的直线段；
Ⅵ、对连接炮点与虚震源点的直线段做垂直平分线，三维坐标系中根据该垂直平分线与隧道掘进中轴线之间夹角确定断层界面的倾角、方位角，该垂直平分线的位置与隧道掘进中轴线之间交点即为断层界面位置。
与现有技术相比，本发明先在巷道同一帮布设16个三分量检波器和一个炮点，建立观测***，在炮点激发地震波后，16个三分量检波器接收到地震波信号并传递给地震仪，然后将地震仪采集到的各个地震波信号整理为三个分量的共炮点道集，在三个分量的共炮点道集中利用拉冬变换提取负速度同相轴，得到每个三分量检波器关于断层界面的反射三分量信号；利用获得的反射三分量信号，对上述各个检波器的反射三分量信号进行时频域主极化方向计算，求取瞬时极化倾角
与方位角φi，利用空间交汇关系确定虚震源点的位置，由于震源位置和虚震源点关于反射界面中心对称，确定虚震源在空间中的坐标后将震源点与虚震源点连接直线段，然后做该直线段的垂直平分线即可准确的找出反射界面位置及倾角关系。本发明直接利用原始反射信号可以避免因为各种预处理方法导致的有效信号丢失，并且利用时频域极化方法可以有效地确定反射信号传播方向，其中的关键是本发明利用时间域几何成像，并不依赖速度模型，而速度模型是超前探测的关键难点，因此本发明无需速度模型的建立，从而能精确确定断层的具体倾角参数以及位置信息。
附图说明
图1是本发明中交汇关系在XOZ面投影的示意图。
具体实施方式
下面将对本发明做进一步说明。
如图所示，以隧道的掘进方向为前方进行专利描述，本发明的具体步骤为：
A、在隧道两帮根据现场施工条件选取其中一帮，在该帮上沿着隧道掘进方向在隧道腰线位置(垂直方向中间)以等间距2m布设16个三分量检波器，由后向前依次编号为R1～R16，在三分量检波器R1的后方10m位置同一高度布置一个炮点，记为S1，在炮点向前激发地震波，16个三分量检波器将接收到的地震波信号传递给地震仪进行记录；
B、根据已知的炮点S1和三分量检波器R1～R16各个位置，建立观测***，并将地震仪采集到的各个地震波信号整理为三个分量的共炮点道集；
C、在三个分量的共炮点道集中利用拉冬变换提取负速度同相轴，得到每个三分量检波器关于断层界面的反射三分量信号；
D、利用步骤C中获得的反射三分量信号，对上述16个三分量检波器的反射三分量信号进行时频域主极化方向计算，求取主极化倾角
与主极化方位角φi，利用空间交汇关系确定断层反射界面倾角与位置，具体确定过程为：
Ⅰ、选取三分量检波器R1关于断层界面的反射三分量信号得出的主极化倾角
与主极化方位角φi，其步骤为：
①反射三分量信号采用解析信号Ai(t)来表示：
式中，
表示采样点优势频率，argAi(t)表示求取Ai(t)的幅角；
②构建协方差矩阵：
式中
表示t时刻的l分量时窗叠加振幅，第l分量在时窗内的均值μl，Tlh(t)表示t时刻
的自适应时窗长度，h为1、2、3；
③通过上述协方差矩阵求取的特征向量归一化后表示为[x1(t2),y1(t2),z1(t2)]；
④根据上述特征向量求取主极化倾角，其表示为：
其中
表示主极化方向与XOZ面之间的夹角，即主极化倾角；
⑤根据上述特征向量求取主极化方位角，可表示为：
其中
表示主极化方向与XOY面之间的夹角，即主极化方位角；
Ⅱ、从而确定该断层界面反射波信号传播的方向
其中(m n p)由主极化倾角
与主极化方位角φi确定，(xR yR zR)表示三分量检波器的空间坐标：
Ⅲ、重复步骤Ⅰ、Ⅱ，确定其他15个三分量检波器关于断层界面的反射三分量信号传播方向
Ⅳ、将各个三分量检波器的采样点优势频率的倒数取平均值，以该平均值的1/4为网格大小，并将三维坐标系网格化，将16个反射信号对应的传播向量
延伸经过各个空间网格，并记录每个网格内经过的射线数目；
Ⅴ、统计各个网格内经过的射线数目，经过射线数目最多的网格即为所有向量的空间交汇点，确定空间交汇点的网格中心为虚震源点，在三维坐标系中做连接炮点与虚震源点的直线段；
Ⅵ、对连接炮点与虚震源点的直线段做垂直平分线，三维坐标系中根据该垂直平分线与隧道掘进中轴线之间夹角确定断层界面的倾角、方位角，该垂直平分线的位置与隧道掘进中轴线之间交点即为断层界面位置。
Claims (1)
1.一种基于三分量反射信号时频域极化参数的断层确定方法，其特征在于，该方法的具体步骤为：
A、在隧道两帮根据现场施工条件选取其中一帮，在该帮上沿着隧道掘进方向在隧道腰线位置以等间距2m布设16个三分量检波器，由后向前依次编号为R1～R16，在三分量检波器R1的后方10m位置同一高度布置一个炮点，记为S1，在炮点向前激发地震波，16个三分量检波器将接收到的地震波信号传递给地震仪进行记录；
B、根据已知的炮点S1和三分量检波器R1～R16各个位置，建立观测***，并将地震仪采集到的各个地震波信号整理为三个分量的共炮点道集；
C、在三个分量的共炮点道集中利用拉冬变换提取负速度同相轴，得到每个三分量检波器关于断层界面的反射三分量信号；
D、利用步骤C中获得的反射三分量信号，对上述16个三分量检波器的反射三分量信号进行时频域主极化方向计算，求取主极化倾角
与主极化方位角φi，利用空间交汇关系确定断层反射界面倾角与位置，具体确定过程为：
Ⅰ、选取三分量检波器R1关于断层界面的反射三分量信号得出的主极化倾角
与主极化方位角φi，其步骤为：
①反射三分量信号采用解析信号Ai(t)来表示：
式中，
表示采样点优势频率，argAi(t)表示求取Ai(t)的幅角；
②构建协方差矩阵：
式中
表示t时刻的l分量时窗叠加振幅，第l分量在时窗内的均值μl，Tlh(t)表示t时刻
的自适应时窗长度，h为1、2、3；
③通过上述协方差矩阵求取的特征向量归一化后表示为[x1(t2),y1(t2),z1(t2)]；
④根据上述特征向量求取主极化倾角，其表示为：
其中
表示主极化方向与XOZ面之间的夹角，即主极化倾角；
⑤根据上述特征向量求取主极化方位角，表示为：
其中
表示主极化方向与XOY面之间的夹角，即主极化方位角；
Ⅱ、确定断层界面反射波信号传播的方向
其中(m n p)由主极化倾角
与主极化方位角φi确定，(xR yR zR)表示三分量检波器的空间坐标：
Ⅲ、重复步骤Ⅰ、Ⅱ，确定其他15个三分量检波器关于断层界面反射波信号传播的方向
Ⅳ、将各个三分量检波器的采样点优势频率的倒数取平均值，以该平均值的1/4为网格大小，并将三维坐标系网格化，将16个断层界面反射波信号传播的方向
延伸经过各个空间网格，并记录每个网格内经过的射线数目；
Ⅴ、统计各个网格内经过的射线数目，经过射线数目最多的网格即为所有向量的空间交汇点，确定空间交汇点的网格中心为虚震源点，在三维坐标系中做连接炮点与虚震源点的直线段；
Ⅵ、对连接炮点与虚震源点的直线段做垂直平分线，三维坐标系中根据该垂直平分线与隧道掘进中轴线之间夹角确定断层界面的倾角、方位角，该垂直平分线的位置与隧道掘进中轴线之间交点即为断层界面位置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910774252.2A CN110531416B (zh) | 2019-08-21 | 2019-08-21 | 一种基于三分量反射信号时频域极化参数的断层确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910774252.2A CN110531416B (zh) | 2019-08-21 | 2019-08-21 | 一种基于三分量反射信号时频域极化参数的断层确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110531416A CN110531416A (zh) | 2019-12-03 |
| CN110531416B true CN110531416B (zh) | 2020-11-20 |
Family
ID=68662491
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910774252.2A Active CN110531416B (zh) | 2019-08-21 | 2019-08-21 | 一种基于三分量反射信号时频域极化参数的断层确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110531416B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114185082B (zh) * | 2021-12-02 | 2023-04-21 | 中国矿业大学 | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5014249A (en) * | 1989-03-17 | 1991-05-07 | Geco A.S. | Method in a real time control system for seismic interference |
| CN103217719A (zh) * | 2013-04-11 | 2013-07-24 | 中国矿业大学 | 基于单一炮检对观测***的煤巷超前探测断失翼煤层的方法 |
| CN104570111A (zh) * | 2015-01-21 | 2015-04-29 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
| CN104678437A (zh) * | 2015-03-18 | 2015-06-03 | 河南理工大学 | 一种矿井槽波信号波场分离方法 |
| CN105652322A (zh) * | 2016-01-07 | 2016-06-08 | 中国科学院地球化学研究所 | 多分量地震数据的t-f-k域极化滤波方法 |
| CN105974472A (zh) * | 2016-05-13 | 2016-09-28 | 中国矿业大学 | 一种基于反射信号的巷道超前探测速度建模方法 |
| CN109164483A (zh) * | 2018-08-29 | 2019-01-08 | 中国科学院地球化学研究所 | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9176243B2 (en) * | 2006-02-24 | 2015-11-03 | Hanner Zueroher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
-
2019
- 2019-08-21 CN CN201910774252.2A patent/CN110531416B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5014249A (en) * | 1989-03-17 | 1991-05-07 | Geco A.S. | Method in a real time control system for seismic interference |
| CN103217719A (zh) * | 2013-04-11 | 2013-07-24 | 中国矿业大学 | 基于单一炮检对观测***的煤巷超前探测断失翼煤层的方法 |
| CN104570111A (zh) * | 2015-01-21 | 2015-04-29 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
| CN104678437A (zh) * | 2015-03-18 | 2015-06-03 | 河南理工大学 | 一种矿井槽波信号波场分离方法 |
| CN105652322A (zh) * | 2016-01-07 | 2016-06-08 | 中国科学院地球化学研究所 | 多分量地震数据的t-f-k域极化滤波方法 |
| CN105974472A (zh) * | 2016-05-13 | 2016-09-28 | 中国矿业大学 | 一种基于反射信号的巷道超前探测速度建模方法 |
| CN109164483A (zh) * | 2018-08-29 | 2019-01-08 | 中国科学院地球化学研究所 | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110531416A (zh) | 2019-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
| CN102395902B (zh) | 使用快速面向目标照明计算的地震成像***及方法 | |
| Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
| CN112578428B (zh) | 一种基于巷道垂向虚拟测线的散射多波超前探测方法 | |
| CN105676279A (zh) | 一种采用同心圆等炮检距的地震反射数据采集方法 | |
| CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
| CN109557582A (zh) | 一种二维多分量地震资料偏移成像方法及*** | |
| CN106199704B (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
| CN104166159B (zh) | 四维微地震监测的裂缝形态处理方法和*** | |
| CN104991268A (zh) | 一种真振幅偏移成像方法 | |
| CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
| CN110632667B (zh) | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 | |
| CN111103620A (zh) | 一种岩巷超前探测三维偏移成像方法 | |
| CN110531416B (zh) | 一种基于三分量反射信号时频域极化参数的断层确定方法 | |
| CN208334667U (zh) | 一种用于软岩隧道长距离三维超前地质预报装置 | |
| Fujimoto et al. | DAS 3DVSP survey at Stratigraphic Test Well (Hydrate-01) | |
| CN110531413B (zh) | 一种小断层超前三维可视化建模方法 | |
| CN110531418B (zh) | 一种基于希尔伯特极化成像的断点三维精细定位方法 | |
| CN109143362A (zh) | 基于共散射角道集的散射波分离方法 | |
| CN103492909A (zh) | 动态地震勘探方法 | |
| CN111158050A (zh) | 数据采集***、方法及隧道地震波超前预报方法 | |
| CN110579799A (zh) | 一种等旅行时间间隔的地震采集观测方法及*** | |
| Ren et al. | Pre-stack elastic reverse time migration in tunnels based on cylindrical coordinates | |
| US3354985A (en) | Seismic prospecting method of locating a suspected fault | |
| CN109375261B (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 |