CN105467440A - 一种全向矢量地震数据处理方法及装置 - Google Patents
一种全向矢量地震数据处理方法及装置 Download PDFInfo
- Publication number
- CN105467440A CN105467440A CN201510716250.XA CN201510716250A CN105467440A CN 105467440 A CN105467440 A CN 105467440A CN 201510716250 A CN201510716250 A CN 201510716250A CN 105467440 A CN105467440 A CN 105467440A
- Authority
- CN
- China
- Prior art keywords
- wave
- data
- shear wave
- vector
- omnidirectional
- 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
- 239000013598 vector Substances 0.000 title claims abstract description 256
- 238000003672 processing method Methods 0.000 title abstract description 3
- 238000003384 imaging method Methods 0.000 claims abstract description 89
- 238000000034 method Methods 0.000 claims abstract description 82
- 230000008569 process Effects 0.000 claims abstract description 35
- 238000000926 separation method Methods 0.000 claims abstract description 11
- 238000012937 correction Methods 0.000 claims description 30
- 230000003068 static effect Effects 0.000 claims description 30
- 238000012545 processing Methods 0.000 claims description 27
- 238000003908 quality control method Methods 0.000 claims description 24
- 238000010276 construction Methods 0.000 claims description 23
- 238000004458 analytical method Methods 0.000 claims description 21
- 238000013508 migration Methods 0.000 claims description 15
- 230000005012 migration Effects 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 abstract description 17
- 238000011161 development Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 abstract description 8
- 238000007781 pre-processing Methods 0.000 abstract 2
- 230000001737 promoting effect Effects 0.000 abstract 1
- 238000011084 recovery Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 14
- 210000000988 bone and bone Anatomy 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 239000003990 capacitor Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000003534 oscillatory effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000004035 construction material Substances 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 201000006549 dyspepsia Diseases 0.000 description 1
- 239000013536 elastomeric material Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000000741 silica gel Substances 0.000 description 1
- 229910002027 silica gel Inorganic materials 0.000 description 1
- 230000002123 temporal effect Effects 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
- G01V1/284—Application of the shear wave component and/or several components of the seismic signal
-
- 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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
- G01V1/181—Geophones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
- G01V1/181—Geophones
- G01V1/184—Multi-component geophones
-
- 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/14—Quality control
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- 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/67—Wave propagation modeling
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
本发明公开了一种全向矢量地震数据处理方法及装置,应用于地震全向矢量检波器。其中,该方法包括:采集地震全向矢量检波器的全向矢量地震数据,对全向矢量地震数据执行预处理操作;对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;对纵波数据和横波数据相继执行空间矢量计算和波场恢复操作、成像操作,再建模得到纵波速度体和横波速度体。通过本发明,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。
Description
技术领域
本发明涉及地震勘探技术领域,尤其涉及一种全向矢量地震数据处理方法及装置。
背景技术
根据波动理论,波动不只表示振动,还表示振动的传播。具体地:1、振动可分解为线振动、旋振动;2、波场散度驱动纵波,波场旋度驱动横波;3、线振动是散度和旋度的共同作用,既包含纵波,也包含横波;4、旋振动只与旋度有关;5、体旋度是完全旋度,面旋度是不完全旋度。因此,只能检测线振动的技术,不能干净分离纵横波。必须是能够检测体旋度或散度的技术,才有可能解出纯横波和纯纵波。
波动的空间运动属性,包含了丰富的信息,可以在波场分离、信噪比、保真度、成像精度、介质属性分析等方面起到重要作用。但是,现有采集技术却只能检测振幅、频率、相位等信息,基本没有检测波动空间运动属性的功能。
目前地震波的检测,是将波振动转换为电信号(电压、电流)或再进行数字转换来实现。将机械运动转换为电信号的方法,不外乎应用电磁检波器、电容检波器、压电检波器和光纤应变检波器。
电磁检波器和电容检波器是线振动型,有工作方向。理想方向滤波效应为cosθ,out(t)=A(t)cosθ。A(t)和θ都是未知的,仅单个器件是求不出真振幅A(t)和夹角θ的。MEMS检波器基本属于电容型检波器,也遵循上述规律。
光纤应变检波器可归于线振动或压强型,单器件也不能确定振动方向。所以,单个电磁检波器、电容检波器和光纤应变检波器,都不能区分波场振动方向,更不具备检测波场散度旋度的功能。
压电检波器是压强型,输出与周围液态介质压强有关,无方向性,不能区分振动方向。液态环境中压强各向同性,等效于波场散度。但在陆上固态环境,即使将其置于液态容器中,也无法实现散度测量。
三分量检波器,以三矢量合成方法,可以求出波场振动方向θ和振幅A(t)。这也是称其为矢量检波器的原因,但只是测量点的线振动矢量,不能检测振动性质、旋度和散度。
综上所述,现有地震勘探技术不能测量地震波场的散度和全旋度,不能解决波场振动方向、传播方向、真振幅和纵横波分离的问题,致使现有地震勘探技术的地层成像、地层岩性和流体检测精度不能完全满足资源勘探需求,具体表现在以下两个方面:
1)对地震波动性质研究分析的限制
由于对地震波场的信息采集不全,导致纵波、横波分解不开,有用信息和干扰信息的界定模糊不清。干扰源的分析基于假设和海量试验,但效果不佳。例如,陆上提高保真度和山地改善信噪比的采集处理方法相互制约,噪音压制效果有待进一步提高。
2)限制了波动理论,尤其是全弹理论的进一步研究
弹性波理论的应用,通过全弹波动方程的研究、数据采集、实验、应用,能够解出介质弹性参数,准确确定介质属性及其空间分布状态,从而提高岩性和流体识别精度。利用地震波动的空间运动属性,例如:振动性质、散度、旋度等介质弹性参数,开展弹性波研究,是提高资源勘探精度的技术发展方向。但是,现有地震勘探技术无法准确测量地震波场的散度数据和旋度数据,从而限制了波动理论,尤其是全弹理论的进一步研究。
现有传统地震勘探技术在野外利用垂直振动的检波器,记录地震波场在垂直方向投影的振幅、到达时间,称之为纵波勘探,在某些情况下利用水平放置的检波器,记录地震波场在水平方向投影的振幅和到达时间,称之为转换波勘探,也称之为横波勘探。这种观测方式存在纵波和横波不纯等问题。对观测到的纵波、横波分别和联合对纵横波进行成像处理,分别和联合对纵横波进行属性拾取和分析、反演。由于观测到的地震波记录中纵横波混杂,相互干涉,造成陆上纵波信噪比低,横波不易识别,地震成像精度受到一定影响,岩性预测和流体预测存在一定的多解性。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明提供了一种全向矢量地震数据处理方法及装置,以至少解决现有地震勘探技术不能测量及处理地震波场的散度数据和旋度数据的问题。
根据本发明的一个方面,提供了一种全向矢量地震数据处理方法,应用于地震全向矢量检波器,其中,所述地震全向矢量检波器包括:八个检波器及支撑结构,所述支撑结构用于支撑所述八个检波器,使得每两个检波器的底面位于一正四面体的其中一个正三角形面上,位于同一正三角形面上的两个检波器的工作轴与该正三角形面的交点均位于该正三角形面的一个角的角平分线上并且相对于该正三角形面的中心对称;
所述方法包括:
采集所述地震全向矢量检波器的全向矢量地震数据,对所述全向矢量地震数据执行预处理操作;对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;对所述纵波数据和所述横波数据执行空间矢量计算和波场恢复操作;其中,所述波场恢复操作包括:静校正操作和反褶积操作;基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到纵波速度体和横波速度体。
优选地,对所述全向矢量地震数据执行预处理操作,包括:利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;归类存储符合质控条件的全向矢量地震数据。
优选地,对执行预处理操作后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据,包括:基于所述符合质控条件的全向矢量地震数据,求出总线振动量;根据所述地震全向矢量检波器的旋量数据求出横波线振动分量,根据所述横波线振动分量得到横波数据;将所述总线振动量减去所述横波线振动分量,得到纵波线振动分量,根据所述纵波线振动分量得到纵波数据。
优选地,对所述纵波数据和所述横波数据执行空间矢量计算,包括:在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;将所述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。
优选地,对所述纵波数据和所述横波数据执行波场恢复操作,包括:基于所述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;基于所述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;分方向根据所述纵波地表方向一致性静校正量、所述纵波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的纵波数据;分方向根据所述横波地表方向一致性静校正量、所述横波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的横波数据。
优选地,基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,包括:基于执行所述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;基于执行所述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。
优选地,基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到所述纵波速度体和所述横波速度体,包括:基于所述纵波成像数据在时间域上建立纵波速度界面,基于所述横波成像数据在时间域上建立横波速度界面;在所述纵波速度界面中输入纵波初始速度,得到纵波速度模型;在所述横波速度界面中输入横波初始速度,得到横波速度模型;基于所述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于所述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;利用所述纵波数据CIP道集对执行所述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用所述横波数据CIP道集对执行所述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
根据本发明的另一个方面,提供了一种全向矢量地震数据处理装置,应用于地震全向矢量检波器,其特征在于,所述地震全向矢量检波器包括:八个检波器及支撑结构,所述支撑结构用于支撑所述八个检波器,使得每两个检波器的底面位于一正四面体的其中一个正三角形面上,位于同一正三角形面上的两个检波器的工作轴与该正三角形面的交点均位于该正三角形面的一个角的角平分线上并且相对于该正三角形面的中心对称;
所述装置包括:
数据采集模块,用于采集所述地震全向矢量检波器的全向矢量地震数据全向矢量地震数据;预处理模块,用于对所述全向矢量地震数据执行预处理操作;纵横波分离模块,用于对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;空间矢量计算模块,用于对所述纵波数据和所述横波数据执行空间矢量计算;波场恢复模块,用于对执行空间矢量计算之后的纵波数据和横波数据执行波场恢复操作;其中,所述波场恢复操作包括:静校正操作和反褶积操作;成像模块,用于基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;建模模块,用于基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到纵波速度体和横波速度体。
优选地,所述预处理模块包括:质控单元,用于利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;存储单元,用于归类存储符合质控条件的全向矢量地震数据。
优选地,所述纵横波分离模块包括:横波分离单元,用于基于所述符合质控条件的全向矢量地震数据,求出总线振动量;根据所述地震全向矢量检波器的旋量数据求出横波线振动分量,根据所述横波线振动分量得到横波数据;纵波分离单元,用于将所述总线振动量减去所述横波线振动分量,得到纵波线振动分量,根据所述纵波线振动分量得到纵波数据。
优选地,所述空间矢量计算模块包括:纵波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;横波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;将所述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。
优选地,所述波场恢复模块包括:纵波计算单元,用于基于所述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;分方向根据所述纵波地表方向一致性静校正量、所述纵波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的纵波数据;横波计算单元,用于基于所述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;分方向根据所述横波地表方向一致性静校正量、所述横波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的横波数据。
优选地,所述成像模块包括:纵波成像单元,用于基于执行所述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;横波成像单元,用于基于执行所述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。
优选地,所述建模模块包括:界面建立单元,用于基于所述纵波成像数据在时间域上建立纵波速度界面,基于所述横波成像数据在时间域上建立横波速度界面;模型建立单元,用于在所述纵波速度界面中输入纵波初始速度,得到纵波速度模型;在所述横波速度界面中输入横波初始速度,得到横波速度模型;叠前深度偏移操作单元,用于基于所述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于所述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;速度体生成单元,用于利用所述纵波数据CIP道集对执行所述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用所述横波数据CIP道集对执行所述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
本发明提供了一种全向矢量地震数据处理方法及装置,应用地震全向矢量检波器采集全向矢量地震数据为基础,进行地震波的散度和旋度分离,通过散度和旋度处理分析,最终形成不含横波振动信号(旋度)的纯纵波信息和不含胀缩振动信号(散度)的纯横波信息,分别对纵波和横波进行成像处理,进而进行纵、横波联合反演、解释等工作,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的限定。在附图中:
图1是根据本发明实施例的地震全向矢量检波器的结构示意图;
图2是根据本发明实施例的地震全向矢量检波器的第一种结构示意图;
图3是根据本发明实施例的地震全向矢量检波器的第二种结构示意图;
图4是根据本发明实施例的地震全向矢量检波器的第一种结构的工作矢量示意图;
图5是根据本发明实施例的地震全向矢量检波器的第二种结构的工作矢量示意图;
图6是根据本发明实施例的平行矢量示意图;
图7是根据本发明实施例的地震全向矢量检波器的空间矢量关系图;
图8是根据本发明实施例的全积分等效高斯散度结构示意图;
图9是根据本发明实施例的斯托克斯积分正反旋等效结构示意图;
图10是根据本发明实施例的地震全向矢量检波器实测输出示意图;
图11是根据本发明实施例的地震全向矢量检波器的球形壳体外观示意图;
图12是根据本发明实施例的地震全向矢量检波器的外观示意图;
图13是根据本发明实施例的全向矢量地震数据处理方法的流程图;
图14是根据本发明实施例的全向矢量地震数据处理装置的结构示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
传统技术中的地震波检测流程如下:
根据波动方程,介质质点在力场中的运动关系如下式(1):
其中,λ、μ是拉梅常数,ρ是密度,U是位移向量,t是时间。
是外力向量,
是拉普拉斯算子,
θ是体变系数,
对公式(1)求散度,得到以下公式(2):
其中,div是散度,Vp是纵波传播速度。
对公式(1)求旋度,得到以下公式(3)有:
其中,rot是旋度。w=rot(U),Vs是横波传播速度。
在传统技术中,地震波的检测只是采集到在检波器的竖直工作方向的投影,根本就不是之后就根据公式(2)演化各种方程解纵波;采集到在检波器的水平工作方向的投影,之后就根据公式(3)演化各种方程解横波。在这些误差的基础上进行联合求解、联合反演。因此,致使最终的数据处理结果存在一定误差。
本发明的技术方案就是利用地震全向矢量检波器采集质点振动的旋度,即在采集点的响应,采集线振动矢量,就是在采集点的响应,之后求出散度,就是在采集点的响应。上述是本发明技术方案的核心。下面首先对上述地震全向矢量检波器的结构进行介绍,再结合实施例对本发明基于地震全向矢量检波器进行全向矢量地震数据处理的技术方案进行介绍。
实施例一
本发明提出的地震全向矢量检波器,其在空间结构上是各向均等平衡的,这是现有技术中检波器都不具备的结构优势,利用这个结构优势可实现地震波场全信息的准确检测。
图1是根据本发明实施例的地震全向矢量检波器的结构示意图,如图1所示,该地震全向矢量检波器包括:八个检波器及支撑结构,将八个检波器分成四组,两个检波器为一组,该支撑结构用于支撑上述八个检波器,使得每一组的两个检波器的底面位于一正四面体的其中一个正三角形面上,位于同一正三角形面上的两个检波器的工作轴与该正三角形面的交点均位于该正三角形面的一个角的角平分线上并且相对于该正三角形面的中心对称。
需要说明的是,上述工作轴为虚拟轴,用以说明位置关系。在实际应用中,检波器的工作轴并不一定位于检波器的正中央位置(即中心轴)。上述正四面体并非为实体结构,而是虚拟的空间结构,旨在描述八个检波器之间的位置关系。
本实施例根据场论的散度和旋度公式,设计特定的空间运动全矢量检测结构,以实现对地震波场的频率、振幅、相位、振动方向、尤其是波动力场的散度和旋度的全信息检测,形成全新的地震全向矢量检波器技术。
地震全向矢量检波器的支撑结构可以基于需求设计成任意形状,本发明对此不做限定,只要能保证空间矢量形态即可。优选地,可以设计为正四面体形状,正四面体支撑结构的尺寸一般大于上述虚拟的正四面体。优选地,正四面体支撑结构的每个面可以改变为任意曲面,从而满足不同需求。支撑结构的材料不局限于固态,支撑结构的刚度与测量结果相关,因此可基于需求选择不同材料。
在本实施例中,正四面体的每个正三角形面上有两个检波器,这两个检波器在正三角形面的摆放位置,需要符合特定要求,才能保证地震全向矢量检波器的正常工作。对于每个正三角形面上两个检波器的位置,至少有两种摆放方式,下面分别进行介绍。
1)图2是根据本发明实施例的地震全向矢量检波器的第一种结构示意图,图2以支撑结构为正四面体形状为例进行介绍,如图2所示,设正四面体的其中一个正三角形面为主三角形面100,主三角形面100其中一个角定义为正四面体的第一顶角a,主三角形面100上的两个检波器的工作轴与主三角形面100的交点均位于第一顶角a的角平分线上;
第一顶角a对应的底边所在的正四面体的另一面定义为正四面体的第一侧面101,该底边的中点在主三角形面100逆时针方向滑动,遇到的第一侧面101上的角定义为第二顶角b,第一侧面101上的两个检波器的工作轴与第一侧面101的交点均位于第二顶角b的角平分线上;
需要说明的是,上述逆时针方向是在主三角形面100水平朝上放置时,人面对第一侧面101时的逆时针方向,下述的逆时针方向及顺时针方向均采用此方式。
第二顶角b对应的底边所在的正四面体的另一面定义为正四面体的第二侧面102,该底边的中点在第一侧面101顺时针方向滑动,遇到的第二侧面102上的角定义为第三顶角c,第二侧面102上的两个检波器的工作轴与第二侧面102的交点均位于第三顶角c的角平分线上;
第三顶角c对应的底边所在的正四面体的另一面定义为正四面体的第三侧面103,该底边的中点在第二侧面102逆时针方向滑动,遇到的第三侧面103上的角定义为第四顶角d,第三侧面103上的两个检波器的工作轴与第三侧面103的交点均位于第四顶角d的角平分线上。
2)图3是根据本发明实施例的地震全向矢量检波器的第二种结构示意图,图3以支撑结构为正四面体形状为例进行介绍,如图3所示,上述第一顶角a、上述第二顶角b、上述第三顶角c及上述第四顶角d分别对应正四面体的第一顶点A、第二顶点B、第三顶点C及第四顶点D。
主三角形面100上的两个检波器的工作轴与主三角形面100的交点均位于第一顶角a的角平分线上;
第一侧面101上的两个检波器的工作轴与第一侧面101的交点均位于第三顶点C在第一侧面101上的角的角平分线上;
第二侧面102上的两个检波器的工作轴与第二侧面102的交点均位于第四顶点D在第二侧面102上的角的角平分线上;
第三侧面103上的两个检波器的工作轴与第三侧面103的交点均位于第二顶点B在第三侧面103上的角的角平分线上。
对于上述两种摆放方式,还可以采取下述方法实现:
1)在正四面体的四个正三角形面上,分别设置连接一个角、面几何中心和对边棱(即对应的底边)中点的唯一的一条直线。正四面体的一个面朝上设置(即上述主三角形面),体几何中心与此面的几何中心的连线竖直向上定义为地震全向矢量检波器的竖直轴(也可以称为主轴),也是第一线矢量轴,朝上设置的面上设置的上述直线的方向,定义为地震全向矢量检波器的自北方向(也可以称为第一方向)。沿与自北方向垂直相交的棱(即底边)逆时针滑动遇到的第一侧面上的角,指向对棱中点的方向定义为第二方向,沿与第二方向垂直相交的棱顺时针滑动遇到的第二侧面上的角,指向对棱中点的方向定义为第三方向,沿与第三方向垂直相交的棱逆时针滑动遇到的第三侧面上的角,指向对棱中点的方向定义为第四方向。
沿上述四个方向的方向线上,在每个正三角面的中心向角和棱分别移动a(其中,0<a<λ/4,λ是波长)各安装一个工作方向垂直于正三角面的高灵敏度宽频带全倾角检波器。
2)在上述朝上设置的面上,沿与上述自北方向垂直相交的棱顺时针滑动遇到的角,指向对棱中点的方向定义为第五方向,沿与第五方向垂直相交的棱逆时针滑动遇到的角,指向对棱中点的方向定义为第六方向,沿与第六方向垂直相交的棱顺时针滑动遇到的角,指向对棱中点的方向定义为第七方向。
沿上述四个方向的方向线上,在每个正三角面的中心向角和棱分别移动a(其中,0<a<λ/4,λ是波长)各安装一个工作方向垂直于正三角面的高灵敏度宽频带全倾角检波器。
本实施例中的检波器至少可以是圆柱电磁型检波器,需满足:各个圆柱电磁型检波器的底面与正四面体的体几何中心的距离为各个检波器互不接触的最小值;还可以是扁片电容型检波器。当然,也可以是其他类型的检波器,只要能够在地震全向矢量检波器中正常工作即可,本发明对此不做限制。
图4是根据本发明实施例的地震全向矢量检波器的第一种结构的工作矢量示意图,图5是根据本发明实施例的地震全向矢量检波器的第二种结构的工作矢量示意图,图4和图5中的每个面上的箭头表示工作矢量,与主三角形面100垂直的中心轴设置为主轴200,主三角形面100上第一顶角a的角平分线的方向设置为自北方向。
图6是根据本发明实施例的平行矢量示意图,如图6所示,在光滑连续波场中,间距远小于波长的两个平行矢量与波场的点乘之和除以2,等于中间位置的矢量与波场的点乘。
图7是根据本发明实施例的地震全向矢量检波器的空间矢量关系图,如图7所示,图中的各个箭头即表示每个面上的工作矢量。图7中的圆表示地震全向矢量检波器的八个空间矢量所围成的内切圆球。
图8是根据本发明实施例的全积分等效高斯散度结构示意图,图9是根据本发明实施例的斯托克斯积分正反旋等效结构示意图,结合图8和图9所示的矢量关系可知:
根据散度高斯定理:可得本发明实施例的地震全向矢量检波器的散度公式:
根据旋度斯托克斯公式:可得本实施例的正旋度公式:以及,反旋度公式:
其中,A为波场函数,l为检波器响应矢量,i为检波器的内道序号,为采集到的第i道内部道数据,在上述全积分等效高斯散度结构上:为第一组,i=1-4,为第二组,i=5-8。div为散度,rot为旋度,V为体积,dv为体积微分,S为面积,dS为面积微分,m为正m面体的个数;ni为第i面法线上的检波器响应矢量。根据三角和差化积公式,在本特定结构上,
需要说明的是,本发明提供的地震全向矢量检波器是基于虚拟正四面体的四根轴(即正四面体的体中心与四个正三角面的中心连线构成的立体星型轴)形成的四均分立体空间平衡结构。各个检波器的底面距体中心距离d,d为单元器件互不接触的最小值,在正四面体的每个正三角形面上分别与轴平行安装两个高灵敏度宽频带全倾角检波器,其振幅响应与自身工作轴和波场法线间夹角θ的关系为out(t)=A(t)cosθ,其中,A(t)为检波器在θ=0时对波场的响应。
在地震全向矢量检波器中,每个正三角形面上靠近角的检波器为第一组,靠近棱(即对边)的检波器为第二组,第一组或第二组的检波器在空间中构成符合斯托克斯旋度公式的封闭环路。第一组输出为正旋度,第二组输出为负旋度。
位于同一正三角形面上的两个检波器的算术和,在a小于λ/4的情况下,等效于位于每个正三角面的中心处的检波器,形成四个指向正四面体的体中心的测量矢量,与高斯散度公式的要求完全吻合。
本实施例完美地将斯托克斯旋度结构和高斯散度结构结合于一体,形成了旋散结合的全矢量地震测量空间结构。
地震全向矢量检波器中包括八个检波器,第一组和第二组的检波器共计8道数据输出。图10是根据本发明实施例的地震全向矢量检波器实测输出示意图,如图10所示,四个上述第一电路的输出结果,对应于图10中所示的第一组旋量,四个上述第二电路的输出结果,对应于图10中所示的第二组旋量。第一组旋量与第二组旋量的和为0,第一组旋量与第二组旋量的差为双旋量。地震全向矢量检波器的矢量方向由1旋向2右手法则确定,总旋矢量方向可经矢量运算求出。振动线矢量可以由冗余矢量计算求出。等同于传统三分量检波器计算出的振动线矢量。应用旋度求解横波振动线矢量的方法,可以求出纯横波线矢量,从总振动线矢量中减去纯横波线矢量,可得纯纵波振动线矢量。
前面已经对地震全向矢量检波器的具体结构进行了详细描述,在实际使用过程中,为了方便放置,地震全向矢量检波器的支撑结构上连接一圆椎形尾椎,该圆椎形尾椎的轴线的延长线经过正四面体的一个顶点后,垂直穿过正四面体的其中一个面的中心。在正常使用过程中,该圆锥形尾椎竖直朝下插入地面,从而便于固定。
地震全向矢量检波器还可以设置球形壳体,分为上半球形壳体和下半球形壳体,八个检波器及支撑结构放置在球形壳体内部,下半球形壳体的底部设置尾椎孔,上述尾椎穿过下半球形壳体的尾椎孔。图11是根据本发明实施例的地震全向矢量检波器的球形壳体外观示意图,如图11所示,地震全向矢量检波器的球形壳体分为上半球形壳体20和下半球形壳体21,地震全向矢量检波器的支撑结构上连接的圆椎形尾椎30穿过下半球形壳体21的底部设置的尾椎孔31。通过该优选实施方式,可以保护检波器,稳固地震全向矢量检波器。
为了将地震全向矢量检波器的支撑结构及其上的八个检波器稳固放在球形壳体中,可以将支撑结构的形状设计为易于放入球形壳体的形状。图12是根据本发明实施例的地震全向矢量检波器的外观示意图,如图12所示,将正四面体形状的支撑结构的各个顶点削掉一块,以便于安置在球形壳体中,确保地震全向矢量检波器采集数据的准确性。
考虑到上半球形壳体和下半球形壳体的固定问题,本实施例提供了一种优选实施方式,即:上半球形壳体和下半球形壳体的边缘分别设置相互配合的突出部,上半球形壳体和下半球形壳体的突出部通过固定组件固定。例如,通过螺丝和螺丝孔的方式固定。从而保证地震全向矢量检波器在实际操作过程中的稳定性。
对于地震全向矢量检波器的信号输出线,可以在球形壳体上设置信号线孔,使检波器的信号输出线穿过该信号线孔,从而便于信号输出线输出到壳体外面,有效获取数据。
为了避免地震全向矢量检波器进水影响使用,可以在球形壳体的上半球形壳体和下半球形壳体之间的接合缝隙、尾椎孔和信号线孔,均密封防水,例如,可以选择硅胶或橡胶材料进行密封,从而提高地震全向矢量检波器的密封性。
本发明提供的地震全向矢量检波器可以适用于以下方面:陆上人工地震勘探、天然地震探测。本发明的应用前景如下:由于基于本发明的地震全向矢量检波器可检测到地震波的旋度,并可求出波场的振动方向及真振幅,可分解出纯净的纵横波,因此可以获得以往从未有过的丰富的地震波信息,为地球科学的探索研究奠定了全新的数据基础。
地震全向矢量检波器的散度和旋度特征,为刚性支撑结构的检波器放置于地表自由表面时的特征。当检波器放置于地下或是介质内部时,这些特征将由放置环境、检波器支撑结构材料及其与外界的接触关系而确定。可随具体应用目的的变化而设计。本发明不做详细说明。
前面对地震全向矢量检波器的结构和工作原理进行了详细介绍,下面介绍基于地震全向矢量检波器进行全向矢量地震数据处理的技术方案。
实施例二
本实施例提供了一种全向矢量地震数据处理方法,应用于上述地震全向矢量检波器,图13是根据本发明实施例的全向矢量地震数据处理方法的流程图,如图13所示,该流程包括以下步骤(步骤S102-步骤S110):
步骤S102,采集上述地震全向矢量检波器的全向矢量地震数据,对上述全向矢量地震数据执行预处理操作;
步骤S104,对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;
步骤S106,对上述纵波数据和上述横波数据执行空间矢量计算和波场恢复操作;其中,上述波场恢复操作包括:静校正操作和反褶积操作;
步骤S108,基于执行上述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;
步骤S110,基于执行上述波场恢复操作之后的纵波数据和横波数据,以及上述纵波成像数据和上述横波成像数据进行建模,得到纵波速度体和横波速度体。
本实施例应用地震全向矢量检波器采集全向矢量地震数据为基础,对地震波的纵波数据和横波数据进行处理,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。
一实施例中,在步骤S102中,对全向矢量地震数据执行预处理操作,具体过程是:利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;归类存储符合质控条件的全向矢量地震数据。基于此,可根据需求设置质控条件,从而对全向矢量地震数据进行筛选,利于后续的数据处理。
例如:质控条件可以是:位于正四面体的每个面上的两个检波器中,将靠近其所在的角平分线对应的角的检波器归为第一组,远离其所在的角平分线对应的角的检波器归为第二组。设第一组的四个检波器的内道序号分别为:1、2、3、4,第二组的四个检波器的内道序号分别为:5、6、7、8。按照内道序号1、2、3、4、5、6、7、8的顺序,检测检波器的内方位角是否依次为:0°、0°、120°、-120°、0°、0°、120°、-120°;其中,内方位角为自北方向与自然北极方向的夹角,自北方向是正四面体的水平朝上设置的面上的两个检波器所在角平分线的方向。按照内道序号1、2、3、4、5、6、7、8的顺序,检测检波器的内倾角是否依次为:90°、-30°、-30°、-30°、90°、-30°、-30°、-30°;其中,内倾角为检波器所在平面与自然水平面的夹角。
一实施例中,在步骤S104中,对执行预处理操作后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据,具体过程是:基于符合质控条件的全向矢量地震数据,求出总线振动量;根据上述地震全向矢量检波器的旋量数据求出横波线振动分量,根据上述横波线振动分量得到横波数据;将上述总线振动量减去上述横波线振动分量,得到纵波线振动分量,根据上述纵波线振动分量得到纵波数据。基于此,可准确分离纵波数据和横波数据,解决了现有地震勘探技术不能测量地震波场的散度数据和全旋度数据的问题。
一实施例中,在步骤S106中,对纵波数据和横波数据执行空间矢量计算,具体过程是:在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;将上述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。基于此,对纵波数据和横波数据进行空间矢量计算,为后续执行波场恢复操作提供基础。
一实施例中,在步骤S106中,对纵波数据和横波数据执行波场恢复操作,具体过程是:基于上述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;基于上述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;分方向根据上述纵波地表方向一致性静校正量、上述纵波地表方向一致性反褶积因子,生成执行上述波场恢复操作之后的纵波数据;分方向根据上述横波地表方向一致性静校正量、上述横波地表方向一致性反褶积因子,生成执行上述波场恢复操作之后的横波数据。基于此,对纵波数据和横波数据进行波场恢复处理,为后续执行成像操作提供基础。
一实施例中,在步骤S108中,基于执行波场恢复操作之后的纵波数据和横波数据执行成像操作,具体过程是:基于执行上述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;基于执行上述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。基于此,为后续建模提供数据支撑。
一实施例中,在步骤S110中,基于执行波场恢复操作之后的纵波数据和横波数据,以及纵波成像数据和横波成像数据进行建模,得到纵波速度体和横波速度体,具体过程是:基于纵波成像数据在时间域上建立纵波速度界面,基于横波成像数据在时间域上建立横波速度界面;在纵波速度界面中输入纵波初始速度,得到纵波速度模型;在横波速度界面中输入横波初始速度,得到横波速度模型;基于上述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于上述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;利用上述纵波数据CIP道集对执行上述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用上述横波数据CIP道集对执行上述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
本实施例应用地震全向矢量检波器采集全向矢量地震数据为基础,进行地震波的散度和旋度分离,通过散度和旋度处理分析,最终形成不含横波振动信号(旋度)的纯纵波信息和不含胀缩振动信号(散度)的纯横波信息,分别对纵波和横波进行成像处理,进而进行纵、横波联合反演、解释等工作,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。
实施例三
基于同一发明构思,本实施例中还提供了一种全向矢量地震数据处理装置,可以用于实现上述实施例所描述的方法,如下面的实施例所述。由于全向矢量地震数据处理装置解决问题的原理与全向矢量地震数据处理方法相似,因此全向矢量地震数据处理装置的实施可以参见全向矢量地震数据处理方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图14是根据本发明实施例的全向矢量地震数据处理装置的结构示意图,该装置应用于地震全向矢量检波器,前面已经对地震全向矢量检波器的结构进行了介绍,在此不再赘述。如图14所示,该装置包括:数据采集模块12、预处理模块14、纵横波分离模块16、空间矢量计算模块18、波场恢复模块22、成像模块24、建模模块26。下面对各个模块的功能进行介绍。
数据采集模块12,用于采集上述地震全向矢量检波器的全向矢量地震数据全向矢量地震数据;
预处理模块14,连接至数据采集模块12,用于对上述全向矢量地震数据执行预处理操作;
纵横波分离模块16,连接至预处理模块14,用于对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;
空间矢量计算模块18,连接至纵横波分离模块16,用于对上述纵波数据和上述横波数据执行空间矢量计算;
波场恢复模块22,连接至空间矢量计算模块18,用于对执行空间矢量计算之后的纵波数据和横波数据执行波场恢复操作;其中,上述波场恢复操作包括:静校正操作和反褶积操作;
成像模块24,连接至波场恢复模块22,用于基于执行上述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;
建模模块26,连接至成像模块24,用于基于执行上述波场恢复操作之后的纵波数据和横波数据,以及上述纵波成像数据和上述横波成像数据进行建模,得到纵波速度体和横波速度体。
本实施例应用地震全向矢量检波器采集全向矢量地震数据为基础,对地震波的纵波数据和横波数据进行处理,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。
在本实施例中,上述预处理模块14包括:质控单元,用于利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;存储单元,用于归类存储符合质控条件的全向矢量地震数据。基于此,可根据需求设置质控条件,从而对全向矢量地震数据进行筛选,利于后续的数据处理。
在本实施例中,上述纵横波分离模块16包括:横波分离单元,用于基于上述符合质控条件的全向矢量地震数据,求出总线振动量;根据上述地震全向矢量检波器的旋量数据求出横波线振动分量,根据上述横波线振动分量得到横波数据;纵波分离单元,用于将上述总线振动量减去上述横波线振动分量,得到纵波线振动分量,根据上述纵波线振动分量得到纵波数据。基于此,可准确分离纵波数据和横波数据,解决了现有地震勘探技术不能测量地震波场的散度数据和全旋度数据的问题。
在本实施例中,上述空间矢量计算模块18包括:纵波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;横波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;将上述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。基于此,对纵波数据和横波数据进行空间矢量计算,为后续执行波场恢复操作提供基础。
在本实施例中,上述波场恢复模块22包括:纵波计算单元,用于基于上述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;分方向根据上述纵波地表方向一致性静校正量、上述纵波地表方向一致性反褶积因子,生成执行上述波场恢复操作之后的纵波数据;横波计算单元,用于基于上述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;分方向根据上述横波地表方向一致性静校正量、上述横波地表方向一致性反褶积因子,生成执行上述波场恢复操作之后的横波数据。基于此,对纵波数据和横波数据进行波场恢复处理,为后续执行成像操作提供基础。
在本实施例中,上述成像模块24包括:纵波成像单元,用于基于执行上述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;横波成像单元,用于基于执行上述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。基于此,为后续建模提供数据支撑。
在本实施例中,上述建模模块26包括:界面建立单元,用于基于上述纵波成像数据在时间域上建立纵波速度界面,基于上述横波成像数据在时间域上建立横波速度界面;模型建立单元,用于在上述纵波速度界面中输入纵波初始速度,得到纵波速度模型;在上述横波速度界面中输入横波初始速度,得到横波速度模型;叠前深度偏移操作单元,用于基于上述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于上述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;速度体生成单元,用于利用上述纵波数据CIP道集对执行上述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用上述横波数据CIP道集对执行上述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
当然,上述模块划分只是一种示意划分,本发明并不局限于此。只要能实现本发明的目的的模块划分,均应属于本发明的保护范围。
为了完整的采集地震波场的振幅、频率、相位、振动方向、传播方向、质点力场的散度和旋度,开展地震波场全息研究,本发明提出了全矢量地震数据采集、处理、成像的地震全向矢量勘探技术,可达到以下目的:
1、在传统采集技术的基础上,增加采集地震波的散度、旋度、线矢量、旋矢量的功能。
2、比传统处理技术提高纵横波分离的纯净度。
3、提高陆上,尤其是复杂山地采集的地震资料的信噪比。
4、提高地震资料的保真度。
5、以全矢量地震数据为输入,为应用全弹性波地震成像和弹性参数提取奠定基础。
6、为全信息速度反演(Fullwaveinformationinversion,FWI)、岩性等反演技术应用奠定基础。
7、实现全矢量地震属性的提取,提高岩石、流体预测的精度和可靠性。
传统地震勘探利用地下介质弹性和密度的差异,在地表以人工方法激发地震波,在向地下传播时,遇有介质性质不同的岩层分界面,地震波将发生反射与折射,在地表或井中用检波器接收这种地震波。通过对地震波记录进行处理和解释,可以推断地下岩层的性质和形态。源于光学反射理论,从上世纪二十年代至今,经历了二维、三维观测的技术发展之路。
本发明的技术方案源于对地震波场复杂性的考虑,传统地震观测方式难以解决横波分离、难以记录地震波场的全析信息,为了提高陆上复杂地表地震资料的信噪比,获得更丰富的地震波场信息,从数据采集源头上记录波的散度、旋度、方向、大小、时间信息,是基于传统地震勘探技术,又对传统地震勘探技术革命性跨越式发展。
本发明的关键特征是在野外利用全向检波器记录全向波场信息,通过特殊的地震信号处理手段,进行纵波、横波分离,分别和联合对纵横波进行成像处理,提高地震信号的信噪比、灵敏度、保真度,提高构造成像精度,分别和联合对纵横波进行属性拾取和分析、反演,提高地震属性的物理信息量,提高岩性和流体预测精度。
本发明应用地震全向矢量检波器采集全向矢量地震数据为基础,进行地震波的散度和旋度分离,通过散度和旋度处理分析,最终形成不含横波振动信号(旋度)的纯纵波信息和不含胀缩振动信号(散度)的纯横波信息,分别对纵波和横波进行成像处理,进而进行纵、横波联合反演、解释等工作,解决了现有地震勘探技术不能测量及处理地震波场的散度数据和全旋度数据的问题,从而提高了构造、岩性、流体勘探精度和可靠性,推动了地震勘探从构造勘探向岩性勘探和流体勘探方向发展。本发明的理论基础、数据性质、结构、格式、信息量等和现有地震勘探技术有很大不同,是一套全新的地震勘探技术系列,代表了地震勘探技术的发展方向。
本发明适用于陆上各类资源的人工地震勘探、天然地震监测、自然灾害监测等领域。本发明的应用前景如下:由于可检测到地震波的旋度,可求出波场的振动方向及真振幅,分解出纯净的纵横波,可以获得以往从未有过的丰富的地震波信息,还原波场的本质面貌,不仅提高地震信号的信噪比、灵敏度、保真度,还提供了更加丰富的地震属性,必将提高资源勘探的精度、天然气地震强度、烈度监测的可靠性,并可在矿山资源生产安全监测中发挥作用。该技术的进一步的研究和应用,必将推广地球物理技术革命性的发展,为我国复杂地下矿产资源勘探、自然灾害动态监测、生产安全监测发挥更大的作用。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种全向矢量地震数据处理方法,应用于地震全向矢量检波器,其特征在于,所述地震全向矢量检波器包括:八个检波器及支撑结构,所述支撑结构用于支撑所述八个检波器,使得每两个检波器的底面位于一正四面体的其中一个正三角形面上,位于同一正三角形面上的两个检波器的工作轴与该正三角形面的交点均位于该正三角形面的一个角的角平分线上并且相对于该正三角形面的中心对称;
所述方法包括:
采集所述地震全向矢量检波器的全向矢量地震数据,对所述全向矢量地震数据执行预处理操作;
对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;
对所述纵波数据和所述横波数据执行空间矢量计算和波场恢复操作;其中,所述波场恢复操作包括:静校正操作和反褶积操作;
基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;
基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到纵波速度体和横波速度体。
2.根据权利要求1所述的全向矢量地震数据处理方法,其特征在于,对所述全向矢量地震数据执行预处理操作,包括:
利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;
归类存储符合质控条件的全向矢量地震数据。
3.根据权利要求2所述的全向矢量地震数据处理方法,其特征在于,对执行预处理操作后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据,包括:
基于所述符合质控条件的全向矢量地震数据,求出总线振动量;
根据所述地震全向矢量检波器的旋量数据求出横波线振动分量,根据所述横波线振动分量得到横波数据;
将所述总线振动量减去所述横波线振动分量,得到纵波线振动分量,根据所述纵波线振动分量得到纵波数据。
4.根据权利要求1所述的全向矢量地震数据处理方法,其特征在于,对所述纵波数据和所述横波数据执行空间矢量计算,包括:
在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;
在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;
将所述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。
5.根据权利要求4所述的全向矢量地震数据处理方法,其特征在于,对所述纵波数据和所述横波数据执行波场恢复操作,包括:
基于所述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;
基于所述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;
分方向根据所述纵波地表方向一致性静校正量、所述纵波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的纵波数据;
分方向根据所述横波地表方向一致性静校正量、所述横波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的横波数据。
6.根据权利要求1所述的全向矢量地震数据处理方法,其特征在于,基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,包括:
基于执行所述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;
基于执行所述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。
7.根据权利要求1所述的全向矢量地震数据处理方法,其特征在于,基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到所述纵波速度体和所述横波速度体,包括:
基于所述纵波成像数据在时间域上建立纵波速度界面,基于所述横波成像数据在时间域上建立横波速度界面;
在所述纵波速度界面中输入纵波初始速度,得到纵波速度模型;在所述横波速度界面中输入横波初始速度,得到横波速度模型;
基于所述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于所述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;
利用所述纵波数据CIP道集对执行所述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用所述横波数据CIP道集对执行所述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
8.一种全向矢量地震数据处理装置,应用于地震全向矢量检波器,其特征在于,所述地震全向矢量检波器包括:八个检波器及支撑结构,所述支撑结构用于支撑所述八个检波器,使得每两个检波器的底面位于一正四面体的其中一个正三角形面上,位于同一正三角形面上的两个检波器的工作轴与该正三角形面的交点均位于该正三角形面的一个角的角平分线上并且相对于该正三角形面的中心对称;
所述装置包括:
数据采集模块,用于采集所述地震全向矢量检波器的全向矢量地震数据全向矢量地震数据;
预处理模块,用于对所述全向矢量地震数据执行预处理操作;
纵横波分离模块,用于对执行预处理操作之后的全向矢量地震数据执行纵横波分离操作,得到纵波数据和横波数据;
空间矢量计算模块,用于对所述纵波数据和所述横波数据执行空间矢量计算;
波场恢复模块,用于对执行空间矢量计算之后的纵波数据和横波数据执行波场恢复操作;其中,所述波场恢复操作包括:静校正操作和反褶积操作;
成像模块,用于基于执行所述波场恢复操作之后的纵波数据和横波数据执行成像操作,得到纵波成像数据和横波成像数据;
建模模块,用于基于执行所述波场恢复操作之后的纵波数据和横波数据,以及所述纵波成像数据和所述横波成像数据进行建模,得到纵波速度体和横波速度体。
9.根据权利要求8所述的全向矢量地震数据处理装置,其特征在于,所述预处理模块包括:
质控单元,用于利用全向矢量观测系统,检查采集到的全向矢量地震数据是否符合质控条件;
存储单元,用于归类存储符合质控条件的全向矢量地震数据。
10.根据权利要求9所述的全向矢量地震数据处理装置,其特征在于,所述纵横波分离模块包括:
横波分离单元,用于基于所述符合质控条件的全向矢量地震数据,求出总线振动量;根据所述地震全向矢量检波器的旋量数据求出横波线振动分量,根据所述横波线振动分量得到横波数据;
纵波分离单元,用于将所述总线振动量减去所述横波线振动分量,得到纵波线振动分量,根据所述纵波线振动分量得到纵波数据。
11.根据权利要求8所述的全向矢量地震数据处理装置,其特征在于,所述空间矢量计算模块包括:
纵波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的纵波数据;
横波数据计算单元,用于在矢量方向上,根据方位角及倾角分选输出不同方向的横波数据;将所述不同方向的横波数据减去地面横波数据,得到不同方向的有效横波数据。
12.根据权利要求11所述的全向矢量地震数据处理装置,其特征在于,所述波场恢复模块包括:
纵波计算单元,用于基于所述不同方向的纵波数据,分方向计算纵波地表方向一致性静校正量及纵波地表方向一致性反褶积因子;分方向根据所述纵波地表方向一致性静校正量、所述纵波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的纵波数据;
横波计算单元,用于基于所述不同方向的有效横波数据,分方向计算横波地表方向一致性静校正量及横波地表方向一致性反褶积因子;分方向根据所述横波地表方向一致性静校正量、所述横波地表方向一致性反褶积因子,生成执行所述波场恢复操作之后的横波数据。
13.根据权利要求8所述的全向矢量地震数据处理装置,其特征在于,所述成像模块包括:
纵波成像单元,用于基于执行所述波场恢复操作之后的纵波数据,形成纵波成像剖面或纵波成像三维体;
横波成像单元,用于基于执行所述波场恢复操作之后的横波数据,形成横波成像剖面或横波成像三维体。
14.根据权利要求8所述的全向矢量地震数据处理装置,其特征在于,所述建模模块包括:
界面建立单元,用于基于所述纵波成像数据在时间域上建立纵波速度界面,基于所述横波成像数据在时间域上建立横波速度界面;
模型建立单元,用于在所述纵波速度界面中输入纵波初始速度,得到纵波速度模型;在所述横波速度界面中输入横波初始速度,得到横波速度模型;
叠前深度偏移操作单元,用于基于所述纵波速度模型进行叠前深度偏移操作,得到纵波数据CIP道集;基于所述横波速度模型进行叠前深度偏移操作,得到横波数据CIP道集;
速度体生成单元,用于利用所述纵波数据CIP道集对执行所述波场恢复操作之后的纵波数据进行迭代速度分析,得到纵波速度体;利用所述横波数据CIP道集对执行所述波场恢复操作之后的横波数据进行迭代速度分析,得到横波速度体。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510716250.XA CN105467440B (zh) | 2015-10-28 | 2015-10-28 | 一种全向矢量地震数据处理方法及装置 |
PCT/CN2016/099825 WO2017071438A1 (zh) | 2015-10-28 | 2016-09-23 | 全向矢量地震数据处理方法、装置、计算机可读存储介质及设备 |
EP16858874.7A EP3370091B1 (en) | 2015-10-28 | 2016-09-23 | Omni-directional vector geophone seismic data processing method and apparatus |
US15/965,503 US10578757B2 (en) | 2015-10-28 | 2018-04-27 | Omnidirectional vector seismic data processing method and apparatus, computer-readable storage medium and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510716250.XA CN105467440B (zh) | 2015-10-28 | 2015-10-28 | 一种全向矢量地震数据处理方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105467440A true CN105467440A (zh) | 2016-04-06 |
CN105467440B CN105467440B (zh) | 2018-02-02 |
Family
ID=55605337
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510716250.XA Active CN105467440B (zh) | 2015-10-28 | 2015-10-28 | 一种全向矢量地震数据处理方法及装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10578757B2 (zh) |
EP (1) | EP3370091B1 (zh) |
CN (1) | CN105467440B (zh) |
WO (1) | WO2017071438A1 (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017071438A1 (zh) * | 2015-10-28 | 2017-05-04 | 中国石油天然气股份有限公司 | 全向矢量地震数据处理方法、装置、计算机可读存储介质及设备 |
CN109100800A (zh) * | 2018-06-29 | 2018-12-28 | 安徽万泰地球物理技术有限公司 | 一种基于微地震成像的围岩稳定性预警方法和系统 |
CN110187008A (zh) * | 2019-04-29 | 2019-08-30 | 中国科学院地理科学与资源研究所 | 一种无损声波全弹性参数探测装置及方法 |
CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110824427A (zh) * | 2019-10-21 | 2020-02-21 | 中国人民解放军陆军勤务学院 | 一种斜三棱锥声压测试装置及其空间次声源定向定位方法 |
CN111025386A (zh) * | 2019-12-13 | 2020-04-17 | 成都理工大学 | 一种无分离假象的纵横波分离方法 |
US10718875B2 (en) | 2015-10-28 | 2020-07-21 | Petrochina Company Limited | Omnidirectional vector electrostatic levitation geophone |
US10739477B2 (en) | 2015-10-28 | 2020-08-11 | Petrochina Company Limited | Omnidirectional vector geophone |
CN117574245A (zh) * | 2024-01-16 | 2024-02-20 | 电子科技大学 | 一种应用于山地勘探的智能检波器指标自检方法及系统 |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3670830B1 (en) | 2016-04-07 | 2021-08-11 | BP Exploration Operating Company Limited | Detecting downhole events using acoustic frequency domain features |
US11530606B2 (en) | 2016-04-07 | 2022-12-20 | Bp Exploration Operating Company Limited | Detecting downhole sand ingress locations |
CA3058256C (en) | 2017-03-31 | 2023-09-12 | Bp Exploration Operating Company Limited | Well and overburden monitoring using distributed acoustic sensors |
WO2019038401A1 (en) | 2017-08-23 | 2019-02-28 | Bp Exploration Operating Company Limited | DETECTION OF SAND INPUT LOCATIONS AT THE BOTTOM OF A HOLE |
CN111771042A (zh) | 2017-10-11 | 2020-10-13 | 英国石油勘探运作有限公司 | 使用声学频域特征来检测事件 |
CN110879413A (zh) * | 2018-09-05 | 2020-03-13 | 中国石油化工股份有限公司 | 一种射线参数域转换波静校正方法及系统 |
US20210389486A1 (en) | 2018-11-29 | 2021-12-16 | Bp Exploration Operating Company Limited | DAS Data Processing to Identify Fluid Inflow Locations and Fluid Type |
GB201820331D0 (en) | 2018-12-13 | 2019-01-30 | Bp Exploration Operating Co Ltd | Distributed acoustic sensing autocalibration |
CN112327353A (zh) * | 2019-07-30 | 2021-02-05 | 中国石油天然气集团有限公司 | 一种横波折射层建模方法及装置 |
EP4045766A1 (en) | 2019-10-17 | 2022-08-24 | Lytt Limited | Fluid inflow characterization using hybrid das/dts measurements |
WO2021073740A1 (en) | 2019-10-17 | 2021-04-22 | Lytt Limited | Inflow detection using dts features |
WO2021093974A1 (en) | 2019-11-15 | 2021-05-20 | Lytt Limited | Systems and methods for draw down improvements across wellbores |
CN113138407B (zh) * | 2020-01-20 | 2024-05-03 | 中国石油天然气集团有限公司 | 一种深层页岩气的多尺度断裂地震预测方法及系统 |
CN111597677A (zh) * | 2020-04-01 | 2020-08-28 | 王茂发 | 一种测震图像的矢量化方法、装置、设备及存储介质 |
WO2021249643A1 (en) | 2020-06-11 | 2021-12-16 | Lytt Limited | Systems and methods for subterranean fluid flow characterization |
CA3182376A1 (en) | 2020-06-18 | 2021-12-23 | Cagri CERRAHOGLU | Event model training using in situ data |
CN111948710A (zh) * | 2020-08-05 | 2020-11-17 | 河海大学 | 一种利用海底多分量地震记录进行弹性矢量波成像的方法 |
CN114137609B (zh) * | 2021-11-09 | 2023-12-01 | 长江地球物理探测(武汉)有限公司 | 线性微动数据校正方法及装置 |
US12000264B2 (en) | 2021-12-14 | 2024-06-04 | Saudi Arabian Oil Company | Method for quadrimodal fault prediction using strain tensor cyclides |
CN114814944B (zh) * | 2022-03-15 | 2022-10-21 | 长江大学 | 基于散度和旋度的弹性纵横波场分离和逆时偏移成像方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6791901B1 (en) * | 1998-09-16 | 2004-09-14 | Schlumberger Technology Corporation | Seismic detection apparatus and related method |
JP2006010520A (ja) * | 2004-06-25 | 2006-01-12 | Takamisawa Cybernetics Co Ltd | 感震器及び地震計システム |
US20110063946A1 (en) * | 2004-01-30 | 2011-03-17 | Westerngeco Llc | Marine seismic acquisition system |
CN102053261A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震数据处理方法 |
CN103645497A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 一种基于出射角的多分量波场分离方法 |
Family Cites Families (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE266183C (zh) | 1910-06-13 | |||
JPS6313729A (ja) | 1986-07-05 | 1988-01-21 | Sumitomo Heavy Ind Ltd | ラミネ−ト装置の自動運転システム |
GB8616904D0 (en) * | 1986-07-11 | 1986-08-20 | Texas Instruments Ltd | Motion sensing device |
CN2047790U (zh) | 1989-03-12 | 1989-11-15 | 西安石油学院 | 悬浮式涡流地震检波器 |
GB2291074B (en) | 1994-06-30 | 1998-04-15 | Gec Alsthom Ltd | Static electrical converter arrangement |
FR2733831B1 (fr) | 1995-05-04 | 1997-06-20 | Inst Francais Du Petrole | Capteur de vibrations |
GB0001355D0 (en) * | 2000-01-21 | 2000-03-08 | Geco Prakla Uk Ltd | Wavefield separation and estimation of near surface properties in land seismic |
US6291901B1 (en) | 2000-06-13 | 2001-09-18 | ćEFO NEVRES | Electrical power generating tire system |
CN2453455Y (zh) | 2000-12-08 | 2001-10-10 | 威海双丰电子传感有限公司 | 三分量地震检波器 |
GB2372568B (en) * | 2001-02-26 | 2004-05-12 | Abb Offshore Systems Ltd | Seismic detection |
GB2381314B (en) | 2001-10-26 | 2005-05-04 | Westerngeco Ltd | A method of and an apparatus for processing seismic data |
CN1256597C (zh) | 2003-05-23 | 2006-05-17 | 中国石油化工股份有限公司石油勘探开发研究院南京石油物探研究所 | 三分量数字地震勘探检波器 |
US20100116059A1 (en) | 2004-07-26 | 2010-05-13 | Spider Technologies Security Ltd. | Vibration sensor having a single virtual center of mass |
KR20070062974A (ko) | 2004-07-26 | 2007-06-18 | 스파이더 테크놀로지스 시큐어리티 리미티드 | 진동센서 |
GB2419410B (en) * | 2004-10-20 | 2008-05-21 | Vetco Gray Controls Ltd | Sensor configuration |
CN100449316C (zh) | 2006-08-13 | 2009-01-07 | 重庆大学 | 一种六轴加速度传感器的敏感元件的布局方法 |
CN201043915Y (zh) | 2007-05-23 | 2008-04-02 | 山东大学 | 三分量光纤光栅振动传感器 |
US9304216B2 (en) | 2009-02-05 | 2016-04-05 | Westerngeco L.L.C. | Seismic acquisition system and technique |
US8485027B2 (en) | 2009-09-14 | 2013-07-16 | Schlumberger Technology Corporation | Borehole force measurement |
DE102009046807B4 (de) | 2009-11-18 | 2023-01-05 | Robert Bosch Gmbh | Verfahren zur Empfindlichkeitsbestimmung eines Beschleunigungs- oder Magnetfeldsensors |
CN101893720B (zh) | 2010-07-02 | 2012-09-05 | 中国科学院地质与地球物理研究所 | 一种地震波的矢量波场分离与合成的方法和系统 |
US8040754B1 (en) * | 2010-08-27 | 2011-10-18 | Board Of Regents Of The University Of Texas System | System and method for acquisition and processing of elastic wavefield seismic data |
CN102023309B (zh) | 2010-10-29 | 2013-03-27 | 谭成忠 | 一种磁悬浮电磁感应检波器 |
CN201852944U (zh) | 2010-10-29 | 2011-06-01 | 谭成忠 | 一种磁悬浮电磁感应检波器 |
FR2985039B1 (fr) | 2011-12-21 | 2015-07-03 | Cggveritas Services Sa | Noeud sous-marin couple avec l'eau pour des etudes sismiques |
CN102692646B (zh) | 2012-06-19 | 2013-11-13 | 北京多分量地震技术研究院 | 一种三维三分量矢量波场分离的方法和系统 |
CN103513273B (zh) | 2012-06-28 | 2016-06-01 | 中国科学院地质与地球物理研究所 | 陆用四分量数字地震检波器 |
CN103048683B (zh) | 2012-12-12 | 2015-07-01 | 北京科络捷通讯产品有限责任公司 | 磁悬浮检波器 |
US10048395B2 (en) | 2013-02-01 | 2018-08-14 | Westerngeco L.L.C. | Computing a gradient based on differences of plural pairs of particle motion sensors |
CN103675350B (zh) | 2013-12-17 | 2015-11-25 | 北京航天控制仪器研究所 | 一种硅微静电悬浮加速度计 |
CN104237934A (zh) | 2014-09-29 | 2014-12-24 | 中国地质大学 | 差动式压电地震数字检波器 |
CN104808242B (zh) * | 2014-10-20 | 2017-12-19 | 重庆丰银包装材料有限公司 | 一种地震资料成像处理方法 |
CN105467440B (zh) * | 2015-10-28 | 2018-02-02 | 中国石油天然气股份有限公司 | 一种全向矢量地震数据处理方法及装置 |
CN105259566B (zh) | 2015-10-28 | 2018-02-02 | 中国石油天然气股份有限公司 | 一种地震全向矢量检波器 |
CN105388514B (zh) | 2015-10-28 | 2017-12-05 | 中国石油天然气股份有限公司 | 一种地震全向矢量静电悬浮检波器 |
-
2015
- 2015-10-28 CN CN201510716250.XA patent/CN105467440B/zh active Active
-
2016
- 2016-09-23 WO PCT/CN2016/099825 patent/WO2017071438A1/zh active Application Filing
- 2016-09-23 EP EP16858874.7A patent/EP3370091B1/en active Active
-
2018
- 2018-04-27 US US15/965,503 patent/US10578757B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6791901B1 (en) * | 1998-09-16 | 2004-09-14 | Schlumberger Technology Corporation | Seismic detection apparatus and related method |
US20110063946A1 (en) * | 2004-01-30 | 2011-03-17 | Westerngeco Llc | Marine seismic acquisition system |
JP2006010520A (ja) * | 2004-06-25 | 2006-01-12 | Takamisawa Cybernetics Co Ltd | 感震器及び地震計システム |
CN102053261A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震数据处理方法 |
CN103645497A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 一种基于出射角的多分量波场分离方法 |
Non-Patent Citations (3)
Title |
---|
孙歧峰 等: ""多分量地震数据处理技术研究现状"", 《石油勘探与开发》 * |
李志远 等: ""基于散度和旋度纵横波分离方法的改进"", 《地球物理学报》 * |
李澈: ""多波多分量地震勘探技术进展"", 《资源环境与节能减灾》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10578757B2 (en) | 2015-10-28 | 2020-03-03 | Petrochina Company Limited | Omnidirectional vector seismic data processing method and apparatus, computer-readable storage medium and device |
WO2017071438A1 (zh) * | 2015-10-28 | 2017-05-04 | 中国石油天然气股份有限公司 | 全向矢量地震数据处理方法、装置、计算机可读存储介质及设备 |
US10718875B2 (en) | 2015-10-28 | 2020-07-21 | Petrochina Company Limited | Omnidirectional vector electrostatic levitation geophone |
US10739477B2 (en) | 2015-10-28 | 2020-08-11 | Petrochina Company Limited | Omnidirectional vector geophone |
CN109100800A (zh) * | 2018-06-29 | 2018-12-28 | 安徽万泰地球物理技术有限公司 | 一种基于微地震成像的围岩稳定性预警方法和系统 |
CN110187008A (zh) * | 2019-04-29 | 2019-08-30 | 中国科学院地理科学与资源研究所 | 一种无损声波全弹性参数探测装置及方法 |
CN110398776B (zh) * | 2019-08-21 | 2021-03-30 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110824427A (zh) * | 2019-10-21 | 2020-02-21 | 中国人民解放军陆军勤务学院 | 一种斜三棱锥声压测试装置及其空间次声源定向定位方法 |
CN111025386A (zh) * | 2019-12-13 | 2020-04-17 | 成都理工大学 | 一种无分离假象的纵横波分离方法 |
CN111025386B (zh) * | 2019-12-13 | 2020-11-17 | 成都理工大学 | 一种无分离假象的纵横波分离方法 |
CN117574245A (zh) * | 2024-01-16 | 2024-02-20 | 电子科技大学 | 一种应用于山地勘探的智能检波器指标自检方法及系统 |
CN117574245B (zh) * | 2024-01-16 | 2024-03-26 | 电子科技大学 | 一种应用于山地勘探的智能检波器指标自检方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105467440B (zh) | 2018-02-02 |
US10578757B2 (en) | 2020-03-03 |
WO2017071438A1 (zh) | 2017-05-04 |
US20180246238A1 (en) | 2018-08-30 |
EP3370091A1 (en) | 2018-09-05 |
EP3370091A4 (en) | 2019-06-12 |
EP3370091B1 (en) | 2021-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105467440B (zh) | 一种全向矢量地震数据处理方法及装置 | |
CN102341728B (zh) | 地震获取系统和技术 | |
Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
Ferrari et al. | Rock Mass Rating spatial estimation by geostatistical analysis | |
CN105510993B (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
Akingboye et al. | Insight into seismic refraction and electrical resistivity tomography techniques in subsurface investigations | |
CN107045143B (zh) | 一种预测裂缝发育的方法及装置 | |
RU2014111160A (ru) | Гибридная детерминированная геостатистическая модель геологической среды | |
CN104237940A (zh) | 一种基于动力学特征的绕射波成像方法及装置 | |
CA2828425A1 (en) | Sensitivity kernel-based migration velocity analysis in 3d anisotropic media | |
CN105259566A (zh) | 一种地震全向矢量检波器 | |
CN105388515A (zh) | 一种检波器的质控方法及装置 | |
WO2016001697A1 (en) | Systems and methods for geologic surface reconstruction using implicit functions | |
CN105259567B (zh) | 一种地震全向矢量旋度检波器 | |
Vernon et al. | Near-surface scattering effects observed with a high-frequency phased array at Pinyon Flats, California | |
CN103364834A (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
CN107422375B (zh) | 海底节点的共反射点道集的确定方法和装置 | |
CN105388516A (zh) | 一种地震全向矢量散度检波器 | |
CN112305594B (zh) | 非均质储层的油气分布确定方法及系统 | |
CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
CN107085236B (zh) | 最大炮检距的确定方法和装置 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
CN104375180A (zh) | 一种地震数据处理方法、装置及系统 | |
CN106168678B (zh) | 一种裂隙中传播横波的分离方法及系统 | |
Hekmatnejad | Geostatistical modeling of discrete fracture networks for geomechanical applications in heterogeneous fractured media based on the cox-boolean model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |