CN114137616B - 一种利用量子力学原理进行储层含气性检测的方法 - Google Patents
一种利用量子力学原理进行储层含气性检测的方法 Download PDFInfo
- Publication number
- CN114137616B CN114137616B CN202111463516.6A CN202111463516A CN114137616B CN 114137616 B CN114137616 B CN 114137616B CN 202111463516 A CN202111463516 A CN 202111463516A CN 114137616 B CN114137616 B CN 114137616B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- gas
- wave function
- quantum
- seismic
- 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 36
- 230000005610 quantum mechanics Effects 0.000 title claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 230000005428 wave function Effects 0.000 claims abstract description 22
- 238000012512 characterization method Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 7
- 238000005381 potential energy Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 15
- 238000003672 processing method Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 125000001183 hydrocarbyl group Chemical group 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
本发明属于油气勘探地球物理处理方法领域。本发明公开了一种利用量子力学原理进行储层含气性检测的方法。该方法通过将地震信号进行量子力学特征表征，获得适当普朗克常数下的波函数，计算每个波函数的高频衰减系数，在势能‑波函数域估计地震子波穿过含气层导致的高频信号的衰减现象，进而实现储层预测并给出含气性解释结果。本发明提供了一种利用量子力学原理进行储层含气性检测的方法，将传统的各种储层含气性检测技术拓展到了量子领域，首次利用势能‑波函数域内地震信号的量子力学表征波函数进行储层含气性检测，避免了传统储层含气性检测方法中的各种假设前提问题，为传统储层含气性检测技术提供了一种新的解决技术。
Description
技术领域
本发明涉及油气勘探地球物理处理方法领域，具体涉及一种利用量子力学原理进行储层含气性检测的方法。
背景技术
储层气藏勘探的关键是寻找优质储层并进行储层信息提取和流体识别。目前储层信息提取算法，尤其是烃类检测技术，发端于20世纪70年代出现的“亮点”技术。“亮点”技术在一些条件简单的地区取得了成功，但该类方法存在更严重的多解性。随着勘探新思想、新方法和新技术的发展，已出现了很多的储层信息提取方法，目前在实践应用中应用的储层信息提取方法技术概括起来主要有：AVO异常检测技术，叠前弹性波反演和地震属性衰减技术等方法。AVO异常检测技术，是利用地震纵波反射振幅随炮检距（或入射角）的变化特征进行地下地层岩性分析和油气识别的一种方法，其理论基础是Zoeppritz方程。很多学者根据不同的假设，通过一定近似相继推导出多个 Zoeppritz 方程的近似公式并将它们应用于岩性分析和油气识别。基于AVO方法的储层信息提取的主要难点在于预测结果受地震资料品质影响较大；地震道集数据受观测系统和环境噪声及压制噪声方法的影响，一般小偏移距和大偏移距处数据振幅发生畸变，影响AVO分析效果；已有叠前道集数据大都表现为中、小偏移距振幅强而远偏移距弱，对储层信息提取和油气检测有误导或错误指示。叠前弹性波反演技术产生于20 世纪 90 年代末。综合了 AVO 反演和叠后波阻抗反演二者的优点，利用多个部分角度叠加数据的弹性波阻抗反演方法成为地球物理工作者研究的热点。但是叠前反演面临的问题是横波数据少，成本高，对数据品质要求高。地震属性衰减技术主要利用时频分析方法进行衰减属性提取。近年来通过采用小波变换、希尔伯特黄变换等信号分析技术，增加了时频分析的分辨率，并被用于频谱参数的计算和含油气检测。利用频率信息预测含油气主要是根据地震子波穿过油气层所导致高频信号的衰减现象。该类油气检测方法受地层结构、处理精度、时频分析方法等多种因素制约，针对不同地区进行油气检测的可靠性极不稳定。
量子理论是物理世界的数学模型。量子计算利用量子力学原理进行信息处理，作为一种新型的计算理论模型，正逐渐成为信息研究领域的关注焦点。利用量子机制深度挖掘地震信号深层次非线性关系并提高运算速度，进行储层含气性检测已成为突破传统地震勘探处理解释方法的一条重要的途径。
本发明的目的在于解决上述现有技术中存在的难题，提供一种新的利用量子力学原理进行储层含气性检测的高精度方法，采用的是通过将地震信号进行量子力学特征表征，在势能-波函数域估计地震子波穿过含气层导致的高频信号的衰减现象，进而实现储层预测并给出含气性解释结果。
发明内容
一种利用量子力学原理进行储层含气性检测的方法，包括以下步骤：
（1）对于叠后地震数据，综合利用地质、测井以及合成地震记录等资料准确标定层位，确定研究的目标区。
（2）对于目标区的地震数据，确定输入势能表达式，构造哈密顿矩阵，解算量子力学的薛定谔方程，计算哈密顿矩阵的特征值和特征向量，将地震数据在势能-波函数域进行特征表征。
（3）对于目标区的地震数据产生的波函数，逐道逐点结合最小二乘法计算高频衰减系数，对获得的衰减系数剖面进行高斯平滑，归一化，得到最终的含气性检测结果。
本发明的一种利用量子力学原理进行储层含气性检测的技术，其核心问题是将地震数据进行量子力学特征表征，获得适当普朗克常数下的波函数，计算每个波函数的高频衰减系数，利用波函数的高频衰减系数实现储层含气层检测。
本发明的具体实现原理如下：
1. 对于叠后地震数据，综合利用地质、测井以及合成地震记录等资料准确标定层位，确定研究的目标区。
2. 对于目标区的地震数据，确定输入势能表达式，构造哈密顿矩阵，解算量子力学的薛定谔方程，计算哈密顿矩阵的特征值和特征向量，将地震数据在势能-波函数域进行特征表征。
假设势能由矢量形式的信号样本/>表示，即/>。其中，/>为空间位置，/>大小为N，波函数/>属于具有可积功能的希尔伯特空间，势能/>处具有能量/>的平稳量子颗粒的概率满足下式：
。 （1）
量子力学的薛定谔方程可以表示为：
， （2）
其中，哈密顿算子，/>是普朗克常数，/>是量子粒子的质量。/>是梯度算子。
对于矢量信号，其数值微分和数值拉普拉斯变换是
， （3）
从而方程（2）进一步表示为
, （4）
其中，。
整理得：
. （5）
从而有:
. （6）
其中，表示哈密顿矩阵中第/>个元素。
这里，对于储层含气性检测，希尔伯特空间是在有限数量的点中被离散化的。从而，有限矩阵运算符使等式（2）的解等同于对角度化哈密顿矩阵。
构造哈密顿矩阵满足下式：
， （7）
其中，是一个/>矩阵。
计算哈密顿矩阵的特征值和特征向量，记/>，其中，/>为哈密顿矩阵的特征向量，/>为哈密顿矩阵的特征值，/>表示求满秩矩阵特征值和特征向量运算。
地震信号的量子力学特征表征体现在上，即势能-波函数域分解产生的波函数上，并且每个波函数与相应的特征值即能量相关联。
3. 对于目标区的地震数据产生的波函数，逐道逐点结合最小二乘法计算高频衰减系数，对获得的衰减系数剖面进行高斯平滑，归一化，得到最终的含气性检测结果。
对每一道地震信号产生的N条波函数，逐条利用最小二乘法进行高频段拟合，即计算下式
, （8）
其中，为一次近似曲线在点/>处的偏差。拟合范围为选取从波函数最大值开始到大于波函数最大值点的第一个局域最小值处。
通过式（7）求得近似曲线的拟合系数，记为/>，其中/>分别为常数项和一次项拟合系数。从而，对于第/>条波函数/>，高频衰减系数为：
。 （9）
对目标区整个地震数据，有高频衰减系数矩阵：
。 （10）
其中，为地震道道数。
对高频衰减系数矩阵进行高斯平滑，在进行归一化，得到最终的含气性检测剖面。
本发明的一种利用量子力学原理进行储层含气性检测的方法，具有如下特点，主要表现为：
（1）提供了一种地震信号量子力学特征表征方式，利用地震信号量子力学特征推导了储层含气性检测的算法，更适合非线性非平稳地震信号的处理，保证了计算结果的精确性。
（2）将传统的各种储层含气性检测技术拓展到了量子领域，首次利用势能-波函数域内地震信号的量子力学表征波函数进行储层含气性检测。
（3）避免了传统储层含气性检测方法中的各种假设前提问题，为传统储层含气性检测技术提供了一种新的解决技术。
（4）算法运行速度快，适合大批量地震信号处理。
附图说明
图1 为算法流程图。
图2 为依据四川盆地某碳酸盐岩储层测井数据和地震数据建立的模型参数。
图3 为根据模型参数建立的地质模型。
图4 为地质模型的地震响应。
图5 为过含气层的第250道地震道。
图6 为第250道地震道对应的量子域波函数。
图7 为地震响应剖面利用本技术检测的含气性剖面。
具体实施方式
(1) 图1为算法流程图。
(2) 图 2为依据四川盆地某碳酸盐岩储层测井数据和地震数据建立的模型参数。其中，为速度，/>为密度，/>为弥散系数，/>为粘滞系数，Q为品质因子。
(3) 图3为根据模型参数建立的地质模型。采样频率1000Hz。子波频率25Hz。层④为含气层，紧邻的层③为干层。含气层厚度为15m。
(4) 图4为地质模型的地震响应。
(5) 过含气层的第250道地震道。
(6) 为第250道地震道对应的量子域波函数，可以看到，第250道地震道，对应产生320条波函数。
(7) 图5 为地震响应剖面利用本技术检测的含气性剖面。从图中可以看到，本技术很好的检测到了含气层。
Claims (4)
1.一种利用量子力学原理进行储层含气性检测的方法，其特征在于采用以下步骤：
(1)对于叠后地震数据，综合利用地质、测井以及合成地震记录准确标定层位，确定研究的目标区；
(2)对于目标区的地震数据，确定输入势能表达式，构造哈密顿矩阵，解算量子力学的薛定谔方程，计算哈密顿矩阵的特征值和特征向量，将地震数据在势能-波函数域进行特征表征；其中，输入势能V(y)表达式为V(y)＝x(y)ψ(y) (1)
其中，x为矢量形式的信号样本，y为空间位置，x大小为N；ψ(y)为波函数;哈密顿矩阵H具有如下形式：
其中，H是一个N×N矩阵，h是普朗克常数，m是量子粒子的质量；
计算哈密顿矩阵H的特征值和特征向量，记[vectP,valP]＝eig(H)，其中，vectP为哈密顿矩阵的特征向量，valP为哈密顿矩阵的特征值，eig表示求满秩矩阵特征值和特征向量运算；则地震信号的量子力学特征表征体现在vectP上，即势能-波函数域分解产生的波函数上，并且每个波函数与相应的特征值即能量相关联;
(3)对于目标区的地震数据产生的波函数，逐道逐点结合最小二乘法计算高频衰减系数，对获得的衰减系数剖面进行高斯平滑，归一化，得到最终的含气性检测结果。
2.根据权利要求1所述的一种利用量子力学原理进行储层含气性检测的方法，其特征在于：步骤1中准确标定目标区。
3.根据权利要求2所述的一种利用量子力学原理进行储层含气性检测的方法，其特征在于：确定输入势能表达式，选择合适的普朗克常数构造哈密顿矩阵，计算哈密顿矩阵的特征值和特征向量，提取特征向量作为地震数据在势能-波函数域的特征表征函数。
4.根据权利要求3所述的一种利用量子力学原理进行储层含气性检测的方法，其特征在于：利用产生的波函数，逐道逐点在拟合区间结合最小二乘法计算高频衰减系数，对获得的衰减系数剖面进行高斯平滑，归一化，得到最终的含气性检测结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463516.6A CN114137616B (zh) | 2021-12-03 | 2021-12-03 | 一种利用量子力学原理进行储层含气性检测的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463516.6A CN114137616B (zh) | 2021-12-03 | 2021-12-03 | 一种利用量子力学原理进行储层含气性检测的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114137616A CN114137616A (zh) | 2022-03-04 |
| CN114137616B true CN114137616B (zh) | 2023-08-11 |
Family
ID=80387407
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111463516.6A Active CN114137616B (zh) | 2021-12-03 | 2021-12-03 | 一种利用量子力学原理进行储层含气性检测的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114137616B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115061189A (zh) * | 2022-06-22 | 2022-09-16 | 北京世纪金道石油技术开发有限公司 | 基于量子测量的地震波采集方法与系统 |
| CN117538346A (zh) * | 2024-01-08 | 2024-02-09 | 四川阳光墨新科技有限公司 | 基于量子探测技术的混凝土结构裂缝探测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4447020A1 (de) * | 1994-12-30 | 1996-07-04 | Rainer Gottlieb Dr Plaga | Verfahren zum Austausch von Informationen und Energie zwischen quantenmechnanischen Parallelwelten |
| US7383235B1 (en) * | 2000-03-09 | 2008-06-03 | Stmicroelectronic S.R.L. | Method and hardware architecture for controlling a process or for processing data based on quantum soft computing |
| CN102893183A (zh) * | 2010-03-05 | 2013-01-23 | 瓦勒吉有限责任公司 | 改进油气田发现与表征的可预测性的有源噪声注入计算 |
| CN112381232A (zh) * | 2020-11-16 | 2021-02-19 | 成都信息工程大学 | 一种量子模糊机器学习对抗攻击模型方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6819628B2 (en) * | 2003-04-07 | 2004-11-16 | Paradigm Geophysical (Luxembourg) S.A.R.L. | Wave migration by a krylov space expansion of the square root exponent operator, for use in seismic imaging |
-
2021
- 2021-12-03 CN CN202111463516.6A patent/CN114137616B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4447020A1 (de) * | 1994-12-30 | 1996-07-04 | Rainer Gottlieb Dr Plaga | Verfahren zum Austausch von Informationen und Energie zwischen quantenmechnanischen Parallelwelten |
| US7383235B1 (en) * | 2000-03-09 | 2008-06-03 | Stmicroelectronic S.R.L. | Method and hardware architecture for controlling a process or for processing data based on quantum soft computing |
| CN102893183A (zh) * | 2010-03-05 | 2013-01-23 | 瓦勒吉有限责任公司 | 改进油气田发现与表征的可预测性的有源噪声注入计算 |
| CN112381232A (zh) * | 2020-11-16 | 2021-02-19 | 成都信息工程大学 | 一种量子模糊机器学习对抗攻击模型方法 |
Non-Patent Citations (1)
| Title |
|---|
| 量子遗传算法及其在地球物理反演中的应用研究;罗红明;中国优秀硕士论文辑;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114137616A (zh) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111208561B (zh) | 基于时变子波与曲波变换约束的地震声波阻抗反演方法 | |
| CN114137616B (zh) | 一种利用量子力学原理进行储层含气性检测的方法 | |
| CN105093294B (zh) | 基于可变模态分解的地震波衰减梯度估计方法 | |
| Yanhu et al. | A method of seismic meme inversion and its application | |
| CN104090302B (zh) | 工区地下介质频率域异常分析的方法 | |
| Iturrarán-Viveros | Smooth regression to estimate effective porosity using seismic attributes | |
| CN109375265B (zh) | 一种基于变相位雷克子波匹配追踪的理想地震谱分解方法 | |
| CN112162316B (zh) | 一种avo波形数据驱动的高分辨率井震融合叠前反演方法 | |
| CN103675911B (zh) | 一种基于纵波和转换横波联合反演截距和梯度的方法 | |
| CN108020863A (zh) | 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法 | |
| CN107678064B (zh) | 一种声波时差实时提取方法 | |
| CN105301644B (zh) | 基于多参数梯度向量和海色矩阵的油气检测方法及装置 | |
| Zhong et al. | Statistical analysis of background noise in seismic prospecting | |
| CN103728662A (zh) | 一种基于地震信号包络峰值的地层介质品质因子估计方法 | |
| CN106772586A (zh) | 一种基于地震信号奇异性的隐蔽性断裂检测方法 | |
| US10132945B2 (en) | Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume | |
| US20110166842A1 (en) | Layer stripping method | |
| CN114089416B (zh) | 一种利用薛定谔方程进行地震波衰减梯度估计的方法 | |
| CN114415234B (zh) | 基于主动源面波频散和h/v确定浅地表横波速度的方法 | |
| CN113109875B (zh) | 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 | |
| CN109143345B (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
| CN103984013B (zh) | 一种小波域叠前地震道集吸收衰减参数估计算法 | |
| CN114152981B (zh) | 一种势能波函数域地震资料品质因子估计方法 | |
| CN111679326B (zh) | 基于地质探测的瞬态宽带电磁探测信号处理方法 | |
| CN111856564B (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 |