CN113466943B - 变偏移距vsp地表一致性振幅补偿方法和装置 - Google Patents
变偏移距vsp地表一致性振幅补偿方法和装置 Download PDFInfo
- Publication number
- CN113466943B CN113466943B CN202010242100.0A CN202010242100A CN113466943B CN 113466943 B CN113466943 B CN 113466943B CN 202010242100 A CN202010242100 A CN 202010242100A CN 113466943 B CN113466943 B CN 113466943B
- Authority
- CN
- China
- Prior art keywords
- reservoir section
- angle
- offset vsp
- determining
- variable offset
- 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 64
- 230000005540 biological transmission Effects 0.000 claims abstract description 134
- 230000008859 change Effects 0.000 claims abstract description 82
- 230000010287 polarization Effects 0.000 claims abstract description 51
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 101710194948 Protein phosphatase PhpP Proteins 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 16
- HWGNBUXHKFFFIH-UHFFFAOYSA-I pentasodium;[oxido(phosphonatooxy)phosphoryl] phosphate Chemical compound [Na+].[Na+].[Na+].[Na+].[Na+].[O-]P([O-])(=O)OP([O-])(=O)OP([O-])([O-])=O HWGNBUXHKFFFIH-UHFFFAOYSA-I 0.000 claims description 16
- 235000019832 sodium triphosphate Nutrition 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 28
- 238000012986 modification Methods 0.000 description 9
- 230000004048 modification Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 239000000523 sample Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000004321 preservation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (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
本发明提供了一种变偏移距VSP地表一致性振幅补偿方法和装置，该方法包括：根据零偏移距VSP，确定储层段弹性参数模型；根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。有效解决了变偏移距VSP地表一致性振幅补偿的问题。
Description
技术领域
本发明涉及地球物理勘探中地震数据处理技术领域，尤其涉及一种变偏移距VSP地表一致性振幅补偿方法和装置。
背景技术
近年来井中地震技术发展迅速，变偏移距VSP(即Walkaway VSP)已实现工业化应用。近地表条件的变化，如风化层速度、厚度、地表高程、炮点的变化等，导致地震反射波振幅差异很大。
这种地表非一致性的振幅差异必须在叠前去除，否则影响后续的反褶积、速度分析、剩余静校正、偏移成像等，进而影响叠前、叠后反演的正确性。
现有技术中还没有关于变偏移距VSP地表一致性振幅补偿的技术方案。
发明内容
本发明实施例提供一种变偏移距VSP地表一致性振幅补偿方法，有效解决了变偏移距VSP地表一致性振幅补偿的问题，该方法包括：
根据零偏移距VSP，确定储层段弹性参数模型；
根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。
本发明实施例还提供一种变偏移距VSP地表一致性振幅补偿装置，包括：
储层段弹性参数模型确定模块，用于根据零偏移距VSP，确定储层段弹性参数模型；
变偏移距VSP初至波极化角确定模块，用于根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
储层段理论透射系数相对变化率确定模块，用于根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
储层段实测透射能量相对变化率确定模块，用于统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
地表一致性振幅补偿系数确定模块，用于根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
多检波点统计共炮点地表一致性振幅补偿系数模块，用于根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述变偏移距VSP地表一致性振幅补偿方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现上述变偏移距VSP地表一致性振幅补偿方法的计算机程序。
本发明实施例提供的一种变偏移距VSP地表一致性振幅补偿方法和装置，利用零偏移距VSP建立储层段弹性参数模型，以变偏移距VSP初至波极化角为入射角，利用透射能量相对变化率与理论透射系数相对变化率的差异计算地表一致性振幅补偿系数，有效解决了变偏移距VSP地表一致性振幅补偿的问题。为变偏移距VSP后续处理、成像提供相对振幅保持的输入数据，为叠前反演保留振幅随角度变化信息。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法示意图。
图2为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段纵波速度模型示意图。
图3为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段密度模型示意图。
图4为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的变偏移距VSP初至波极化角示意图。
图5为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段理论透射系数示意图。
图6为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段实测透射能量示意图。
图7为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段多个检波点统计的共炮点地表一致性振幅补偿系数示意图。
图8为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段地表一致性校正后的变偏移距VSP透射能量示意图。
图9为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的地表一致性振幅补偿前的共检波点道集示意图。
图10为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的地表一致性振幅补偿后的共检波点道集示意图。
图11为运行本发明实施的一种变偏移距VSP地表一致性振幅补偿方法的计算机装置示意图。
图12为本发明实施例一种变偏移距VSP地表一致性振幅补偿装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
如图1本发明实施例一种变偏移距VSP地表一致性振幅补偿方法示意图所示，本发明实施例提供一种变偏移距VSP地表一致性振幅补偿方法，有效解决了变偏移距VSP地表一致性振幅补偿的问题，该方法包括：
步骤101：根据零偏移距VSP，确定储层段弹性参数模型；
步骤102：根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
步骤103：根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
步骤104：统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
步骤105：根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
步骤106：根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。
本发明实施例提供的一种变偏移距VSP地表一致性振幅补偿方法，利用零偏移距VSP建立储层段弹性参数模型，以变偏移距VSP初至波极化角为入射角，利用透射能量相对变化率与理论透射系数相对变化率的差异计算地表一致性振幅补偿系数，有效解决了变偏移距VSP地表一致性振幅补偿的问题。为变偏移距VSP后续处理、成像提供相对振幅保持的输入数据，为叠前反演保留振幅随角度变化信息。
本发明实施例提供一种变偏移距VSP地表一致性振幅补偿方法，包括：
根据零偏移距VSP，确定储层段弹性参数模型；根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。
图2为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段纵波速度模型示意图，横坐标为纵波速度(单位：m/s)；纵坐标为深度(单位：m)。如图2所示，在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，储层段弹性参数模型，包括储层段纵波速度模型和储层段密度模型；
按照如下方式，建立储层段纵波速度模型：
其中，Hi是第i个检波器深度，tp,i是第i个纵波初至，Offset是零偏移距VSP炮点偏离井口的距离，tp,vi是第i个检波器的纵波垂直单程时；
Vp,i＝(Hi-Hi-1)/(tp,vi-tp,vi-1)
其中，Hi-1是第i-1个检波器深度，Hi是第i个检波器深度，tp,vi-1是第i-1个检波器的纵波垂直单程时，tp,vi是第i个检波器的纵波垂直单程时，Vp,i是第i层的纵波速度。
前述提到的建立储层段纵波速度模型的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
图3为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段密度模型示意图，横坐标为密度(单位：g/cm3)；纵坐标为深度(单位：m)。如图3所示，在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，根据声波测井的密度资料重采样成VSP深度采样，建立储层段密度模型。
图4为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的变偏移距VSP初至波极化角示意图，横坐标为偏移距(单位：m)；纵坐标为极化角(单位：°)。如图4所示，在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，根据变偏移距VSP记录，确定变偏移距VSP初至波极化角，包括：
三分量变偏移距VSP记录，根据原始ZXY坐标系与PRT坐标系的关系，确定变偏移距VSP初至波极化角。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，按如下方式，确定变偏移距VSP初至波极化角：
H＝Xcosφ+Ysinφ
T＝Xsinφ-Ycosφ
其中，φ是水平旋转角，N是初至时窗内的采样点数，HT是XY水平旋转后的水平径向分量和水平切向分量；
P＝Zcosθ+Hsinθ
R＝Zsinθ-Hcosθ
其中，θ是变偏移距VSP初至波极化角，N是初至时窗内的采样点数，PR是ZH极化旋转后的垂直径向分量和垂直切向分量。
前述提到的确定变偏移距VSP初至波极化角的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率，包括：
根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数；
根据储层段理论透射系数，确定储层段理论透射系数相对变化率。
图5为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段理论透射系数示意图，横坐标为入射角(单位：°)；纵坐标为透射系数(单位：无)。如图5所示，在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，按照如下方式，确定储层段理论透射系数：
其中，Vp是储层段纵波速度模型，Vp,1、Vp,2是界面上、下的纵波速度，ρ是储层段密度模型，ρ1、ρ2是界面上、下的密度，Tpp(θ)是入射角为θ的理论透射系数，θ采用变偏移距VSP初至波极化角代替计算理论透射系数。
前述提到的确定储层段理论透射系数的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，按照如下方式，确定储层段理论透射系数相对变化率：
STPPi＝Tpp(θi)/Tpp(0)
其中，Tpp(θi)是入射角θi的理论透射系数，Tpp(0)是入射角为0的理论透射系数，STPPi是入射角θi的理论透射系数相对变化率。
前述提到的确定储层段理论透射系数相对变化率的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
图6为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段实测透射能量示意图，横坐标为偏移距(单位：m)；纵坐标为振幅能量(单位：无)，如图6所示，在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，统计储层段实测透射能量：
其中，N是初至时窗内的采样点数，Ampi是初至时窗内第i采样点的振幅，TEng是初至时窗的透射波能量。
前述提到的统计储层段实测透射能量的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，确定储层段实测透射能量相对变化率：
MTPPi＝TEng(θi)/TEng(0)
其中，TEng(θi)是入射角θi的实测透射能量，TEng(0)是入射角为0的实测透射能量，MTPPi是入射角θi的实测透射能量相对变化率。
前述提到的确定储层段实测透射能量相对变化率的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，按如下方式，确定地表一致性振幅补偿系数：
SCAi＝MTPPi/STPPi
其中，MTPPi是入射角θi的实测透射能量相对变化率，STPPi是入射角θi的理论透射系数相对变化率，SCAi是入射角θi的地表一致性振幅补偿系数，入射角θi对应着第i个炮点，SCAi即是第i个炮点的地表一致性振幅补偿系数。
前述提到的确定地表一致性振幅补偿系数的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿方法时，在一个实施例中，按如下方式，多检波点统计共炮点地表一致性振幅补偿系数：
其中，N是参与统计的检波点个数，SCAij是第i个炮点的第j个检波点的地表一致性振幅补偿系数，SCASi是第i个炮点的N个检波点统计共炮点地表一致性振幅补偿系数。
前述提到的多检波点统计共炮点地表一致性振幅补偿系数的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
图7为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段多个检波点统计的共炮点地表一致性振幅补偿系数示意图；横坐标为偏移距(单位：m)；纵坐标为振幅补偿系数(单位：无)。图8为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的储层段地表一致性校正后的变偏移距VSP透射能量示意图；横坐标为偏移距(单位：m)；纵坐标为振幅能量(单位：无)。图9为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的地表一致性振幅补偿前的共检波点道集示意图；横坐标为偏移距(单位：m)；纵坐标为时间(单位：ms)。图10为本发明实施例一种变偏移距VSP地表一致性振幅补偿方法的地表一致性振幅补偿后的共检波点道集示意图；横坐标为偏移距(单位：m)；纵坐标为时间(单位：ms)。
如图11所示，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述变偏移距VSP地表一致性振幅补偿方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现上述变偏移距VSP地表一致性振幅补偿方法的计算机程序。
本发明实施例中还提供了一种变偏移距VSP地表一致性振幅补偿装置，如下面的实施例所述。由于该装置解决问题的原理与一种变偏移距VSP地表一致性振幅补偿方法相似，因此该装置的实施可以参见一种变偏移距VSP地表一致性振幅补偿方法的实施，重复之处不再赘述。
如图12所示，本发明实施例还提供一种变偏移距VSP地表一致性振幅补偿装置，包括：
储层段弹性参数模型确定模块1201，用于根据零偏移距VSP，确定储层段弹性参数模型；
变偏移距VSP初至波极化角确定模块1202，用于根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
储层段理论透射系数相对变化率确定模块1203，用于根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
储层段实测透射能量相对变化率确定模块1204，用于统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
地表一致性振幅补偿系数确定模块1205，用于根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
多检波点统计共炮点地表一致性振幅补偿系数模块1206，用于根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段弹性参数模型，包括储层段纵波速度模型和储层段密度模型；
储层段弹性参数模型确定模块，还用于按照如下方式，建立储层段纵波速度模型：
其中，Hi是第i个检波器深度，tp,i是第i个纵波初至，Offset是零偏移距VSP炮点偏离井口的距离，tp,vi是第i个检波器的纵波垂直单程时；
Vp,i＝(Hi-Hi-1)/(tp,vi-tp,vi-1)
其中，Hi-1是第i-1个检波器深度，Hi是第i个检波器深度，tp,vi-1是第i-1个检波器的纵波垂直单程时，tp,vi是第i个检波器的纵波垂直单程时，Vp,i是第i层的纵波速度。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段弹性参数模型确定模块，还用于根据声波测井的密度资料重采样成VSP深度采样，建立储层段密度模型。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，变偏移距VSP初至波极化角确定模块，具体用于：
三分量变偏移距VSP记录，根据原始ZXY坐标系与PRT坐标系的关系，确定变偏移距VSP初至波极化角。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，变偏移距VSP初至波极化角确定模块，还用于按如下方式，确定变偏移距VSP初至波极化角：
H＝Xcosφ+Ysinφ
T＝Xsinφ-Ycosφ
其中，φ是水平旋转角，N是初至时窗内的采样点数，HT是XY水平旋转后的水平径向分量和水平切向分量；
P＝Z cosθ+Hsinθ
R＝Z sinθ-Hcosθ
其中，θ是变偏移距VSP初至波极化角，N是初至时窗内的采样点数，PR是ZH极化旋转后的垂直径向分量和垂直切向分量。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段理论透射系数相对变化率确定模块，具体用于：
根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数；
根据储层段理论透射系数，确定储层段理论透射系数相对变化率。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段理论透射系数相对变化率确定模块，还用于按照如下方式，确定储层段理论透射系数：
其中，Vp是储层段纵波速度模型，Vp,1、Vp,2是界面上、下的纵波速度，ρ是储层段密度模型，ρ1、ρ2是界面上、下的密度，Tpp(θ)是入射角为θ的理论透射系数，θ采用变偏移距VSP初至波极化角代替计算理论透射系数。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段理论透射系数相对变化率确定模块，还用于按照如下方式，确定储层段理论透射系数相对变化率：
STPPi＝Tpp(θi)/Tpp(0)
其中，Tpp(θi)是入射角θi的理论透射系数，Tpp(0)是入射角为0的理论透射系数，STPPi是入射角θi的理论透射系数相对变化率。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段实测透射能量相对变化率确定模块，具体用于按如下方式，统计储层段实测透射能量：
其中，N是初至时窗内的采样点数，Ampi是初至时窗内第i采样点的振幅，TEng是初至时窗的透射波能量。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，储层段实测透射能量相对变化率确定模块，还用于按如下方式，确定储层段实测透射能量相对变化率：
MTPPi＝TEng(θi)/TEng(0)
其中，TEng(θi)是入射角θi的实测透射能量，TEng(0)是入射角为0的实测透射能量，MTPPi是入射角θi的实测透射能量相对变化率。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，地表一致性振幅补偿系数确定模块，具体用于按如下方式，确定地表一致性振幅补偿系数：
SCAi＝MTPPi/STPPi
其中，MTPPi是入射角θi的实测透射能量相对变化率，STPPi是入射角θi的理论透射系数相对变化率，SCAi是入射角θi的地表一致性振幅补偿系数，入射角θi对应着第i个炮点，SCAi即是第i个炮点的地表一致性振幅补偿系数。
在具体实施本发明实施例的一种变偏移距VSP地表一致性振幅补偿装置时，在一个实施例中，多检波点统计共炮点地表一致性振幅补偿系数模块，具体用于按如下方式，多检波点统计共炮点地表一致性振幅补偿系数：
其中，N是参与统计的检波点个数，SCAij是第i个炮点的第j个检波点的地表一致性振幅补偿系数，SCASi是第i个炮点的N个检波点统计共炮点地表一致性振幅补偿系数。
综上，本发明实施例提供的一种变偏移距VSP地表一致性振幅补偿方法和装置，利用零偏移距VSP建立储层段弹性参数模型，以变偏移距VSP初至波极化角为入射角正演模拟检波器处的理论透射系数，以零度透射系数为基准计算各个角度透射系数的理论变化率。统计各个炮点的初至透射波能量，以井口道为基准计算初至透射波的能量变化率，多检波点统计能量变化率与理论变化率的差异，计算各个炮点的地表一致性振幅补偿系数，有效解决了变偏移距VSP地表一致性振幅补偿的问题。为变偏移距VSP后续处理、成像提供相对振幅保持的输入数据，为叠前反演保留振幅随角度变化信息。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.一种变偏移距VSP地表一致性振幅补偿方法，其特征在于，包括：
根据零偏移距VSP，确定储层段弹性参数模型；
根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数；
其中，按照如下方式，确定储层段理论透射系数相对变化率：
STPPk＝Tpp(θk)/Tpp(0)
其中，Tpp(θk)是入射角θk的理论透射系数，Tpp(0)是入射角为0的理论透射系数，STPPk是入射角θk的理论透射系数相对变化率；
其中，按如下方式，确定储层段实测透射能量相对变化率：
MTPPk＝TEng(θk)/TEng(0)
其中，TEng(θk)是入射角θk的实测透射能量，TEng(0)是入射角为0的实测透射能量，MTPPk是入射角θk的实测透射能量相对变化率；
其中，按如下方式，确定地表一致性振幅补偿系数：
SCAk＝MTPPk/STPPk
其中，MTPPk是入射角θk的实测透射能量相对变化率，STPPk是入射角θk的理论透射系数相对变化率，SCAk是入射角θk的地表一致性振幅补偿系数，入射角θk对应着第k个炮点，SCAk即是第k个炮点的地表一致性振幅补偿系数；
按照如下方式，确定储层段理论透射系数：
其中，Vp是储层段纵波速度模型，Vp,1、Vp,2分别是界面上、下的纵波速度，ρ是储层段密度模型，ρ1、ρ2分别是界面上、下的密度，Tpp(θ)是入射角为θ的理论透射系数，θ采用变偏移距VSP初至波极化角代替，计算理论透射系数。
2.如权利要求1所述的方法，其特征在于，根据声波测井的密度资料重采样成VSP深度采样，建立储层段密度模型。
3.如权利要求1所述的方法，其特征在于，根据变偏移距VSP记录，确定变偏移距VSP初至波极化角，包括：
根据三分量变偏移距VSP记录，根据原始XYZ坐标系与PRT坐标系的关系，确定变偏移距VSP初至波极化角；所述PRT坐标系是由XYZ坐标系水平旋转和极化旋转得到。
4.如权利要求3所述的方法，其特征在于，按如下方式，确定变偏移距VSP初至波极化角：
H＝X cosφ+Y sinφ
T＝X sinφ-Y cosφ
其中，φ是水平旋转角，N是初至时窗内的采样点数，H、T分别是XY水平旋转后的水平径向分量和水平切向分量；
P＝Z cosθ+H sinθ
R＝Z sinθ-H cosθ
其中，θ是变偏移距VSP初至波极化角，N是初至时窗内的采样点数，P、R分别是ZH极化旋转后的垂直径向分量和垂直切向分量；X、Y、Z是XYZ坐标下的三分量，Xi是XYZ坐标系下的X分量的第i个采样点，Yi是XYZ坐标系下的Y分量的第i个采样点，Zi是XYZ坐标系下的Z分量的第i个采样点，P、R、T是水平旋转、极化旋转后的三分量，Hi是水平旋转后的H分量的第i个采样点。
5.一种变偏移距VSP地表一致性振幅补偿装置，其特征在于，包括：
储层段弹性参数模型确定模块，用于根据零偏移距VSP，确定储层段弹性参数模型；
变偏移距VSP初至波极化角确定模块，用于根据变偏移距VSP记录，确定变偏移距VSP初至波极化角；
储层段理论透射系数相对变化率确定模块，用于根据储层段弹性参数模型和变偏移距VSP初至波极化角，确定储层段理论透射系数相对变化率；
储层段实测透射能量相对变化率确定模块，用于统计储层段实测透射能量，根据储层段实测透射能量，确定储层段实测透射能量相对变化率；
地表一致性振幅补偿系数确定模块，用于根据储层段理论透射系数相对变化率和储层段实测透射能量相对变化率，确定地表一致性振幅补偿系数；
多检波点统计共炮点地表一致性振幅补偿系数模块，用于根据地表一致性振幅补偿系数，多检波点统计共炮点地表一致性振幅补偿系数；
其中，储层段理论透射系数相对变化率确定模块，用于按照如下方式，确定储层段理论透射系数相对变化率：
STPPk＝Tpp(θk)/Tpp(0)
其中，Tpp(θk)是入射角θk的理论透射系数，Tpp(0)是入射角为0的理论透射系数，STPPk是入射角θk的理论透射系数相对变化率；
其中，储层段实测透射能量相对变化率确定模块，用于按如下方式，确定储层段实测透射能量相对变化率：
MTPPk＝TEng(θk)/TEng(0)
其中，TEng(θk)是入射角θk的实测透射能量，TEng(0)是入射角为0的实测透射能量，MTPPk是入射角θk的实测透射能量相对变化率；
其中，地表一致性振幅补偿系数确定模块，用于按如下方式，确定地表一致性振幅补偿系数：
SCAk＝MTPPk/STPPk
其中，MTPPk是入射角θk的实测透射能量相对变化率，STPPk是入射角θk的理论透射系数相对变化率，SCAk是入射角θk的地表一致性振幅补偿系数，入射角θk对应着第k个炮点，SCAk即是第k个炮点的地表一致性振幅补偿系数；
按照如下方式，确定储层段理论透射系数：
其中，Vp是储层段纵波速度模型，Vp,1、Vp,2分别是界面上、下的纵波速度，ρ是储层段密度模型，ρ1、ρ2分别是界面上、下的密度，Tpp(θ)是入射角为θ的理论透射系数，θ采用变偏移距VSP初至波极化角代替，计算理论透射系数。
6.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述变偏移距VSP地表一致性振幅补偿方法。
7.一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现权利要求1至4任一项所述变偏移距VSP地表一致性振幅补偿方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010242100.0A CN113466943B (zh) | 2020-03-31 | 2020-03-31 | 变偏移距vsp地表一致性振幅补偿方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010242100.0A CN113466943B (zh) | 2020-03-31 | 2020-03-31 | 变偏移距vsp地表一致性振幅补偿方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113466943A CN113466943A (zh) | 2021-10-01 |
| CN113466943B true CN113466943B (zh) | 2023-04-07 |
Family
ID=77865229
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010242100.0A Active CN113466943B (zh) | 2020-03-31 | 2020-03-31 | 变偏移距vsp地表一致性振幅补偿方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113466943B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2831961B1 (fr) * | 2001-11-07 | 2004-07-23 | Inst Francais Du Petrole | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue |
| CN102608659B (zh) * | 2012-03-27 | 2014-11-19 | 中国科学院地质与地球物理研究所 | 一种耦合透射系数的地震偏移方法 |
| CN104375188B (zh) * | 2014-11-19 | 2017-02-08 | 中国石油天然气集团公司 | 一种地震波透射衰减补偿方法及装置 |
| US10852450B2 (en) * | 2017-05-03 | 2020-12-01 | Saudi Arabian Oil Company | Refraction-based surface-consistent amplitude compensation and deconvolution |
| CN108181649B (zh) * | 2017-11-23 | 2019-09-10 | 中国石油天然气集团公司 | 一种介质分界面振幅补偿方法及装置 |
| CN110031897B (zh) * | 2019-03-21 | 2020-08-28 | 北京信息科技大学 | 多分量地震数据在海水中振幅能量补偿与校正方法及系统 |
-
2020
- 2020-03-31 CN CN202010242100.0A patent/CN113466943B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113466943A (zh) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| CN110133715B (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
| US10705238B2 (en) | Method and apparatus for processing seismic data | |
| US10641916B2 (en) | Processing seismic data by nonlinear stacking | |
| EP2739994B1 (en) | Correcting time lapse seismic data for overburden and recording effects | |
| CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
| CN1404581A (zh) | 估算地震介质特性的系统和方法 | |
| WO2012139082A1 (en) | Event selection in the image domain | |
| EP2321671B1 (en) | Processing seismic data in common group-center gathers | |
| CN107678063A (zh) | 一种基于等级相关分析的多分量转换波裂缝预测方法 | |
| Vlček et al. | Back‐projection stacking of P‐and S‐waves to determine location and focal mechanism of microseismic events recorded by a surface array | |
| RU2339056C2 (ru) | Обобщенное трехмерное прогнозирование кратных волн от поверхности | |
| Zhang et al. | Automated microseismic event location by amplitude stacking and semblance | |
| Hunziker et al. | Bayesian full-waveform inversion of tube waves to estimate fracture aperture and compliance | |
| CN101609163B (zh) | 基于波动理论的多尺度地震资料联合成像方法 | |
| US20180164454A1 (en) | Determining residual statics of survey receivers | |
| CN113466943B (zh) | 变偏移距vsp地表一致性振幅补偿方法和装置 | |
| CN110007349B (zh) | 一种弹性参数反演方法 | |
| CN106707345A (zh) | 一种角度弹性参数岩性识别方法和装置 | |
| CN107816348B (zh) | 一种利用纵波和斯通利波识别气层的方法和装置 | |
| US4592032A (en) | Methods of processing geophysical data and displays of geophysical data | |
| CN112415601B (zh) | 表层品质因子q值的确定方法及装置 | |
| CN113267809B (zh) | I类页岩储层预测方法及装置 | |
| CN108802822B (zh) | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 | |
| CN112285781B (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 |