CN110389377A - 基于波形互相关系数相乘的微震偏移成像定位方法 - Google Patents
基于波形互相关系数相乘的微震偏移成像定位方法 Download PDFInfo
- Publication number
- CN110389377A CN110389377A CN201810366402.1A CN201810366402A CN110389377A CN 110389377 A CN110389377 A CN 110389377A CN 201810366402 A CN201810366402 A CN 201810366402A CN 110389377 A CN110389377 A CN 110389377A
- Authority
- CN
- China
- Prior art keywords
- imaging
- point
- microseism
- correlation coefficient
- cross
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000004807 localization Effects 0.000 title claims abstract description 30
- 238000013508 migration Methods 0.000 title claims abstract description 17
- 230000005012 migration Effects 0.000 title claims abstract description 17
- 238000012544 monitoring process Methods 0.000 claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 8
- FESBVLZDDCQLFY-UHFFFAOYSA-N sete Chemical compound [Te]=[Se] FESBVLZDDCQLFY-UHFFFAOYSA-N 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 6
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000006835 compression Effects 0.000 abstract description 3
- 238000007906 compression Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 230000001788 irregular Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 206010019133 Hangover Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- 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
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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
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
本发明涉及一种基于波形互相关系数相乘的微震偏移成像定位方法，属于油气储层压裂微震监测技术领域。其包括如下步骤：实际微震数据预处理；地下成像点划分，计算所有成像点到地面各个接收点的旅行走时；针对某一个成像点，消除该成像点到不同接收点的到时差，计算时差校正后的微地震记录的互相关系数；计算该成像点的互相关能量值E；重复第三步和第四步获得地下所有成像点的能量值E；设置能量阈值E 0 ，得到大于阈值的能量极大值成像点位置，作为微震源的发生位置；本发明通过使用波形互相关系数相乘的形式，实现了对噪音的双重压制，极大的提高了抗噪能力；避开了初至极性对定位成像结果的干扰，具有更高的计算效率，满足微震的实时监测。
Description
技术领域
本发明涉及一种基于波形互相关系数相乘的微震偏移成像定位方法，属于油气储层压裂微震监测技术领域。
背景技术
在地面微震资料处理中，关键的步骤是要准确的识别出微震有效事件，并且要对识别出的微震有效事件进行精确定位，从而有效地指导水力压裂工作和裂隙预测工作。然而，由于地下介质的吸收作用和地面噪声的干扰，相比于井中监测相比，地面微震记录信噪比较低，致使震相难以识别和初至时间不易拾取，则限制了一些基于震相识别和到时拾取的传统定位方法的应用，使不需震相识别和到时拾取的微震偏移成像成为地面微震监测震源定位的研究热点。
基于Kirchoff偏移成像方法包括震源扫描（source-scanning algorithm，SSA）和地震发射层析成像( seismic-emission-tomography，SET)。SSA定位方法的主要优势是不需要拾取微震到时，也不用计算理论地震图，只需通过系统扫描一定位置范围和时段，就可以获得整个震源序列分布。但是其成像的亮度函数是对波形取绝对值后叠加，在增强有效信号的同时也增强了噪声，对信噪比较低的地震资料处理效果不佳。 SET是地面微地震监测常用的一种震源定位算法，对于低信噪比的地震记录，该方法依然能够较准确的对微地震事件进行定位，但是定位精度易受到初至波极性的影响。针对SET定位方法的缺陷，通过加入震源机制极性校正项消除了初至极性的影响，并在实际应用中取得了较好的效果，但是震源机制反演参数较多，使得计算效率较低。基于互相关法的震源定位方法是利用不同地震道所接收到同震源的地震事件通常在波形上具有相似性，而随机噪音是互不相关的特点，有效地压制了噪音影响，但是同样受到初至波极性的影响，虽然可以对波形互相关系数取绝对值再叠加，既可以消除初至极性影响，又可以不受台站布设方式的影响，但是取绝对值的同时也增强噪音的影响，定位成像效果仍然不理想。因此，有必要寻找一种既不受地面微震资料信噪比低、不规则台站的布设方式和初至极性反转等因素影响，又能提高计算效率，满足地面微震信号的实时监测要求的微震定位方法。
发明内容
针对现有技术存在的上述缺陷，本发明提出了基于波形互相关系数相乘的微震偏移成像定位方法。该方法将传统的叠加成像思想改为相乘成像，将现有的波形互相关能量叠加的成像条件改进为波形互相关能量相乘的成像条件，实现了对噪音的双重压制，较大地提高了压制噪音的能力；并且不需对初至极性进行校正，适合微地震的实时监测；也不受制于台站布设方式的影响，较好地满足了地形复杂地区的微地震监测工作。
本发明所述的基于波形互相关系数相乘的微震偏移成像定位方法，包括如下步骤：
（1）对实际微震数据进行预处理，得到有效地震道的微震数据；
（2）将监测的地下区域在x, y, z三个坐标方向进行三维网格划分，三个方向上的划分尺度不一定相同，将每个网格点设置为成像点，用基于地震走时扰动插值的射线追踪方法计算每个成像点到达地面各个接收点的旅行时间；
（3）利用第二步中得到某一个成像点η到各个接收点i的旅行走时对相应的微地震记录进行走时校正，消除该成像点到不同接收点的时差，然后对时差校正后的微地震记录进行两两互相关，计算对应的互相关系数；
利用第二步中某一个成像点η的旅行走时作为微震信号到达各个地面接收点i的旅行走时t ηi ，并设微震信号的发震时刻T=t 0 ，则到达时刻为t=t ηi + t 0 ，确定时间窗口宽度W，则任意两个地震道的波形互相关系数为：
其中，W为所取时间窗口内的采样点个数，u i 为第i个地面接收点的微震记录幅值；u j 为第j个地面接收点微震记录的幅值；t 0 为发震时刻；t ηi 为地下成像点η到第i个地面接收点的旅行走时；t ηj 为地下成像点η到第j个接收点的旅行走时；t k 为所取时间窗口内的第k个采样点时间；
（4）将步骤（3）中得到的波形互相关系数进行相乘计算该成像点的互相关能量值E，计算公式为：
其中，N为地震道数，E（η，t 0 ）表示成像点η在发震时刻T=t 0 的成像能量值；
（5）针对每一个成像点重复步骤（3）和步骤（4），获得地下所有成像点的互相关能量值E；
（6）设置一个合适的能量阈值E 0 ，将步骤（4）中所得到的所有成像点的能量值E与能量阈值E 0 进行比较，通过搜寻出大于所设能量阈值的能量极大值成像点（E（η，t 0 ）> E 0 ）；即所有符合要求的成像点η为微震事件发生位置，所设的发震时刻t 0 即为微震事件的发震时刻；
本发明的有益效果是：采用本发明所述的基于波形互相关系数相乘的微震偏移成像定位方法，（1）提高了抗噪能力。针对地面微震资料信噪比低，震相难以识别和初至不易拾取的问题，利用微震有效事件波形相似，而随机噪音不相干的特点，通过计算波形互相关函数压制噪音的影响，并且将计算得到的波形互相关系数进行相乘计算，更是对噪音的第二次压制，实现了对随机噪音的双重压制，极大地提高了抗噪能力。（2）避免了震源机制反演和微震极性判定，提高了计算效率。由于对波形互相关系数进行了相乘运算，不存在常用的叠加技术使正负极性微震波形相加抵消的问题，不需要通过反演震源机制修正初至负极性，就可以消除初至极性反转对定位精度的影响，减少了震源机制反演，提高了微震事件定位的计算效率。（3）适合于地面微震监测台站的不规则布设。当地面微震监测台站为不规则布设方式时，监测台站数目较少，常规基于叠加能量成像方法，由于台站数目较少，震源空间成像分辨率较差，利用波形互相关系数相乘的思想将震源处的能量快速“聚焦”，使震源空间成像分辨率明显提高，能够快速和准确的确定出震源位置。
附图说明
图1是本发明的流程框图。
图2是本发明的地面微震监测系统和速度模型，其中实心倒三角为接收点，实心小圆点为在目标成像区域内网格划分的成像点，五角星为真实震源位置。
图3是合成的微震记录图。
图4是加入一定程度的随机噪音后的低信噪比的微震记录图。
图5是基于振幅叠加的微震定位方法在空间的成像结果图。
图6是基于波形互相关系数叠加的微震定位方法在空间的成像结果图。
图7是基于震源机制和S参数联合的微震定位方法在空间成像结果图。
图8是基于波形互相关系数相乘的偏移成像的微震定位方法的空间成像结果图。
图9是西南地区某一油田的地面微震监测台站位置示意图。
图10是实际射孔的微震记录图。
图11是常规振幅叠加的微震定位方法在Z方向上的水平面成像结果图。
图12（a）是常规叠加的微震定位方法在XZ面的成像结果图。
图12（b）是常规叠加的微震定位方法在YZ面的成像结果图。
图13是基于波形互相关系数相乘的微震偏移成像定位方法在Z方向上的水平面成像结果图。
图14（a）是基于波形互相关系数相乘的微震偏移成像定位方法在XZ面的成像结果图。
图14（b）是基于波形互相关系数相乘的微震偏移成像定位方法在YZ面的成像结果图。
具体实施方式
为了使本发明的特点、优势更加明显，下面结合说明书附图和具体实施例对本发明作进一步详细描述。
步骤一：对实际微震数据进行预处理，得到有效地震道的微震数据。
步骤二：将监测的地下区域在x, y, z三个坐标方向进行三维网格划分，三个方向上的划分尺度不一定相同，将每个网格点设置为成像点，用基于地震走时扰动插值的射线追踪方法计算每个成像点到达地面各个接收点的旅行时间。
步骤三：利用第二步中得到某一个成像点η到各个接收点i的旅行走时对相应的微地震记录进行走时校正，消除该成像点到不同接收点的时差，然后对时差校正后的微地震记录进行两两互相关，计算对应的互相关系数。
利用第二步中某一个成像点η的旅行走时作为微震信号到达各个地面接收点i的旅行走时t ηi ，并设微震信号的发震时刻T=t 0 ，则到达时刻为t=t ηi + t 0 ，确定时间窗口宽度W，则任意两个地震道的波形互相关系数为：
其中，W为所取时间窗口内的采样点个数，u i 为第i个地面接收点的微震记录幅值；u j 为第j个地面接收点微震记录的幅值；t 0 为发震时刻；t ηi 为地下成像点η到第i个地面接收点的旅行走时；t ηj 为地下成像点η到第j个接收点的旅行走时；t k 为所取时间窗口内的第k个采样点时间。
步骤四：将步骤三中得到的波形互相关系数进行相乘计算该成像点的互相关能量值E，计算公式为：
其中，N为地震道数，E（η，t 0 ）表示成像点η在发震时刻T=t 0 的成像能量值。
步骤五：针对每一个成像点重复步骤（3）和步骤（4），获得地下所有成像点的互相关能量值E。
步骤六：设置一个合适的能量阈值E 0 ，将步骤（4）中所得到的所有成像点的能量值E与能量阈值E 0 进行比较，通过搜寻出大于所设能量阈值的能量极大值成像点（E（η，t 0 ）> E 0 ）；即所有符合要求的成像点η为微震事件发生位置，所设的发震时刻t 0 即为微震事件的发震时刻。
下面通过具体的实施例对本发明的效果及可靠性作进一步说明。
实例一：模拟实验
如图2所示，模型区域范围为2000m*2000m*4400m；速度模型由5层水平速度层组成，由上到下的速度值依次为1500m/s，2400m/s，3000m/s，4000m/s和5000m/s；深度依次为1200m，2000m，2400m，3200m和4400m；地面共布设了50个不规则分布的接收点，假设理论微震源位置坐标为（950m，950m，2725m），发震时刻t=0s。图3是通过褶积计算得到的合成微震记录，图4是加入一定程度的随机噪声后的微震记录，从图4中可以看出微震记录信噪比较低，微震事件信号识别不出来。对于地下成像点的划分，设置X和Y方向的成像区域范围为500m~1000m，Z方向的成像区域范围为2600m~2800m，X、Y和Z三个方向的成像点间隔均为25m。使用常规叠加的微震定位方法的震源空间成像结果如图5，图6和图7所示。由于微震记录信噪比低，初至极性发生反转和台站布设等因素的影响，常规方法的震源空间成像分辨率较低，出现了“拖尾”现象，导致不易识别出微震源发生位置。根据本发明前述流程，进行微震事件进行偏移成像定位，震源空间成像结果如图8所示。能量极大值的成像点位置明显可见，表明了震源空间成像分辨率较高，且该能量极大值所对应的成像点位置与理论微震源位置完全一致。
实例二：实际资料
图9是某地区油田的地面微震监测台站位置示意图。采用了不规则分布的27个台站进行监测，将无效台站进行剔除后，选用18个有效台站进行射孔定位。图10为第一段射孔信号，其中可识别出射孔信号的台站有第7和第13个台站，其他台站所接收到的射孔信号均被噪音所淹没，射孔地震资料信噪比较低。根据本发明前述流程，进行射孔定位：根据井口位置坐标和射孔深度，划分X方向的成像区域范围为700m~1000m，Y方向的成像区域范围为400m~0.70km，Z方向的成像区域范围为2560m~2680m；X、Y和Z方向的成像区域间隔分别为10m，10m和20m；选用信噪比较高的第5、第7和第13个台站记录作为主地震道。图11为常规叠加的微震定位方法在Z方向上的水平成像结果示意图，从图11可见，归一化能量值在垂直方向上有逐渐集中的趋势，并且在Z=2.62km处的水平切片出现能量极大值，但是成像分辨率较低。图12（a）和图12（b）是常规叠加的微震定位方法在XZ面和YZ面的成像结果图，从图12（a）和图12（b）可得到，能量极大值成像点位置不易分辨，坐标约为（850m，590m，2620m），与射孔坐标（866.58m，588.94m，2620m）较为吻合。图13是基于波形互相关系数相乘的微震偏移成像定位方法在Z方向上的水平成像结果示意图。从图13可见，能量仅仅在Z=2620m的水平切片出现极大值，且能量值极大值与周围空间的能量值相差较大，差值范围约为1，易于确定微震源位置。图14（a）和图14（b）是基于波形互相关系数相乘的微震偏移成像定位方法在XZ面和YZ面的成像结果图。从图14（a）和图14（b）可知，能量极大值成像点清晰可见，易于得到该能量极大值成像点的位置坐标为（860m，590m，2620m），相对于常规叠加的微震扫描定位方法所得到的射孔位置更接近于射孔位置，满足了精度要求，且空间成像分辨率更高。
以上所述仅为本发明的较佳实施例而己，并不以本发明为限制，凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等，均应包含在本发明的专利涵盖范围内。
Claims (5)
1.一种基于波形互相关系数相乘的微震偏移成像定位方法，其特征在于，包括如下步骤：
（1）对实际微震数据进行预处理，得到有效地震道的微震数据；
（2）将监测的地下区域按一定尺度进行网格划分，每个网格设置一个成像点；利用射线追踪方法求得地下所有成像点到各个地面接收点的旅行走时；
（3）利用第二步中得到某一个成像点η到各个接收点i的旅行走时对相应的微地震记录进行走时校正，消除该成像点到不同接收点的时差，然后对时差校正后的微地震记录进行两两互相关，计算对应的互相关系数；
（4）将成像点对应的所有波形互相关系数进行相乘，得到该成像点的互相关能量值E；
（5）针对每一个成像点重复步骤（3）和步骤（4），获得地下所有成像点的互相关能量值E；
（6）设置能量阈值E 0 ，则超过所设阈值的能量极大值成像点所对应的空间位置和时间参数，作为代表微震源的发生位置和发震时刻。
2.根据权利要求1所述的基于波形互相关系数相乘的微震偏移成像定位方法，其特征在于，步骤（2）中网格划分是三维网格化，在x, y, z三个坐标方向的进行划分，三个方向上的划分尺度不一定相同，将每个网格点设置为一个成像点，用基于地震走时扰动插值的射线追踪方法计算每个成像点到达地面各个接收点的旅行时间。
3.根据权利要求1所述的基于波形互相关系数相乘的微震偏移成像定位方法，其特征在于，步骤（3）中得到波形互相关系数，具体做法为：
利用第二步中某一个成像点η的旅行走时作为微震信号到达各个地面接收点i的旅行走时t ηi ，并设微震信号的发震时刻T=t 0 ，则到达时刻为t=t ηi + t 0 ，确定时间窗口宽度W，则任意两个地震道的波形互相关系数为：
其中，W为所取时间窗口内的采样点个数，u i 为第i个地面接收点的微震记录幅值；u j 为第j个地面接收点微震记录的幅值；t 0 为发震时刻；t ηi 为地下成像点η到第i个地面接收点的旅行走时；t ηj 为地下成像点η到第j个接收点的旅行走时；t k 为所取时间窗口内的第k个采样点时间。
4.根据权利要求1所述的基于波形互相关系数相乘的微震偏移成像定位方法，其特征在于，步骤（4）中将得到的波形互相关系数进行相乘计算得到该成像点的互相关能量值E，具体做法为：
其中，N为地震道数，E（η，t 0 ）表示成像点η在发震时刻T=t 0 的成像能量值。
5.根据权利要求1所述的基于波形互相关系数相乘的微震偏移成像定位方法，其特征在于，步骤（6）中确定微震事件发生位置，具体做法为：
设置一个合适的能量阈值E 0 ，将步骤（4）中所得到的所有成像点的能量值E与能量阈值E 0 进行比较，通过搜寻出大于所设能量阈值的能量极大值成像点（E（η，t 0 ）> E 0 ）；即所有符合要求的成像点η为微震事件发生位置，所设的发震时刻t 0 即为微震事件的发震时刻。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810366402.1A CN110389377B (zh) | 2018-04-23 | 2018-04-23 | 基于波形互相关系数相乘的微震偏移成像定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810366402.1A CN110389377B (zh) | 2018-04-23 | 2018-04-23 | 基于波形互相关系数相乘的微震偏移成像定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110389377A true CN110389377A (zh) | 2019-10-29 |
| CN110389377B CN110389377B (zh) | 2020-11-27 |
Family
ID=68284391
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810366402.1A Expired - Fee Related CN110389377B (zh) | 2018-04-23 | 2018-04-23 | 基于波形互相关系数相乘的微震偏移成像定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110389377B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111190740A (zh) * | 2020-01-02 | 2020-05-22 | 江南大学 | 一种用于波形互相关双差定位算法的并行优化方法 |
| CN111257928A (zh) * | 2020-02-12 | 2020-06-09 | 中国科学技术大学 | 一种震源的定位方法及定位系统 |
| CN111538086A (zh) * | 2020-06-05 | 2020-08-14 | 吉林大学 | 一种提高地震数据初至波质量的方法 |
| CN115951403A (zh) * | 2022-12-29 | 2023-04-11 | 南方科技大学 | 微震偏移叠加定位方法、装置、电子设备及可读存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005006020A1 (en) * | 2003-07-05 | 2005-01-20 | Westerngeco Seismic Holdings Limited | Passive seismic event detection |
| CN105242313A (zh) * | 2015-09-06 | 2016-01-13 | 中国科学院地质与地球物理研究所 | 一种弹性波逆时偏移极性反转的校正方法及系统 |
| US20160306058A1 (en) * | 2015-04-16 | 2016-10-20 | The Board Of Trustees Of The Leland Stanford Junior Univerisity | Reverse time migration based on geometric mean for imaging seismic sources |
| CN106353821A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种微地震事件定位方法 |
-
2018
- 2018-04-23 CN CN201810366402.1A patent/CN110389377B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005006020A1 (en) * | 2003-07-05 | 2005-01-20 | Westerngeco Seismic Holdings Limited | Passive seismic event detection |
| US20160306058A1 (en) * | 2015-04-16 | 2016-10-20 | The Board Of Trustees Of The Leland Stanford Junior Univerisity | Reverse time migration based on geometric mean for imaging seismic sources |
| CN106353821A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种微地震事件定位方法 |
| CN105242313A (zh) * | 2015-09-06 | 2016-01-13 | 中国科学院地质与地球物理研究所 | 一种弹性波逆时偏移极性反转的校正方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| CONGCONG YUAN ET AL.: "Automatic microseismic stacking location with a multi-cross-correlation condition", 《SEG 2017 WORKSHOP: MICROSEISMIC TECHNOLOGIES AND APPLICATIONS HEFEI, CHINA》 * |
| 刘翰林等: "地震自动识别及震相自动拾取方法研究进展", 《地球物理学进展》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111190740A (zh) * | 2020-01-02 | 2020-05-22 | 江南大学 | 一种用于波形互相关双差定位算法的并行优化方法 |
| CN111257928A (zh) * | 2020-02-12 | 2020-06-09 | 中国科学技术大学 | 一种震源的定位方法及定位系统 |
| CN111257928B (zh) * | 2020-02-12 | 2021-08-13 | 中国科学技术大学 | 一种震源的定位方法及定位系统 |
| CN111538086A (zh) * | 2020-06-05 | 2020-08-14 | 吉林大学 | 一种提高地震数据初至波质量的方法 |
| CN111538086B (zh) * | 2020-06-05 | 2021-08-10 | 吉林大学 | 一种提高地震数据初至波质量的初至自动拾取方法 |
| CN115951403A (zh) * | 2022-12-29 | 2023-04-11 | 南方科技大学 | 微震偏移叠加定位方法、装置、电子设备及可读存储介质 |
| CN115951403B (zh) * | 2022-12-29 | 2023-08-18 | 南方科技大学 | 微震偏移叠加定位方法、装置、电子设备及可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110389377B (zh) | 2020-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110389377A (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
| CN105093292B (zh) | 一种地震成像的数据处理方法和装置 | |
| CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
| CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
| CN108445533B (zh) | 基于obs的长偏移距构建浅水区速度模型的方法 | |
| CN107315194B (zh) | 地震叠前时间偏移速度的确定方法和装置 | |
| CN107843922A (zh) | 一种基于地震初至波和反射波走时联合的层析成像方法 | |
| CN106249295B (zh) | 一种井中微地震p、s波联合快速定位方法及系统 | |
| CN101923175B (zh) | 波动方程偏移直接产生角道集方法 | |
| CN104483705A (zh) | 一种三维剩余静校正方法 | |
| CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
| CN109856677A (zh) | 一种震电联合获取裂缝信息的定位方法 | |
| CN109444955A (zh) | 三维地震射线追踪的双线性走时扰动插值方法 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN109856679A (zh) | 一种各向异性介质弹性波高斯束偏移成像方法及系统 | |
| CN104570073B (zh) | 一种适用于复杂高陡构造的双反射地震波成像方法 | |
| CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
| CN108303736A (zh) | 各向异性ti介质最短路径射线追踪正演方法 | |
| CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
| CN104614762B (zh) | 疏松砂岩气藏边界确定方法及装置 | |
| MX2010006060A (es) | Estimacion de parametros elasticos sub-superficiales. | |
| CN104216007B (zh) | 二维垂直地震剖面和三维地面数据的空间同步成像处理方法 | |
| CN109358363B (zh) | 基于多次反射等效界面模型的三维鸣震追踪与预测方法 | |
| CN109143362A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201127 Termination date: 20210423 |