CN113960532A - 一种基于假想源的二次定位计算的微地震定位方法 - Google Patents
一种基于假想源的二次定位计算的微地震定位方法 Download PDFInfo
- Publication number
- CN113960532A CN113960532A CN202111221910.9A CN202111221910A CN113960532A CN 113960532 A CN113960532 A CN 113960532A CN 202111221910 A CN202111221910 A CN 202111221910A CN 113960532 A CN113960532 A CN 113960532A
- Authority
- CN
- China
- Prior art keywords
- source
- travel time
- data
- seismic
- positioning
- 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
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000004364 calculation method Methods 0.000 title claims abstract description 27
- 239000003086 colorant Substances 0.000 claims description 2
- 238000012360 testing method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000004807 localization Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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/288—Event detection in seismic signals, e.g. microseismics
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于假想源的二次定位计算的微地震定位方法,其步骤包括首先设定一系列的“假想源”,读取检波器数据,模型数据,微地震事件走时数据,根据地震的互易性原理,计算各个检波器作为源到“假想源”的走时,把上述走时数据作为对比数据集,以实际观察走时数据与数据集进行对比匹配,计算每一个地震事件满足条件的射线数和残差和,迅速确定微地源事件发生的微区域和参考“假想源”,然后再在微区域内以参考“假想源”为基点,采用最小二乘法二次精准定位微地震事件位置。本发明的方法针对一个区域的震源定位,只需要计算一次时间场,不需要重复性计算,避免了反复射线追踪过程的时间损耗,节约大量时间,同时可扩展到三维情况应用计算。
Description
技术领域
本发明属于微地震震源定位技术领域,具体涉及一种基于假想源的二次定位计算的微地震定位方法。
背景技术
地震震源定位可以划分为两个时期,以1912年Geiger的开拓性工作为界线,之前是几何作图法,包括单台资料法和多台资料法,单台资料法包括方位角法和出射角法,多台资料法包括虚波速度法、虚拟法、交切法、引中线法、和达法和高桥法等方法,几何作图法都是通过震中所在轨迹线的相交,确定震中位置,其定位精度低。1912年Geiger开拓性工作之后为数字定位法,1975年Lee等利用Geiger理论,实现了震源定位程序,连续给出了HYPO71-86系列程序。考虑到台站位置以及速度结构对定位的影响,发展了联合定位方法,震源位置和台站校正的联全定位方法(JED,JHD),震源位置与速度结构联合反演方法(SSH),由JED方法还发展出来了主事件定位方法(ATD)。2000年,由Waldhauser和Ellsworth提出了双差定位法,与主事件相对定位法不同,事件对的距离不受限制,大大提高了方法的适用性。在反演方法上逐渐引入非线性定位方法,如牛顿法、全局搜索方法,Powell法、遗传算法、模拟退火算法等。几何作图法具有稳健、效率高的优势,但定位精度低;线性定位法无需速度模型精度,但对初至拾取的精度有较高要求;非线性定位法对初至拾取较为敏感,对速度模型的精度要求较高,但计算量较小。
发明内容
本发明提供一种基于假想源的二次定位计算的微地震定位方法,能够精准定位微地震事件源。
本发明的技术方案是:一种基于假想源的二次定位计算的微地震定位方法,首先设定一系列的“假想源”,读取检波器数据,模型数据,微地震事件走时数据,根据地震的互易性原理,计算各个检波器作为源到“假想源”的走时,把上述走时数据作为对比数据集,以实际观察走时数据与数据集进行对比匹配,计算每一个地震事件满足条件的射线数和残差和,迅速确定微地源事件发生的微区域和参考“假想源”,然后再在微区域内以参考“假想源”为基点,采用最小二乘法二次精准定位微地震事件位置,其步骤具体包括:
步骤1,利用地震的互易性原理,计算时间场值,对观察到的走时数据进行匹配,迅速定位地震发生的小分块和“假想源”,初次匹配确定地震发生域;
1.1、建立分块-节点模型
首先读取检波器数据、模型数据、地震信息,对地质模型采用分块-节点法进行数字化,利用不同灰度块或不同颜色代表不同速度的介质,然后在不同速度块利用节点对模型进行细分,并且节点间距不固定,是可变的,利用复杂模型的数字化读取模型信息;采用井中、井地结合、井间多种方法布设检波器进行微地震信号采集,读取地震信息,即地震事件实际走时;在建立的模型内读取检波器的水平坐标和深度坐标数据;
1.2、自动加密模型获取小分块中点
以1.1中建立的模型的原有纵横坐标节点为依据,在纵横坐标相邻节点的中点设置多个加密点进行加密,各加密点形成的矩形小分块,将小分块中点进行标记作为“假想源”;
1.3、利用互易性计算时间场
利用互易性把检波器作为“震源”,利用射线追踪方法计算从检波器传出来的波到“假想源”的走时场,每个检波器形成1个走时场数据集;
1.4、匹配对比确定微震发生域
将每个检波器中观察的地震事件的实际走时与1.3中获得的走时场数据集进行对比,数据相同时计射线数为1,统计射线数,统计射线数时应该满足以下条件:
其中t实为检波实际观察到的走时数据,t假为假想源的走时,即1.3所得到假想源时间场值,Dx为模型中两个圆形的原始网格节点水平方向的间距,Dh为模型中两个圆形的原始网格节点深度方向的间距,V分为以小分块中点为中心,以圆点为角点的小分块内波传播速度;
当实际震源离假想源较远,在小分块边源时,有可能统计的射线数在两个小分块相同,均为统计的射线数的最大值,此时不利用第二步的二次定位,为了避免此时无法初始匹域震源所在小分块,计匹配射线时数,同时统计走时残差和,其计算公式如下:
其中t和代表残差和,t实和t假分别代表微地震事件实际观察走时和1.3中理论计算假想源的走时,i代表检波器序号,n代表总检波器数量;
射线数最大,残差和最小位置即为初步匹配的假想源位置,该假想源所在的小分块即为微震事件发生的可能区域;
步骤2,利用假想源进行二次定位,准确确定震源位置;
2.1、二次参考事件精准确定震源位置
二次定位是相对定位方法,根据如下计算公式:
2.2、计算上述公式(3)中的导数值
确定实际走时的时间场对坐标的导数值,利用模型中的原始网格节点,加密点,小分块中点的走时数据,建立时间场与坐标系的二次趋势面方程:
ti=a0+a1x+a2h+a3x2+a4xh+a5h2 (4)
其中ti为实际走时,即从检波器到节点处的走时,i代表节点的序号,x代表水平坐标,h代表深度坐标,a0、a1、a2、a3、a4、a5分别是系数;
将数据代入(4)式采用最小二乘法计算获得系数,同时可得对x,对h的导数值:
最终代入公式(3)二次精准定位震源位置。
本发明的优点是:
1)本发明利用地震互易性,一次计算完所有“假想源”的走时场,避免了反复射线追踪过程的时间损耗,能够精准定位微地震事件源;
2)本发明利用射线数和走时残差和,快速匹配锁定震源区,避免了扫描法或非线性方法反复调整源区的过程,节省了大量计算损耗;
3)本发明利用参考“假想源”在小分块中进行二次定位,是相对定位方法,避免了模型速度依赖,且计算量很小;
4)本发明在小分块中二次定位计算时,采用趋势面拟合方法,很好的解决了小范围内时间场导数的求取问题;
5)本发明不仅可适用的二维情况,同时可扩展到三维情况,将由二维的矩形域扩展成三维情况的矩形体。
附图说明
图1是本发明流程示意图;
图2是本发明地质模型采用分块-节点法进行数字化后的示意图;
图3是本发明读取地震信息时的示意图,a为单井采集模式,b为井地结合采集模式,c为双井采集模式;
图4是本发明自动加密后的模型示意图;
图5是本发明利用互易性计算各检波器为源的时间场的数据集;
图6是本发明二次定位原理示意图;
图7是本发明实施例1各向同性均匀介质中定位测试中数据模型化结果和检波器位置分布示意图;
图8是本发明实施例1各向同性均匀介质中定位测试中走时场图;
图9是本发明实施例1各向同性均匀介质中定位测试中两个地震事件射线数与小分块残差和结果图;
图10是本发明实施例1各向同性均匀介质中定位测试中定位结果图;
图11是本发明实施例2层状介质中的测试中模型化结果及检波器分布示意图;
图12是本发明实施例2层状介质中的测试中走时场图;
图13是本发明实施例2层状介质中的测试中两个地震事件射线数与小分块残差和结果图;
图14是本发明实施例2层状介质中的测试中定位结果图;
具体实施方式
下面结合附图对本发明做清楚完整的描述,以使本领域的技术人员在不需要作出创造性劳动的条件下,能够充分实施本发明。
如图1所示,一种基于假想源的二次定位计算的微地震定位方法,利用“假想源”分两步进行精准定位震源位置,其步骤包括:
步骤1,利用地震的互易性原理,计算时间场值,对观察到的走时数据进行匹配,迅速定位地震发生的小分块和“假想源”,初次匹配确定地震发生域;
步骤2,利用假想源进行二次定位,准确确定震源位置。
上述步骤1具体包括:
1.1、建立分块-节点模型
首先读取检波器数据、模型数据、地震信息(即微地震事件走时),如表1-3所示,对地质模型采用分块-节点法进行数字化,如图2所示,不同灰度块代表不同速度的介质,然后在不同速度块利用节点(图2中X符号),对模型进行细分,并且节点间距不固定,是可变的,利用复杂模型的数字化读取模型信息,如表2所示;
在建立的模型内读取检波器的水平坐标和深度坐标数据,如表1所示;
读取地震信息,采用井中、井地结合、井间多种方法布设检波器进行微地震信号采集,读取地震信息,如图3所示,信息数据如表3所示;
表1检波器数据
<u>检波器编号</u> | <u>水平坐标</u> | <u>深度坐标</u> |
<u>1</u> | <u>X1</u> | <u>H1</u> |
<u>2</u> | <u>X2</u> | <u>H2</u> |
<u>…</u> | <u>…</u> | <u>…</u> |
表2模型数据格式
<u>节点编号</u> | <u>水平坐标</u> | <u>深度坐标</u> | <u>速度值</u> |
<u>1</u> | <u>X<sub>1</sub></u> | <u>H<sub>1</sub></u> | <u>V<sub>1</sub></u> |
<u>2</u> | <u>X<sub>1</sub></u> | <u>H<sub>2</sub></u> | <u>V<sub>2</sub></u> |
<u>…</u> | <u>…</u> | <u>…</u> | <u>…</u> |
<u>N</u> | <u>X<sub>1</sub></u> | <u>H<sub>N</sub></u> | <u>V<sub>N</sub></u> |
<u>N+1</u> | <u>X<sub>2</sub></u> | <u>H<sub>1</sub></u> | <u>V<sub>N+1</sub></u> |
<u>N+2</u> | <u>X<sub>2</sub></u> | <u>H<sub>2</sub></u> | <u>V<sub>N+2</sub></u> |
<u>…</u> | <u>…</u> | <u>…</u> | <u>…</u> |
<u>N+N</u> | <u>X<sub>2</sub></u> | <u>H<sub>N</sub></u> | <u>V<sub>N+N</sub></u> |
<u>…</u> | <u>…</u> | <u>…</u> | <u>…</u> |
<u>(M-1)*N+1</u> | <u>X<sub>m</sub></u> | <u>H<sub>1</sub></u> | <u>V<sub>(M-1)*N+1</sub></u> |
<u>(M-1)*N+2</u> | <u>X<sub>m</sub></u> | <u>H<sub>2</sub></u> | <u>V<sub>(M-1)*N+2</sub></u> |
<u>…</u> | <u>…</u> | <u>…</u> | <u>…</u> |
<u>M*N</u> | <u>X<sub>m</sub></u> | <u>H<sub>N</sub></u> | <u>V<sub>M*N</sub></u> |
表3地震信息数据格式
1.2、自动加密模型获取小分块中点
以图2中原有纵横坐标节点为依据,在纵横坐标相邻节点的中点设置多个加密点,自动加密结果如图4,图4中五星代表检波器位置,圆点是图2中原始网格节点位置,三角形是加密之后的小分块中点,即假想源位置,X位置是加密点。
1.3、利用互易性计算时间场
从震源发出来的波到检波器的走时与从检波器发出的波到震源的走时是相同的,正是利用此地震互易性原理,在不知道真实震源位置的情况下,把检波器作为“震源”,利用射线追踪方法计算从检波器传出来的波到“假想源”,即小分块中点的走时场,如图5所示是图4中六个检波器作为源的走时场,每个检波器形成1个走时场数据集。
1.4、匹配对比确定微震发生域
假设某一个“假想源”是真实的震源位置,那震源到6个检波器的6条射线走时,必然分别在图5中的6个数据集中,第一个检波器中观察的实际走时,与数据集1进行对比,震源位置的走时数据与实际走时相同,此时计射线数为1,同样对比第二、三、四、五、六检波器的实际走时与走时数据集,统计其射线数,实际震源位置射线数必然为6,把震源限定在图4中以三角形为中心,以圆点为角点的小分块内,此时统计射线数时应该满足以下条件:
其中t实为检波实际观察到的走时数据,t假为假想源的走时,即1.3计算所得到假想源时间场值,Dx为图4中两个原始网格节点水平方向的间距,Dh为图4中两个原始网格节点深度方向的间距,V分为以三角形为中心,以圆点为角点的小分块内波传播速度。
当实际震源离假想源较远,在小分块边源时,有可能统计的射线数在两个小分块相同,均为统计的射线数的最大值,此时不利于第二步的二次定位,为了避免此时无法初始匹配震源所在小分块,统计匹配射线数时,同时统计走时残差和,其计算公式如下:
其中t和代表残差和,t实和t假分别代表微地震事件实际观察走时和理论计算假想源的走时,i代表检波器序号,n代表总检波器数量。
通过公式(1)和公式(2)计算的射线数和残差和如表4和表5,易知射线数最大,残差和最小位置即为初步匹配的假想源位置,该假想源所在的小分块即为微震事件发生的可能区域。
表4各分块符合条件的射线数
表5各小分块射线残差和
2.1、二次参考事件精准确定震源位置
二次定位是相对定位方法,如图6所示,三角形是小分块中点,为假想地震源,实际震源位置在圆点网格点所围小分块内一点,根据如下计算公式:
2.2、计算上述公式(3)中的导数值
确定实际走时的时间场对坐标的导数值,利用图6中原始网格节点,加密点,小分块中点共9个点的走时数据,建立时间场与坐标系的二次趋势面方程:
ti=a0+a1x+a2h+a3x2+a4xh+a5h2 (4)
其中ti为实际走时,即从检波器到节点处的走时,i代表节点的序号,i最大为9,x代表水平坐标,h代表深度坐标,因此把9个点的走时数据全部代入方程(4)变为9个方程的方程组,此时需要确定(4)式的系数,此方程的求解同样采用最小二乘法,获得系数后,可得对x,对h的导数值:
最终代入公式(3)二次精准定位震源位置。
实施例1:各向同性均匀介质中定位测试
数据模型化结果和检波器位置如图7,设置2个地震事件6个检波器的接收的走时数据,如表6所示。
表6均匀介质中走时数据
根据预设的到时数据与波场走时正演结果值,可以确定两个震源范围位于图7中的矩形区域框内。
读取检波器、模型数据、地震信息后,根据互易性,计算检波器作为震源的走时图,如图8所示,按图4中的方式进行网格加密,匹配表6的走时与图8中小分块中心的走时值,进行满足条件射线数与小分块残差和计算,两个地震结果如图9所示。
依据图9中的射线数最大值和残差和最小值,确定的小分块中心如图10所示,然后进行二次参考事件定位,其定位结果如图10所示。
实施例2:层状介质中的测试,选择了井地联合的检波器布设方式,其模型化结果及检波器分布图11,地震事件数据如表7所示,其震源位置应落在图11的矩形框内。
表7层状介质微地震走时数据
读取检波器、模型数据、地震信息后,根据互易性,计算检波器作为震源的走时图如图12所示。
按图4中的方式进行网格加密,匹配表7的走时与图12中小分块中心的走时值,进行满足条件射线数与小分块残差和计算,两个地震结果如图13所示,a为射线数,b为残差和。
依据图13中的射线数最大值和残差和最小值,确定的小分块中心如图14所示,然后进行二次参考事件定位,其定位结果如图14所示。
以上对本发明的较佳实施例进行了描述,需要指出的是,本发明,并不局限于上述特定实施方式,其中未尽详细描述的设备和结构应该理解为用本领域中的普通方式予以实施;任何熟悉本领域的技术人员,在不脱离本发明技术方案范围情况下,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。
Claims (4)
1.一种基于假想源的二次定位计算的微地震定位方法,其特征在于,利用“假想源”分两步进行精准定位震源位置,其步骤包括:
步骤1,利用地震的互易性原理,计算时间场值,对观察到的走时数据进行匹配,迅速定位地震发生的小分块和“假想源”,初次匹配确定地震发生域;
1.1、建立分块-节点模型
首先读取检波器数据、模型数据、地震信息,对地质模型采用分块-节点法进行数字化,利用不同灰度块或不同颜色代表不同速度的介质,然后在不同速度块利用节点对模型进行细分,并且节点间距不固定,是可变的,利用复杂模型的数字化读取模型信息;采用井中、井地结合、井间多种方法布设检波器进行微地震信号采集,读取地震信息,即地震事件实际走时;在建立的模型内读取检波器的水平坐标和深度坐标数据;
1.2、自动加密模型获取小分块中点
以1.1中建立的模型的原有纵横坐标节点为依据,在纵横坐标相邻节点的中点设置多个加密点进行加密,各加密点形成的矩形小分块,将小分块中点进行标记作为“假想源”;
1.3、利用互易性计算时间场
利用互易性把检波器作为“震源”,利用射线追踪方法计算从检波器传出来的波到“假想源”的走时场,每个检波器形成1个走时场数据集;
1.4、匹配对比确定微震发生域
将每个检波器中观察的地震事件的实际走时与1.3中获得的走时场数据集进行对比,数据相同时计射线数为1,统计射线数,直至所有地震事件统计计算完成,统计射线数时应该满足以下条件:
其中t实为检波实际观察到的走时数据,t假为假想源的走时,即1.3所得到假想源时间场值,Dx为模型中两个圆形的原始网格节点水平方向的间距,Dh为模型中两个圆形的原始网格节点深度方向的间距,V分为以小分块中点为中心,以圆点为角点的小分块内波传播速度;
当实际震源离假想源较远,在小分块边源时,有可能统计的射线数在两个小分块相同,均为统计的射线数的最大值,此时在计匹配射线时数,同时统计走时残差和,其计算公式如下:
其中t和代表残差和,t实和t假分别代表微地震事件实际观察走时和1.3中理论计算假想源的走时,i代表检波器序号,n代表总检波器数量;
射线数最大,残差和最小位置即为初步匹配的假想源位置,该假想源所在的小分块即为微震事件发生的可能区域;
步骤2,利用假想源进行二次定位,准确确定震源位置;
2.1、二次参考事件精准确定震源位置
二次定位是相对定位方法,根据如下计算公式:
2.2、计算上述公式(3)中的导数值
确定实际走时的时间场对坐标的导数值,利用模型中的原始网格节点,加密点,小分块中点的走时数据,建立时间场与坐标系的二次趋势面方程:
ti=a0+a1x+a2h+a3x2+a4xh+a5h2 (4)
其中ti为实际走时,即从检波器到节点处的走时,i代表节点的序号,x代表水平坐标,h代表深度坐标,a0、a1、a2、a3、a4、a5分别是系数;
可得对x,对h的导数值:
最终代入公式(3)二次精准定位震源位置。
2.根据权利要求1所述的一种基于假想源的二次定位计算的微地震定位方法,其特征在于,上述步骤2.2中系数的确定方式是将数据代入(4)式采用最小二乘法计算获得系数。
3.根据权利要求1所述的一种基于假想源的二次定位计算的微地震定位方法,其特征在于,上述步骤2.2中i最大为9。
4.根据权利要求1所述的一种基于假想源的二次定位计算的微地震定位方法,其特征在于,上述步骤1.4中检波器的数量为6个,统计射线数为6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111221910.9A CN113960532B (zh) | 2021-10-20 | 2021-10-20 | 一种基于假想源的二次定位计算的微地震定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111221910.9A CN113960532B (zh) | 2021-10-20 | 2021-10-20 | 一种基于假想源的二次定位计算的微地震定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113960532A true CN113960532A (zh) | 2022-01-21 |
CN113960532B CN113960532B (zh) | 2024-05-24 |
Family
ID=79464880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111221910.9A Active CN113960532B (zh) | 2021-10-20 | 2021-10-20 | 一种基于假想源的二次定位计算的微地震定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113960532B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114721049A (zh) * | 2022-03-11 | 2022-07-08 | 新疆维吾尔自治区地震局 | 一种虚拟地震台方位角参加偏网地震定位方法 |
CN114879251A (zh) * | 2022-06-23 | 2022-08-09 | 中国地震局地球物理研究所 | 一种基于射线追踪技术的地震定位图形方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0330097D0 (en) * | 2003-12-29 | 2004-02-04 | Westerngeco Ltd | Method for monitoring hydraulic stimulation |
CN103823243A (zh) * | 2012-11-15 | 2014-05-28 | 地球物理集团公司 | 用于分离在连续数据采集地震勘测期间记录的数据的方法 |
US20170045634A1 (en) * | 2015-02-10 | 2017-02-16 | Michael D. Collins | Method of using a parabolic equation model for range-dependent seismo-acoustic problems |
CN110133715A (zh) * | 2019-05-29 | 2019-08-16 | 长江大学 | 一种基于初至时差和波形叠加的微地震震源定位方法 |
CN110261903A (zh) * | 2019-06-10 | 2019-09-20 | 中北大学 | 一种基于逆时能量聚焦的地下震源被动定位方法 |
-
2021
- 2021-10-20 CN CN202111221910.9A patent/CN113960532B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0330097D0 (en) * | 2003-12-29 | 2004-02-04 | Westerngeco Ltd | Method for monitoring hydraulic stimulation |
CN103823243A (zh) * | 2012-11-15 | 2014-05-28 | 地球物理集团公司 | 用于分离在连续数据采集地震勘测期间记录的数据的方法 |
US20170045634A1 (en) * | 2015-02-10 | 2017-02-16 | Michael D. Collins | Method of using a parabolic equation model for range-dependent seismo-acoustic problems |
CN110133715A (zh) * | 2019-05-29 | 2019-08-16 | 长江大学 | 一种基于初至时差和波形叠加的微地震震源定位方法 |
CN110261903A (zh) * | 2019-06-10 | 2019-09-20 | 中北大学 | 一种基于逆时能量聚焦的地下震源被动定位方法 |
Non-Patent Citations (1)
Title |
---|
王辉;梁苗;朱梦博;: "基于单纯形-最短路径射线追踪的微震震源混合定位算法", 中国矿业, no. 10, 15 October 2020 (2020-10-15) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114721049A (zh) * | 2022-03-11 | 2022-07-08 | 新疆维吾尔自治区地震局 | 一种虚拟地震台方位角参加偏网地震定位方法 |
CN114721049B (zh) * | 2022-03-11 | 2024-04-09 | 新疆维吾尔自治区地震局 | 一种虚拟地震台方位角参加偏网地震定位方法 |
CN114879251A (zh) * | 2022-06-23 | 2022-08-09 | 中国地震局地球物理研究所 | 一种基于射线追踪技术的地震定位图形方法及系统 |
CN114879251B (zh) * | 2022-06-23 | 2023-02-21 | 中国地震局地球物理研究所 | 一种基于射线追踪技术的地震定位图形方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113960532B (zh) | 2024-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104656141B (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN107167786A (zh) | 卫星激光测高数据辅助提取高程控制点方法 | |
Peters et al. | Application of prediction analysis to hypocenter determination using a local array | |
CN100549731C (zh) | 检波点二次定位方法 | |
CN113960532B (zh) | 一种基于假想源的二次定位计算的微地震定位方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN106094029A (zh) | 利用偏移距矢量片地震数据预测储层裂缝的方法 | |
WO2009034567A2 (en) | Complex analysis of kinematics for non-hyperbolic moveout corrections | |
CN105093299A (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
CN105759311A (zh) | 一种近实时地震震源位置定位方法 | |
CN106249295A (zh) | 一种井中微地震p、s波联合快速定位方法及系统 | |
CN104181593B (zh) | 一种三维无射线追踪回折波层析成像方法及装置 | |
CN113850908A (zh) | 一种考虑路径延长因素的地闪回击定位数据的优化方法 | |
CN110146924B (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
CN105573963A (zh) | 一种电离层水平不均匀结构重构方法 | |
CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
CN109085642B (zh) | 一种各向异性介质微地震事件定位方法 | |
CN106443776B (zh) | 一种基于时间切片法的海底地震仪重定位方法 | |
CN113740915A (zh) | 一种球坐标系重力和接收函数联合反演地壳结构参数的方法 | |
Liu et al. | A new method for OBS relocation using direct water-wave arrival times from a shooting line and accurate bathymetric data | |
CN102478663B (zh) | 一种三维地震观测系统偏移噪声获取方法及装置 | |
CN109212594B (zh) | 一种各向异性介质纵横波联合定位方法 | |
CN108508479B (zh) | 一种空地井立体重磁数据协同目标位置反演方法 | |
CN104155703A (zh) | 评价三维观测系统的方法和设备 | |
CN105572741B (zh) | 一种计算3d高频静校正量的方法 |
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 |