CN112946739B - 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 - Google Patents
裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 Download PDFInfo
- Publication number
- CN112946739B CN112946739B CN202110109196.8A CN202110109196A CN112946739B CN 112946739 B CN112946739 B CN 112946739B CN 202110109196 A CN202110109196 A CN 202110109196A CN 112946739 B CN112946739 B CN 112946739B
- Authority
- CN
- China
- Prior art keywords
- fracture
- carbonate reservoir
- density
- different
- reservoir
- 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
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 682
- 239000011435 rock Substances 0.000 title claims abstract description 290
- 238000000034 method Methods 0.000 title claims abstract description 121
- 238000010276 construction Methods 0.000 title claims abstract description 40
- 239000011148 porous material Substances 0.000 claims abstract description 406
- 239000000945 filler Substances 0.000 claims abstract description 195
- 239000000203 mixture Substances 0.000 claims abstract description 178
- 239000011159 matrix material Substances 0.000 claims abstract description 125
- 238000005266 casting Methods 0.000 claims abstract description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 259
- 238000005260 corrosion Methods 0.000 claims description 66
- 230000007797 corrosion Effects 0.000 claims description 56
- 230000003628 erosive effect Effects 0.000 claims description 29
- 239000012530 fluid Substances 0.000 claims description 21
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 18
- 239000011707 mineral Substances 0.000 claims description 18
- 125000004122 cyclic group Chemical group 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 10
- 230000009977 dual effect Effects 0.000 claims description 7
- 239000010430 carbonatite Substances 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims 1
- 206010017076 Fracture Diseases 0.000 description 497
- 208000010392 Bone Fractures Diseases 0.000 description 460
- 239000007789 gas Substances 0.000 description 90
- 238000010586 diagram Methods 0.000 description 15
- 239000003921 oil Substances 0.000 description 9
- 238000011160 research Methods 0.000 description 8
- 239000000463 material Substances 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 239000011499 joint compound Substances 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000006467 substitution reaction Methods 0.000 description 5
- 150000001875 compounds Chemical class 0.000 description 4
- 238000004090 dissolution Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 3
- 238000002441 X-ray diffraction Methods 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 239000002002 slurry Substances 0.000 description 3
- 238000005056 compaction Methods 0.000 description 2
- 235000009508 confectionery Nutrition 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 1
- 230000002378 acidificating effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000012455 biphasic mixture Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 235000011850 desserts Nutrition 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 229910000037 hydrogen sulfide Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009291 secondary effect Effects 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
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
技术领域
本发明属于地球物理技术领域，特别涉及一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法与基于此的储层参数预测方法。
背景技术
现有的碳酸盐岩岩石物理模型主要针对浅层碳酸盐岩储层，一般含有大量的原生孔隙。而中深层碳酸盐岩储层的孔隙结构与浅层差别较大，其流体与孔隙结构耦合的弹性响应也与中深层差异较大，因此现有的面向中浅层的碳酸盐岩地震岩石物理模板不能很好地反映中深层碳酸盐岩储层的地质与地球物理属性。
碳酸盐岩岩石物理模板可以有效地串联起碳酸盐岩储层岩石物理响应特征与储层地质属性，同时对于定量解释地球物理数据(实验、测井、地震)也具有重要意义。随着对深层超深层碳酸盐岩储层探索的不断深入，基于Xu-White模型建立起的浅层碳酸盐岩岩石物理模型已经不能满足我们对于中深层储层解释的需要。相比于浅层碳酸盐岩储层，中深层致密碳酸盐岩经历过更长时间的沉积压实作用，整体岩层显得十分致密，孔隙度较低，并且受次生作用影响显著，后天发育大量裂隙、溶蚀孔，形成了独特的孔隙系统。此外，中深层储层不同的孔隙充填物和流体属性也是影响中深层碳酸盐岩储层地球物理响应特征的重要因素。因此，考虑符合中深层碳酸盐岩储层孔隙结构非均质性和流体属性的地震岩石物理模型对于地震数据和测井数据的定量地震解释都具有重要的应用意义。
发明内容
本发明的目的在于提供一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法。该方法建立的岩石物理模板能够实现非均质裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层弹性特征刻画，对中深层复杂非均质碳酸盐岩储层的速度建模、流体识别、储层“甜点”的刻画具有重要的指导意义。
本发明的另一目的在于提供一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法，该方法实现了对酸盐岩储层的孔隙填充物组成和/或裂隙密度进行大致预测。
为了实现上述目的，本发明第一方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法，其中，该方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
本发明第二方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法，其中，该方法包括：
采用上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
获取目标储层的纵横波参数比数据和纵波阻抗数据；
利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
本发明第三方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统，其中，该系统包括：
第一参数确定模块：用于获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙(即溶蚀孔)的纵横比；
第二参数确定模块：用于获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
第三参数确定模块：用于获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
弹性参数获取模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
物理模板构建模块：用于将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
本发明第四方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测系统，其中，该系统包括：
上述碳酸盐岩储层地震岩石物理模板构建系统：用于建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
目标储层弹性参数获取模块：用于获取目标储层的纵横波参数比数据和纵波阻抗数据；
预测模块：用于利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
本发明第五方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建装置，包括处理器及存储器；其中，
存储器，用于存放计算机程序；
处理器，用于执行存储器上所存放的程序时，实现上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法的步骤。
本发明第六方面实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测装置，包括处理器及存储器；其中，
存储器，用于存放计算机程序；
处理器，用于执行存储器上所存放的程序时，实现上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法的步骤。
本发明第七方面实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法的步骤。
本发明第八方面实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法的步骤。
本发明提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法构建得到的物理模板可以用于裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层测井弹性数据和叠前地震反演结果的定量解释，可以用于指示裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层的孔隙充填物属性和裂隙密度，对储层甜点识别和开发优选具有重要应用价值。
附图说明
图1为本发明一实施例提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法的流程示意图。
图2为本发明一实施例提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法的流程示意图。
图3为本发明一实施例提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统的结构示意图。
图4为本发明一实施例提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测系统的结构示意图。
图5为本发明一实施例提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建装置的结构示意图。
图6为裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层双孔系统的典型铸体切片图。
图7为本发明实施例1提供的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法流程示意图。
图8为本发明实施例1建立的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明的保护范围。
下面参考本发明的若干代表性实施方式，详细阐述本发明的原理和精神。
参见图1，为了实现上述目的，本发明一实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法，其中，该方法包括：
步骤S1：获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
步骤S2：获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
步骤S3：获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
步骤S4：在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
步骤S5：将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
其中，将裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层看做裂隙-溶蚀孔双孔系统，换言之，将碳酸盐岩储层的储集空间分成裂隙、孔隙(溶蚀孔)。对于中深层致密碳酸盐岩储层来说，在长期的压实胶结过程中，大多数的原生孔隙已经消失，受多个地质时期的地层抬升、风化剥蚀等次生作用影响，裂隙和溶蚀孔发育尤为显著。图6展示了根据致密碳酸盐岩储层典型的孔隙系统分布特征抽象出的：裂隙(微裂缝)-溶蚀孔双孔系统模型。由此，中深层复杂非均质碳酸盐储层可抽象成裂隙-溶蚀孔双孔系统。
其中，碳酸盐岩储层中的孔隙填充物优选为泥水、油水和/或油气；此时，不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，和/或，不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，和/或，不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同泥水组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度确定的条件下，获取不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同油水组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度确定的条件下，获取不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同油气组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物泥、水和油的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥、水和油的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同泥水组成比例、不同裂隙密度和不同油水组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油气组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或气)的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同泥水组成比例、不同裂隙密度和不同油气组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物水、油和气的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物水、油和气的体积模量、剪切模量与密度确定的条件下，获取不同油水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油气组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或气)的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同油水组成比例、不同裂隙密度和不同油气组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
在一实施方式中，裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油气组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或气)的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同泥水组成比例、不同裂隙密度以及不同油气组成比例、不同裂隙密度以及不同油水组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
进一步地，获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取泥水比例范围和裂隙密度范围；
获取泥水比例迭代步长以及裂隙密度迭代步长；
在泥水比例范围和裂隙密度范围内，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代，获取各泥水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代过程中，每当泥水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
例如，获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
1)令储层含水饱和度为100％即含泥饱和度为0％(或令储层含水饱和度为0％即含泥饱和度为100％)；
2)令裂隙密度为最大值(或最小值)；
3)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度确定的条件下，确定在此时的含水饱和度即含泥饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
4)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；
5)重复步骤3)-步骤4)直至新的裂隙密度不大于最小值(或不小于最大值)；
6)给定一个含水饱和度的迭代步长即含泥饱和度的迭代步长，将含水饱和度减去(或加上)含水饱和度的迭代步长作为新的含水饱和度即将含泥饱和度加上(或减去)含泥饱和度的迭代步长作为新的含泥饱和度；
7)重复步骤2)-步骤6)直至新的含水饱不大于零即含泥饱和度不小于100％(或含水饱不小于100％即含泥饱和度不大于零)，从而完成确定不同含水饱和度即含泥饱和度、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同泥水组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
再例如，获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
a)令储层含水饱和度为100％即含泥饱和度为0％(或令储层含水饱和度为0％即含泥饱和度为100％)；
b)令裂隙密度为最大值(或最小值)；
c)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度确定的条件下，确定在此时的含水饱和度即含泥饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
d)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；若新的裂隙密度大于等于最小值(或小于等于最大值)，则执行步骤e)；若新的裂隙密度小于最小值(或小于等于最大值)，则执行步骤f)；
e)重复步骤c)确定在此时的含水饱和度即含泥饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗以及步骤d)；
f)给定一个含水饱和度的迭代步长即含泥饱和度的迭代步长，将含水饱和度减去(或加上)含水饱和度的迭代步长作为新的含水饱和度即将含泥饱和度加上(或减去)含泥饱和度的迭代步长作为新的含泥饱和度；若新的含水饱和度大于等于零即含泥饱和度小于等于100％(或含水饱和度小于等于100％即含泥饱和度大于等于零)，则执行步骤g)；若新的含水饱和度小于零即含泥饱和度大于100％(或含水饱和度大于100％即含泥饱和度小于零)，则完成确定不同含水饱和度即含泥饱和度、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同泥水组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
g)重复步骤b)-步骤f)。
进一步地，获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油水比例范围和裂隙密度范围；
获取油水比例迭代步长以及裂隙密度迭代步长；
在油水比例范围和裂隙密度范围内，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代，获取各油水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代过程中，每当油水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
例如，获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
1)令储层含水饱和度为100％即含油饱和度为0％(或令储层含水饱和度为0％即含油饱和度为100％)；
2)令裂隙密度为最大值(或最小值)；
3)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度确定的条件下，确定在此时的含水饱和度即含油饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
4)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；
5)重复步骤3)-步骤4)直至新的裂隙密度不大于最小值(或不小于最大值)；
6)给定一个含水饱和度的迭代步长即含油饱和度的迭代步长，将含水饱和度减去(或加上)含水饱和度的迭代步长作为新的含水饱和度即含油饱和度加上(或减去)含油饱和度的迭代步长作为新的含油饱和度；
7)重复步骤2)-步骤6)直至新的含水饱不大于零即含油饱和度不小于100％(或含水饱不小于100％即含油饱和度不大于零)，从而完成确定不同含水饱和度即含油饱和度、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同油水组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
再例如，获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
a)令储层含水饱和度为100％即含油饱和度为0％(或令储层含水饱和度为0％即含油饱和度为100％)；
b)令裂隙密度为最大值(或最小值)；
c)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度确定的条件下，确定在此时的含水饱和度即含油饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
d)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；若新的裂隙密度大于等于最小值(或小于等于最大值)，则执行步骤e)；若新的裂隙密度小于最小值(或小于等于最大值)，则执行步骤f)；
e)重复步骤c)确定在此时的含水饱和度即含油饱和度、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗以及步骤d)；
f)给定一个含水饱和度的迭代步长即含油饱和度的迭代步长，将含水饱和度减去(或加上)含水饱和度的迭代步长作为新的含水饱和度即将含油饱和度加上(或减去)含油饱和度的迭代步长作为新的含油饱和度；若新的含水饱和度大于等于零即含油饱和度小于等于100％(或含水饱和度小于等于100％即含油饱和度大于等于零)，则执行步骤g)；若新的含水饱和度小于零即含油饱和度大于100％(或含水饱和度大于100％即含油饱和度小于零)，则完成确定不同含水饱和度即含油饱和度、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同油水组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
g)重复步骤b)-步骤f)。
进一步地，获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油气比例范围和裂隙密度范围；
获取油气比例迭代步长以及裂隙密度迭代步长；
在油气比例范围和裂隙密度范围内，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代，获取各油气比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代过程中，每当油气比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度。
例如，获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
1)令储层气-油比为最小值(或令储层气-油比为最大值)；
2)令裂隙密度为最大值(或最小值)；
3)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度确定的条件下，确定在此时的气-油比、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
4)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；
5)重复步骤3)-步骤4)直至新的裂隙密度不大于最小值(或不小于最大值)；
6)给定一个气-油比的迭代步长，将气-油比加上(或减去)气-油比的迭代步长作为新的气-油比作为新的气-油比；
7)重复步骤2)-步骤6)直至新的气-油比不小于最大值(或不大于最小值)，从而完成确定不同气-油比、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同油气组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
再例如，获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
a)令储层气-油比为最小值(或令储层气-油比为最大值)；
b)令裂隙密度为最大值(或最小值)；
c)在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度确定的条件下，确定在此时的气-油比、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
d)给定一个裂隙密度的迭代步长，将裂隙密度减去(或加上)裂隙密度的迭代步长作为新的裂隙密度；若新的裂隙密度大于等于最小值(或小于等于最大值)，则执行步骤e)；若新的裂隙密度小于最小值(或小于等于最大值)，则执行步骤f)；
e)重复步骤c)确定在此时的气-油比、裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗以及步骤d)；
f)给定一个气-油比的迭代步长，将气-油比加上(或减去)气-油比的迭代步长作为新的气-油比作为新的气-油比；若新的气-油比小于等于最大值(或大于等于最小值)，则执行步骤g)；若新的气-油比大于最大值(或小于最小值)，则完成确定不同气-油比、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗即不同油气组成比例、不同裂隙密度条件下的碳酸盐岩储层的纵横波参数比和纵波阻抗；
g)重复步骤b)-步骤f)。
在一实施方式中，裂隙密度范围为0-0.3；例如，裂隙密度的最小值为0、最大值为0.3。
在一实施方式中，气-油比的范围为0-2000；例如，气-油比的最小值为0、最大值为2000。
其中，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗优选通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
进一步地，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中各类孔隙填充物的体积模量与剪切模量，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中各类孔隙填充物的密度，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
例如，获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度，确定不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
进一步地，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物泥和水的体积模量、剪切模量与密度，确定不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物泥和水的体积模量与剪切模量，确定不同泥水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同泥水组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中泥和水的密度，确定不同泥水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
例如，获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度，确定不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
进一步地，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物油和水的体积模量、剪切模量与密度，确定不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物油和水的体积模量与剪切模量，确定不同油水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同油水组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中油和水的密度，确定不同油水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
例如，获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度，确定不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
进一步地，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中孔隙填充物油和气的体积模量、剪切模量与密度，确定不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物油和气的体积模量与剪切模量，确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同油气组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中油和气的密度，确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
其中，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量优选包括：
基于裂隙的纵横比确定不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度；
基于碳酸盐岩储层溶蚀孔孔隙度以及不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度；
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，结合不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度、不同裂隙密度条件下的碳酸盐岩储层总孔隙度，利用微分等效介质模型，确定不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量。
在一实施方式中，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗通过下述公式进行：
式中，α为纵横波参数比；Vp为纵波速度；Vs为横波速度；K为体积模量；μ为剪切模量；ρ为密度；Ip为纵波阻抗；Is为横波阻抗。
在一实施方式中，确定不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度利用碳酸盐岩储层裂隙的纵横比和裂隙密度确定得到，具体可以通过下述公式进行：
式中，
为裂隙密度为cd时的碳酸盐岩储层裂隙孔隙度；cd为裂隙密度(无量纲)；αcrack为裂隙的纵横比。
在一实施方式中，不同裂隙密度条件下的碳酸盐岩储层总孔隙度通过下述公式确定的到：
式中，
为特定裂隙密度时的碳酸盐岩储层总孔隙度；
为碳酸盐岩储层溶蚀孔孔隙度；
为特定裂隙密度时的碳酸盐岩储层裂隙孔隙度；
例如，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度利用碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比和裂隙密度确定得到，具体通过下述公式进行：
式中，
为裂隙密度为cd时的碳酸盐岩储层总孔隙度；
为碳酸盐岩储层溶蚀孔孔隙度；
为裂隙密度为cd时的碳酸盐岩储层裂隙孔隙度；cd为裂隙密度；αcrack为裂隙的纵横比。
在一实施方式中，所述微分等效介质模型为DEM微分等效介质模型，能够刻画含有不同矿孔隙类型的碳酸盐岩的等效模量，通过往固体矿物中相中逐渐加入包含物相来模拟双相混合物；
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，结合不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度、不同裂隙密度条件下的碳酸盐岩储层总孔隙度，利用DEM微分等效介质模型，确定不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，采用常规方法进行即可；
具体而言，DEM微分等效介质模型为：
初始条件为：K*(0)＝K1，μ*(0)＝μ1
式中，K*代表等效体积模量(计算时对应干岩石的等效体积模量)；μ*代表等效剪切模量(计算时对应干岩石的等效剪切模量)；K1为初始主相相1材料的体积模量；μ1为初始主相相1材料的剪切模量；K2为逐渐加入的包含物相2的体积模量；μ2为逐渐加入的包含物相2的剪切模量；y为相2的含量即孔隙度；P，Q为几何因数；(*2)指此几何因数是针对具有等效模量K*和μ*的背景介质中的包含物材料相2；
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，结合不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度、不同裂隙密度条件下的碳酸盐岩储层总孔隙度，利用微分等效介质模型，确定不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量通过下述公式确定：
式中，Kdry为裂隙孔隙度
对应的裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量；μdry为裂隙孔隙度
对应的裂隙密度条件下的碳酸盐岩储层干岩石的等效剪切模量；K1为储层岩石基质的体积模量；μ1为储层岩石基质的剪切模量；αcrack为裂隙的纵横比；αpore为孔隙的纵横比；Kpore为空气的弹性模量；μpore为空气的的剪切模量；
为裂隙孔隙度
对应的裂隙密度条件下的碳酸盐岩储层总孔隙度；
为碳酸盐岩储层溶蚀孔孔隙度；
为裂隙孔隙度；其中Kpore，μpore可以取值0。
在一实施方式中，所述基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中各类孔隙填充物的体积模量与剪切模量，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量利用Gassmann流体替代公式进行；将干岩石模型模量替换为特定孔隙填充物组成比例条件下的岩石物理模型模量；
当孔隙充填物为泥水时，优选通过下述公式确定不同泥水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量：
其中，Kif＝(Kwater·(1-smud)+Kmud·smud)
μif＝μmud·smud
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；μ1为储层岩石基质的剪切模量，Gpa；Kwater为孔隙充填物水的体积模量，Gpa；Kmud为孔隙充填物泥的体积模量，Gpa；μmud为孔隙充填物泥的剪切模量，Gpa；smud为储层含泥饱和度，％；Kif为孔隙充填物泥水的等效体积模量，Gpa；μif为孔隙充填物泥水的等效剪切模量，Gpa；
当孔隙充填物为油水时，优选通过下述公式确定不同油水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量：
μsat＝μdry
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；Koil为孔隙充填物油的体积模量，Gpa；
当孔隙充填物为油气时，优选通过下述公式确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量：
μsat＝μdry
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；Koil(GOR)为气-油比GOR的孔隙充填物油气的体积模量，Gpa；
其中，不同气-油比GOR的孔隙充填物油气的弹性模量可以由实验时测量数据标定获得(在实验室测量标定时所使用的油气选择与孔隙充填物油气弹性模量及密度相同的模拟油气)。
在一实施方式中，基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中油和气的密度，确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度，采用常规方法进行即可；
例如，当孔隙充填物为泥水时，优选通过下述公式确定不同泥水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度：
式中，ρsat为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρwater为孔隙充填物水的密度，g/cm3；ρmud为孔隙充填物泥的密度，g/cm3；smud为储层含泥饱和度，％；
当孔隙充填物为油水时，优选通过下述公式确定不同油水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度：
式中，ρmodel为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρwater为孔隙充填物水的密度，g/cm3；ρoil为孔隙充填物油的密度，g/cm3；swater为储层含水饱和度，％；
当孔隙充填物为油气时，优选通过下述公式确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度：
式中，ρd为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρoil(GOR)为气-油比GOR的孔隙充填物油气的密度，g/cm3；
其中，不同气-油比GOR的孔隙充填物油气的密度可以由实验时测量数据标定获得(在实验室测量标定时所使用的油气选择与孔隙充填物油气弹性模量及密度相同的模拟油气)。
其中，获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度优选包括：
获取碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度；
基于碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
进一步地，基于碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度利用Voigt-Reuss-Hill平均理论进行，具体可以通过下述公式确定：
其中，fi为第i种矿物的体积分数，％；Mi为第i种矿物的体积模量或剪切模量；M为碳酸盐岩储层岩石基质的体积模量或剪切模量，Gpa；ρi为第i种矿物的密度，g/cm3；ρ为碳酸盐岩储层岩石基质的密度，g/cm3。
在一实施方式中，裂隙的纵横比为0.001-0.01。
在一实施方式中，孔隙(溶蚀孔)的纵横比为0.4-0.7。
其中，将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板优选包括：
将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；
拟合各孔隙填充物组成比例下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同孔隙填充物组成比例的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；从而完成碳酸盐岩储层地震岩石物理模板建立。
在一实施方式中，将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板包括：
将不同泥水组成比例(例如以含泥饱和度计)、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；拟合各泥水组成比例(例如以含泥饱和度计)下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同泥水组成比例(例如以含泥饱和度计)的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；
和/或，
将不同油水组成比例(例如以含油饱和度计)、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；拟合各油水组成比例(例如以含油饱和度计)下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同油水组成比例(例如以含油饱和度计)的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；
和/或，
将不同油气组成比例(例如以气-油比计)、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；拟合各油气组成比例(例如以气-油比计)下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同油气组成比例(例如以气-油比计)的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；
从而完成裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板建立。
参见图2，本发明一实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法，其中，该方法包括：
步骤S21：采用上述的方法实施例所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
步骤S22：获取目标储层的纵横波参数比数据和纵波阻抗数据；
步骤S23：利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
在一实施方式中，将实际测井数据或叠前地震反演结果确定的纵横波参数比数据和纵波阻抗数据投射到建立的碳酸盐岩储层地震岩石物理模板上，即可以得到各个数据点所对应的孔隙填充物属性(充填物类别和流体饱和度)以及裂隙发育程度(裂隙密度)。
实施例1
该实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法。
该方法用于针对某研究区储层建立适用于该研究区的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板，其流程示意图如图7所示，通过对该研究区的测井、实验室数据分析，对其岩芯切片扫面观察，可以得到基质矿物以及孔隙填充物(流体、泥质等)模量、矿物百分比、孔隙度、油气比(GOR)、溶孔与裂隙的纵横比等参数，利用VRH平均理论估算得到基质的模量，接着利用DEM理论建立“双孔”系统的干岩石地震岩石物理模板，最后用Gassmann流体或者固体替代方程，获得不同流体填充物的“双孔”系统中深层碳酸盐岩地震岩石物理综合模板。将研究区测井或者叠前地震数据投射到本实施例建立的岩石物理模板上，可以对储层裂隙发育程度和孔隙填充物性质进行进一步分析。
具体而言，该方法包括：
1、获取研究区参数
包括：获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；获取研究区碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；获取研究区碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度；
具体而言：通过目标岩心薄片扫描获取裂隙的纵横比和孔隙的纵横比；通过实验时测量、XRD与测井数据分析结合目标岩心薄片扫描获取孔隙填充物的体积模量、剪切模量与密度并获取溶蚀孔孔隙度孔隙度，在本实施例中获取泥质、水和不同气-油比GOR时油气混合流体的体积模量、剪切模量与密度；通过实验时测量、XRD与测井数据分析获取岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，然后利用Voigt-Reuss-Hill平均理论确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
研究区域的岩芯铸体切片如图6所示，从中不难看出，在中深层碳酸盐岩储层为复杂多孔裂隙介质，裂隙、溶蚀孔发育尤其明显，从中抽象出来其对应的孔隙系统：裂隙-孔隙的“双孔”系统；该碳酸盐岩储层含有大量的酸性腐蚀气体硫化氢，发育较多溶蚀孔，而很多微小溶蚀缝进一步联通了这些溶蚀孔，在整体超压情况下保持了这样的孔隙系统，形成了“双孔”系统；获取得到的研究区中深层碳酸盐岩基本参数，溶蚀孔孔隙度
溶蚀孔纵横比αpore＝0.4，裂隙纵横比αcrack＝0.001，其他如表1、表2所示。
表1
| 体积模量(GPa) | 体积模量(GPa) | 密度(g/cc) | |
| 方解石 | 76.8 | 32 | 2.71 |
| 泥质 | 15 | 4 | 2.79 |
| 水 | 2.25 | 0 | 1 |
| 油(GOR＝0) | 0.9629 | 0 | 0.6678 |
表2
2、获取考虑孔隙流体和填充物的“双孔隙系统”模板的弹性模量
包括：在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙(溶蚀孔)的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中孔隙填充物泥、水、油和气的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油水组成比例、不同裂隙密度时(假定此时孔隙填充物仅为泥和/或水)的碳酸盐岩储层的纵横波参数比和纵波阻抗，获取不同油气组成比例、不同裂隙密度时(假定此时孔隙填充物仅为油和/或气)的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
具体而言：
2.1确定“双孔隙系统”干岩石模板的弹性模量
基于裂隙的纵横比确定不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度；基于碳酸盐岩储层溶蚀孔孔隙度以及不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度；基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，结合不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度、不同裂隙密度条件下的碳酸盐岩储层总孔隙度，利用微分等效介质模型，确定不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
其中，不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度通过下述公式确定：
式中，
为裂隙密度为cd时的碳酸盐岩储层总孔隙度，％；
为碳酸盐岩储层溶蚀孔孔隙度，％；
为裂隙密度为cd时的碳酸盐岩储层裂隙孔隙度，％；cd为裂隙密度；αcrack为裂隙的纵横比；
2.2确定泥水填充模板的“双孔隙系统”的弹性模量
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物泥和水的体积模量与剪切模量，利用Gassmann流体替代公式确定不同含泥饱和度条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同泥水组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；具体通过下述公式进行：
其中，Kif＝(Kwater·(1-smud)+Kmud·smud)
μif＝μmud·smud
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；μ1为储层岩石基质的剪切模量，Gpa；Kwater为孔隙充填物水的体积模量，Gpa；Kmud为孔隙充填物泥的体积模量，Gpa；μmud为孔隙充填物泥的剪切模量，Gpa；smud为储层含泥饱和度，％；Kif为孔隙充填物泥水的等效体积模量，Gpa；μif为孔隙充填物泥水的等效剪切模量，Gpa；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中泥和水的密度，确定不同泥水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度；具体通过下述公式进行：
式中，ρsat为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρwater为孔隙充填物水的密度，g/cm3；ρmud为孔隙充填物泥的密度，g/cm3；smud为储层含泥饱和度，％；
基于不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
2.3确定油水填充模板的“双孔隙系统”的弹性模量
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物油和水的体积模量与剪切模量，利用Gassmann流体替代公式确定不同含税饱和度条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同油水组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；具体通过下述公式进行：
μsat＝μdry
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；Koil为孔隙充填物油的体积模量，Gpa；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中油和水的密度，确定不同油水组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度；具体通过下述公式进行：
式中，ρmodel为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρwater为孔隙充填物水的密度，g/cm3；ρoil为孔隙充填物油的密度，g/cm3；swater为储层含水饱和度，％；
基于不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
2.4确定油气填充模板的“双孔隙系统”的弹性模量
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中孔隙填充物油和气的体积模量与剪切模量，利用Gassmann流体替代公式确定不同气-油比条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同油气组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；具体通过下述公式进行：
μsat＝μdry
式中，Ksat为碳酸盐岩储层的等效体积模量，Gpa；μsat为碳酸盐岩储层的等效剪切模量，Gpa；
为碳酸盐岩储层总孔隙度，％；Kdry为碳酸盐岩储层干岩石的等效体积模量，Gpa；μdry为碳酸盐岩储层干岩石的等效剪切模量，Gpa；K1为储层岩石基质的体积模量，Gpa；Koil(GOR)为气-油比GOR的孔隙充填物油气的体积模量，Gpa；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中油和气的密度，确定不同油气组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度；具体通过下述公式进行：
式中，ρd为碳酸盐岩储层的等效密度，g/cm3；
为碳酸盐岩储层总孔隙度，％；ρ1为储层岩石基质的密度，g/cm3；ρoil(GOR)为气-油比GOR的孔隙充填物油气的密度，g/cm3；
基于不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，裂隙密度cd范围为0-0.15；气-油比的范围为0-2000；含泥饱和度范围为0-1；含水饱和度范围为0-1。
3、构建“双孔”系统中深层碳酸盐岩储层地震岩石物理综合模板
将不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗以及不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同泥水组成比例、不同裂隙密度以及不同油气组成比例、不同裂隙密度以及不同油水组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
结果如图8所示。
将研究区目标储层的测井数据、叠前地震反演结果确定的纵横波参数比数据和纵波阻抗数据投点到本实施例建立的述碳酸盐岩储层地震岩石物理模板上，即可以得到各个数据点所对应的孔隙填充物属性(充填物类别和流体饱和度)以及裂隙发育程度(裂隙密度)，进而推动储层甜点预测、开发优选工作的进行。
本发明实施例还提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统，优选地，该系统用于实现上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法实施例。
图3是根据本发明实施例的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统的结构框图，如图3所示，该系统包括：
第一参数确定模块31：用于获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙(即溶蚀孔)的纵横比；
第二参数确定模块32：用于获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
第三参数确定模块33：用于获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
弹性参数获取模块34：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
物理模板构建模块35：用于将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板。
其中，碳酸盐岩储层中的孔隙填充物为泥水、油水和/或油气时，此时，弹性参数获取模块34优选包括：
第一弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
和/或，
第二弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
和/或，
第三弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗。
其中，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗优选通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗。
其中，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度优选通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中各类孔隙填充物的体积模量与剪切模量，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中各类孔隙填充物的密度，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度；
进一步地，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量优选包括：
基于裂隙的纵横比确定不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度；
基于碳酸盐岩储层溶蚀孔孔隙度以及不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度；
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，结合不同裂隙密度条件下的碳酸盐岩储层裂隙孔隙度、不同裂隙密度条件下的碳酸盐岩储层总孔隙度，利用微分等效介质模型，确定不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量。
其中，第二参数确定模块32优选包括：
矿物参数获取子模块：用于获取碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度；
岩石基质参数确定子模块：用于基于碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度。
其中，物理模板构建模块35优选包括：
投点子模块：用于将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；
拟合子模块：用于拟合各孔隙填充物组成比例下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同孔隙填充物组成比例的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；从而完成碳酸盐岩储层地震岩石物理模板建立。
本发明实施例还提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测系统，优选地，该系统用于实现上述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法实施例。
图4是根据本发明实施例的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测系统的结构框图，如图4所示，该系统包括：
上述实施例的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统41：用于建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
目标储层弹性参数获取模块42：用于获取目标储层的纵横波参数比数据和纵波阻抗数据；
预测模块43：用于利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
图5是根据本发明实施例的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建装置的示意图。图5所示的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建装置为通用数据处理装置，其包含通用的计算机硬件结构，其至少包含处理器1000、存储器1111；所述处理器1000用于执行所述存储器中存储的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建程序，以实现各方法实施例所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法(具体方法参见上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法实施例的描述，在此不再赘述)。
本发明一实施例提供了一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测装置，其包含通用的计算机硬件结构，其至少包含处理器、存储器；所述处理器用于执行所述存储器中存储的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测程序，以实现各方法实施例所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法(具体方法参见上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法实施例的描述，在此不再赘述)。
本发明实施例还提供了一种计算机可读存储介质，所述存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现各方法实施例所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法(具体方法参见上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法实施例的描述，在此不再赘述)。
本发明实施例还提供了一种计算机可读存储介质，所述存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现各方法实施例所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法(具体方法参见上述裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法实施例的描述，在此不再赘述)。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (17)
1.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法，其中，该方法包括：
获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板；
其中，碳酸盐岩储层中的孔隙填充物为泥水、油水和/或油气；此时，不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，和/或，不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，和/或，不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取泥水比例范围和裂隙密度范围；
获取泥水比例迭代步长以及裂隙密度迭代步长；
在泥水比例范围和裂隙密度范围内，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代，获取各泥水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代过程中，每当泥水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油水比例范围和裂隙密度范围；
获取油水比例迭代步长以及裂隙密度迭代步长；
在油水比例范围和裂隙密度范围内，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代，获取各油水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代过程中，每当油水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油气比例范围和裂隙密度范围；
获取油气比例迭代步长以及裂隙密度迭代步长；
在油气比例范围和裂隙密度范围内，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代，获取各油气比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代过程中，每当油气比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度。
2.根据权利要求1所述的构建方法，其中，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗。
3.根据权利要求2所述的构建方法，其中，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中各类孔隙填充物的体积模量与剪切模量，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中各类孔隙填充物的密度，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
4.根据权利要求1所述的构建方法，其中，获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度包括：
获取碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度；
基于碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度。
5.根据权利要求1-4任一项所述的构建方法，其中，将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板包括：
将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；
拟合各孔隙填充物组成比例下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同孔隙填充物组成比例的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；从而完成碳酸盐岩储层地震岩石物理模板建立。
6.根据权利要求1所述的构建方法，其中，
裂隙的纵横比为0.001-0.01；
孔隙的纵横比为0.4-0.7。
7.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法，其中，该方法包括：
采用权利要求1-6任一项所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
获取目标储层的纵横波参数比数据和纵波阻抗数据；
利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
8.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统，其中，该系统包括：
第一参数确定模块：用于获取碳酸盐岩储层溶蚀孔孔隙度，裂隙的纵横比和孔隙的纵横比；
第二参数确定模块：用于获取碳酸盐岩储层岩石基质的体积模量、剪切模量与密度；
第三参数确定模块：用于获取碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度；
弹性参数获取模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；其中，所述纵横波参数比为纵横波阻抗比或者纵横波速度比；
物理模板构建模块：用于将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中，建立不同孔隙填充物组成比例、不同裂隙密度的纵横波参数比值与纵波阻抗值的查询图版即为所述碳酸盐岩储层地震岩石物理模板；
其中，碳酸盐岩储层中的孔隙填充物为泥水、油水和/或油气；此时，弹性参数获取模块包括：
第一弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同泥水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
和/或，
第二弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同油水组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
和/或，
第三弹性参数获取子模块：用于在碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比确定，岩石基质的体积模量、剪切模量与密度确定，碳酸盐岩储层中各孔隙填充物的体积模量、剪切模量与密度确定的条件下，获取不同油气组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取泥水比例范围和裂隙密度范围；
获取泥水比例迭代步长以及裂隙密度迭代步长；
在泥水比例范围和裂隙密度范围内，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代，获取各泥水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同泥水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照泥水比例迭代步长以及裂隙密度迭代步长对泥水比例以及裂隙密度进行循环迭代过程中，每当泥水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油水比例范围和裂隙密度范围；
获取油水比例迭代步长以及裂隙密度迭代步长；
在油水比例范围和裂隙密度范围内，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代，获取各油水比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油水比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油水比例迭代步长以及裂隙密度迭代步长对油水比例以及裂隙密度进行循环迭代过程中，每当油水比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度；
获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗包括：
获取油气比例范围和裂隙密度范围；
获取油气比例迭代步长以及裂隙密度迭代步长；
在油气比例范围和裂隙密度范围内，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代，获取各油气比例、各裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗，从而完成获取不同油气比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗；
其中，按照油气比例迭代步长以及裂隙密度迭代步长对油气比例以及裂隙密度进行循环迭代过程中，每当油气比例迭代一个步长，以裂隙密度迭代步长循环迭代裂隙密度。
9.根据权利要求8所述的系统，其中，获取不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度；
基于不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度，确定不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗。
10.根据权利要求9所述的系统，其中，基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量、剪切模量与密度，和，碳酸盐岩储层中各类孔隙填充物的体积模量、剪切模量与密度，确定不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量、等效剪切模量以及等效密度通过下述方式实现：
基于碳酸盐岩储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比，储层岩石基质的体积模量与剪切模量，确定不同裂隙密度条件下的碳酸盐岩储层总孔隙度、碳酸盐岩储层裂隙孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度以及不同裂隙密度条件下的碳酸盐岩储层干岩石的等效体积模量与等效剪切模量，结合碳酸盐岩储层中各类孔隙填充物的体积模量与剪切模量，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效体积模量和等效剪切模量即不同孔隙填充物组成比例、不同裂隙密度条件下，碳酸盐岩储层的等效体积模量与等效剪切模量；
基于不同裂隙密度条件下的碳酸盐岩储层总孔隙度，结合储层岩石基质密度碳酸盐岩储层中各类孔隙填充物的密度，确定不同孔隙填充物组成比例条件下的不同裂隙密度的碳酸盐岩储层的等效密度。
11.根据权利要求8所述的系统，其中，第二参数确定模块包括：
矿物参数获取子模块：用于获取碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度；
岩石基质参数确定子模块：用于基于碳酸盐岩储层岩石基质中各种类矿物的含量、剪切模量、体积模量和密度，确定碳酸盐岩储层岩石基质的体积模量、剪切模量与密度。
12.根据权利要求8所述的系统，其中，物理模板构建模块包括：
投点子模块：用于将不同孔隙填充物组成比例、不同裂隙密度时的碳酸盐岩储层的纵横波参数比和纵波阻抗投点到坐标轴为纵横波参数比和纵波阻抗的坐标系统中；
拟合子模块：用于拟合各孔隙填充物组成比例下不同裂隙密度的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；和/或，拟合各裂隙密度下不同孔隙填充物组成比例的碳酸盐岩储层的纵横波参数比与纵波阻抗的关系趋势线；从而完成碳酸盐岩储层地震岩石物理模板建立。
13.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测系统，其中，该系统包括：
权利要求8-12任一项所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建系统：用于建立与目标储层溶蚀孔孔隙度、裂隙的纵横比、孔隙的纵横比、流体和岩石基质相同或相似的碳酸盐岩储层地震岩石物理模板；
目标储层弹性参数获取模块：用于获取目标储层的纵横波参数比数据和纵波阻抗数据；
预测模块：用于利用建立的碳酸盐岩储层地震岩石物理模板，基于目标储层的纵横波参数比数据和纵波阻抗数据将目标储层投点到建立的碳酸盐岩储层地震岩石物理模板上，从而对目标储层的孔隙填充物组成和/或裂隙密度进行预测。
14.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建装置，包括处理器及存储器；其中，存储器，用于存放计算机程序；
处理器，用于执行存储器上所存放的程序时，实现权利要求1-6任一项所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法的步骤。
15.一种裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测装置，包括处理器及存储器；其中，
存储器，用于存放计算机程序；
处理器，用于执行存储器上所存放的程序时，实现权利要求7所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法的步骤。
16.一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现权利要求1-6任一项所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法的步骤。
17.一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现权利要求7所述的裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层储层参数预测方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110109196.8A CN112946739B (zh) | 2021-01-27 | 2021-01-27 | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110109196.8A CN112946739B (zh) | 2021-01-27 | 2021-01-27 | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112946739A CN112946739A (zh) | 2021-06-11 |
| CN112946739B true CN112946739B (zh) | 2022-12-02 |
Family
ID=76237606
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110109196.8A Active CN112946739B (zh) | 2021-01-27 | 2021-01-27 | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112946739B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113868979B (zh) * | 2021-09-29 | 2023-03-24 | 中国石油化工股份有限公司 | 浅层气钻井风险的确定方法、装置和服务器 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130013209A1 (en) * | 2010-03-11 | 2013-01-10 | Yaping Zhu | Predicting anisotropic source rock properties from well data |
| EP2591381A4 (en) * | 2010-07-08 | 2015-12-02 | Geco Technology Bv | METHODS AND DEVICES FOR TRANSFORMING COLLECTED DATA TO IMPROVE VISUALIZATION CAPACITY |
| CN103760081B (zh) * | 2013-12-31 | 2016-01-06 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
| CN106054248B (zh) * | 2016-07-15 | 2017-07-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
| CN106597544B (zh) * | 2016-11-25 | 2019-01-18 | 中国石油天然气股份有限公司 | 致密油气藏脆性的预测方法及装置 |
| CN110968928B (zh) * | 2018-09-29 | 2024-03-22 | 北京阳光杰科科技股份有限公司 | 一种低渗致密砂岩等效三孔隙岩石物理建模方法及其应用 |
| CN112147689B (zh) * | 2020-08-27 | 2024-01-30 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层解释方法及解释量版建立方法 |
-
2021
- 2021-01-27 CN CN202110109196.8A patent/CN112946739B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112946739A (zh) | 2021-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Vik et al. | Evaluation of representative elementary volume for a vuggy carbonate rock—Part: Porosity, permeability, and dispersivity | |
| Tobin et al. | Reservoir quality modeling of tight-gas sands in Wamsutter field: Integration of diagenesis, petroleum systems, and production data | |
| Haldorsen | Simulator parameter assignment and the problem of scale in reservoir engineering | |
| US7676349B2 (en) | Integrated anisotropic rock physics model | |
| CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
| CN106154351A (zh) | 一种低孔渗储层渗透率的估算方法 | |
| Malki et al. | Underlying mechanisms and controlling factors of carbonate reservoir characterization from rock physics perspective: A comprehensive review | |
| RU2455483C2 (ru) | Способ оценки степени трещиноватости карбонатных пород через параметр диффузионно-адсорбционной активности | |
| CN110456412B (zh) | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 | |
| CN112946739B (zh) | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 | |
| Naseryan Moghadam et al. | Relative permeability and residual gaseous CO2 saturation in the Jurassic Brentskardhaugen Bed sandstones, Wilhelmøya Subgroup, western central Spitsbergen, Svalbard | |
| Chhatre et al. | A blind study of four digital rock physics vendor laboratories on porosity, absolute permeability, and primary drainage capillary pressure data on tight outcrops | |
| Malki et al. | The impact of thermal maturity on the organic-rich shales properties: A case study in Bakken | |
| US20220056792A1 (en) | Method and system of evaluating hydrocarbon in heterogeneous formation | |
| Matava et al. | Observing maturing source rocks on seismic reflection data | |
| Ringrose et al. | The property model | |
| CN110954942A (zh) | 粘土含量识别的图版制作方法及系统 | |
| RU2389875C1 (ru) | Способ определения геологических свойств терригенной породы | |
| Li | Optimization of multistage hydraulic fracturing treatment for maximization of the tight gas productivity | |
| CN113688515B (zh) | 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法 | |
| Zhang et al. | Missouri University of Science and Technology | |
| Lueck | Characterizing two carbonate formations for CO₂-EOR and carbon geosequestration: applicability of existing rock physics models and implications for feasibility of a time lapse monitoring program in the Wellington Oil Field, Sumner County, Kansas. | |
| Guo et al. | A graphical model for haloanhydrite components and P-wave velocity: a case study of haloanhydrites in Amu Darya Basin | |
| CN116413791A (zh) | 一种横波速度的估算方法和装置 | |
| CN114910499A (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 |