CN117316329B - 天然气水合物饱和度声电测井联合智能反演方法及系统 - Google Patents
天然气水合物饱和度声电测井联合智能反演方法及系统 Download PDFInfo
- Publication number
- CN117316329B CN117316329B CN202311570368.7A CN202311570368A CN117316329B CN 117316329 B CN117316329 B CN 117316329B CN 202311570368 A CN202311570368 A CN 202311570368A CN 117316329 B CN117316329 B CN 117316329B
- Authority
- CN
- China
- Prior art keywords
- hydrate
- resistivity
- pore
- model
- phase
- 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 114
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 title claims abstract description 106
- 239000011435 rock Substances 0.000 claims abstract description 57
- 238000012549 training Methods 0.000 claims abstract description 31
- 230000005012 migration Effects 0.000 claims abstract description 11
- 238000013508 migration Methods 0.000 claims abstract description 11
- 239000011148 porous material Substances 0.000 claims description 142
- 239000012071 phase Substances 0.000 claims description 125
- 239000002245 particle Substances 0.000 claims description 123
- 239000012530 fluid Substances 0.000 claims description 92
- 230000008569 process Effects 0.000 claims description 36
- 239000007790 solid phase Substances 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 32
- 239000007787 solid Substances 0.000 claims description 31
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 27
- 239000013049 sediment Substances 0.000 claims description 25
- 239000007789 gas Substances 0.000 claims description 24
- 239000011248 coating agent Substances 0.000 claims description 21
- 238000000576 coating method Methods 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 239000004568 cement Substances 0.000 claims description 16
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 15
- 239000011707 mineral Substances 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 14
- 239000010453 quartz Substances 0.000 claims description 13
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 13
- 150000004677 hydrates Chemical class 0.000 claims description 11
- 238000011835 investigation Methods 0.000 claims description 11
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000035699 permeability Effects 0.000 claims description 7
- 230000004913 activation Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- 239000002023 wood Substances 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims 1
- 238000012360 testing method Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 8
- 238000013528 artificial neural network Methods 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000013527 convolutional neural network Methods 0.000 description 5
- 238000011176 pooling Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000013526 transfer learning Methods 0.000 description 2
- 238000005056 compaction Methods 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0464—Convolutional networks [CNN, ConvNet]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/096—Transfer learning
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Computational Linguistics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Databases & Information Systems (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- Medical Informatics (AREA)
- Geophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地球物理测井技术领域，本发明涉及一种天然气水合物饱和度声电测井联合智能反演方法及系统，所述方法的步骤为：构建理论储层模型，确定储层物性参数；基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；构建深度残差网络模型，基于声学参数和电阻率建立的训练集训练深度残差网络模型；将训练后的深度残差网络模型对由目标工区储层的声学参数和电阻率构建的待测试集进行迁移学习，得到目标工区储层的天然气水合物饱和度。本发明能够快速、精确地预测天然气水合物饱和度，降低了解释结果的不确定性。
Description
技术领域
本发明属于地球物理测井技术领域，具体地说，涉及一种天然气水合物饱和度声电测井联合智能反演方法及系统。
背景技术
自然界中的天然气水合物（以下简称“水合物”）广泛存在于永久冻土带和陆缘外围的海底沉积物中，尤其后者具有分布广、资源量巨大、密度高等特点，自20世纪60年底至今一直被许多国家重点研究，被誉为清洁、高效、储量丰富的新型清洁能源。地球物理测井（以下简称为“测井”）是水合物识别的主要信息来源之一，可以揭示水合物及伴生游离气的原位性质和评价水合物储层参数，其中水合物饱和度的估算是其核心关键任务，对于水合物甜点区的优选和后续勘探有着至关重要的作用。
现有水合物饱和度评价方法普遍沿用针对压实成岩作用的油气储层模型，如阿尔奇公式（简称：Archie）或其衍生变形、威里平均时间公式（简称：Wyllie）、多相Wood方程等，对具有“弱胶结、细粒度、低孔径”未成岩特征的海底含水合物沉积物，上述饱和度计算模型估算结果往往与实际情况存在较大差异，难以满足水合物高效勘探开发需求。显然，上述模型不适用于水合物饱和度计算主要源于地质沉积特征的显著差异，给水合物饱和度评价和甜点预测带来了巨大挑战，因此，亟需发展新的、科学有效的水合物饱和度估算方法。
发明内容
本发明针对现有技术存在的等上述问题，提供一种天然气水合物饱和度声电测井联合智能反演方法及系统，能够快速、准确地预测天然气水合物饱和度，降低了解释结果的不确定性。
本发明第一方面，提供了一种天然气水合物饱和度声电测井联合智能反演方法，其具体步骤为：
S1、根据天然气水合物储层的文献调研结果，构建理论储层模型，确定储层物性参数；
S2、基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；
S3、基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；
S4、以声学参数和电阻率为输入、以天然气水合物饱和度为输出构建深度残差网络模型，根据步骤S2得到的声学参数和步骤S3得到的电阻率建立训练集，将训练集输入深度残差网络模型中进行训练得到训练后的深度残差网络模型；
S5、根据目标工区储层的声学参数和电阻率构建待测试集，将待测试集输入至训练后的深度残差网络模型中进行迁移学习，得到目标工区储层的天然气水合物饱和度。
在一些实施例中，所述步骤S3中，所述电阻率三相物理模型包括孔隙填充水合物的电阻率三相物理模型和胶结水合物的电阻率三相物理模型，其中：所述孔隙填充水合物的电阻率三相物理模型表示为：
其中：
式中，R T_p为孔隙填充水合物沉积物的电阻率，R f为孔隙流体的电阻率，ϕ为孔隙度，R m为将水合物颗粒与固体颗粒等效为一种颗粒的等效电阻率，R S为固体颗粒的电阻率，R H为天然气水合物颗粒的电阻率，F为电流偏离比例，G为平均几何因子，β S为固体颗粒的体积分数，β H为天然气水合物颗粒的体积分数，S H为天然气水合物饱和度；
所述胶结水合物的电阻率三相物理模型表示为：
式中，R T_c为胶结水合物沉积物的电阻率，R 2T_c为两相介质的等效电阻率，R f为孔隙流体的电阻率，β为孔隙水或流体相所占的体积分数，β=ϕ(1-S H)。
在一些实施例中，所述孔隙填充水合物的电阻率三相物理模型的构建方法为：
对于孔隙填充水合物，水合物颗粒相当于孔隙流体中的孤立颗粒，将水合物颗粒与固体颗粒等效为一种颗粒，这种颗粒的等效电阻率表示为：
则将三相介质转化为两相介质，得到孔隙填充水合物的电阻率上限R HS+_p和孔隙填充水合物的最大几何等效电阻率R geo_p表示为：
定义电流偏离比例为F，F=(1-ϕ)，根据电流偏离比例F、Hashin-Shtrikman电导率上边界和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
结合电导率与电阻率的反比关系，将孔隙填充水合物的电阻率上限和最大几何等效电阻率的表达式代入两相介质的真实电导率的表达式，得到孔隙填充水合物沉积物的真实电阻率表示为：
所述孔隙填充水合物沉积物的真实电阻率的表达式即为所述孔隙填充水合物的电阻率三相物理模型。
在一些实施例中，所述胶结水合物的电阻率三相物理模型的构建方法为：
对于胶结水合物，水合物颗粒堵塞相邻孔隙之间的连接，含水合物沉积物中只存在固体颗粒和水合物颗粒两相介质，即水合物颗粒代替了孔隙流体，则两相介质的电阻率上限R geo_c和电阻率下限R HS+_c分别表示为：
其中：
式中，R geo_c为两相介质的电阻率上限，R HS+_c为两相介质的电阻率下限，β H1为两相介质中水合物的体积分数；
定义电流偏离比例为F，F=(1-ϕ)，根据电流偏离比例F、Hashin-Shtrikman电导率上边界和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
结合电导率与电阻率的反比关系，将两相介质的电阻率上限和两相介质的电阻率下限的表达式代入两相介质的真实电导率的表达式，得到两相介质的等效电阻率表示为：
将孔隙水随机添加到水合物与固体颗粒组成的两相介质中，电流直接流经孔隙水，即流体相完全隔离，固体相全部连接，此时，Hashin-Shtrikman电导率下边界σ HS_可表示为：
式中，σ S和σ f分别表示固体相和流体相的电导率，β为孔隙水或流体相所占的体积分数；
结合电导率与电阻率的反比关系及胶结水合物固体相电导率σ S=1/R 2T_c，将两相介质的等效电阻率R 2T_c代入Hashin-Shtrikman电导率下边界σ HS_中，计算得到胶结型水合物沉积物的真实电阻率R T_c，所述胶结型水合物沉积物的真实电阻率表示为：
所述胶结型水合物沉积物的真实电阻率的表达式即为所述胶结型水合物的电阻率三相物理模型。
在一些实施例中，所述步骤S4中，所述深度残差网络模型通过短路的连接方式，直接把输入传到输出作为初始结果，输出结果表示为：
式中，为输出结果，/>为残差，X为输入。
在一些实施例中，所述步骤S4中，所述深度残差网络模型采用ResNet-18网络结构，在卷积层与激活层之间放置BN层，在卷积层的残差结构中采用ReLU函数作为激活函数。
在一些实施例中，所述步骤S2中，所述声学多尺度岩石物理模型的构建方法为：
S21、将喷射流模型加入到Kuster-Toksöz模型中，计算水合物颗粒的等效体积模量和水合物等效剪切模量；
S22、根据Wood公式计算含水合物沉积物的孔隙流体的等效体积模量和等效剪切模量分别表示为：
式中，K F为孔隙流体的等效体积模量，G F为孔隙流体的等效体积模量，S W为孔隙水饱和度，S G为游离气饱和度，S H为天然气水合物饱和度，K W为孔隙水的体积模量，K G为游离气的体积模量；
根据Hill平均方程计算固体相的等效体积模量和等效剪切模量表示为：
式中，K s为固体相的等效体积模量，G s为固体相的剪切模量；m为固体相中存在的矿物成分的种类数量，K i为固体相中第i种矿物成分的体积模量，G i为固体相中第i种矿物成分的剪切模量，f i为固体相中第i种矿物成分的体积分数；
S23、将步骤S21得到的水合物颗粒的等效体积模量和等效剪切模量代入胶结砂岩模型中，得到接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入胶结砂岩模型中，得到颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将步骤S22得到的固体相的体积模量和剪切模量代入Hashin-Shtrikman模型，得到颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入Hashin-Shtrikman模型中，得到孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量；
S24、联合步骤S22中得到的孔隙流体的体积模量和剪切模量以及步骤S23中得到的干岩石骨架的体积模量和剪切模量，根据Biot-Rayleigh双孔理论模型，得到接触胶结水合物、颗粒涂层水合物、颗粒支撑水合物、孔隙填充水合物四种不同赋存形态的含水合物沉积物的声波速度和声波衰减，即声学多尺度岩石物理模型。
在一些实施例中，在所述步骤S24中：
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到接触胶结水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒涂层水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒支撑水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到孔隙填充水合物沉积物的声波速度和声波衰减。
在一些实施例中，所述Biot-Rayleigh双孔理论模型的波动方程为：
式中，N、A、Q i、R i为刚性系数，i=1,2，u为固体平均粒子位移，、/>分别为u的一阶导数和二阶导数，/>、/>、/>为固体、孔1流体与孔2流体的位移场散度；/>,i=1,2,3，j=1,2,3是密度参数，/>是主相流体的密度，/>、/>分别是水饱和孔隙和水合物饱和孔隙的孔隙度，/>、/>是两个区域的局部孔隙度，/>,i=1,2是耗散参数，/>是局部流动过程引起的流体应变增量，/>、/>分别是/>的一阶导数和二阶导数，/>是主相流体的粘度，/>是主相流体的渗透率；若水合物作为孔隙流体，/>是孔隙流体为地层水时的位移，/>是孔隙流体为天然气水合物/游离气的位移；若水合物作为岩石骨架，/>是孔隙流体为在宿主骨架中的流体时的位移，/>是孔隙流体为在水合物包裹体中的流体时的位移；
根据平面波分析法求解上式，带入平面纵波的解析解，得到：
其中，表示接触胶结水合物或颗粒涂层水合物或颗粒支撑水合物或孔隙填充水合物的波数；/>、/>为方程的系数矩阵，i=1,2,3，j=1,2,3；
则声波速度和声波衰减的计算公式表示为：
其中，为含水合物沉积物的声波衰减；v为含水合物沉积物的声波速度；/>和/>分别取k的实部和虚部。
本发明第二方面，提供了一种天然气水合物饱和度声电测井联合智能反演系统，用于实现本发明第一方面所述天然气水合物饱和度声电测井联合智能反演方法，所述系统包括：
数据获取模块，用于获取天然气水合物储层的文献调研结果以及目标工区储层的声学参数和电阻率；
储层模型生成模块，用于根据天然气水合物储层的文献调研结果生成理论储层模型；
参数确定模块，用于确定储层物性参数；
声学参数计算模块，基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；
电阻率计算模块，基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；
数据集生成模块，根据声学参数计算模块得到的声学参数和电阻率计算模块得到的电阻率生成训练集，并将目标工区储层的声学参数和电阻率进行迁移学习得到待测试集；
模型生成及训练模块，用于以声学参数和电阻率为输入、以天然气水合物饱和度为输出生成深度残差网络模型，并根据训练集训练深度残差网络模型得到训练后深度残差网络模型；
饱和度计算模块，通过训练后深度残差网络模型对待测试集进行迁移学习，得到目标工区储层的天然气水合物饱和度。
与现有技术相比，本发明的优点和积极效果在于：
（1）本发明天然气水合物饱和度声电测井联合智能反演方法及系统，联合声学与电学两种属性的测井参数构建深度残差网络模型，通过深度残差网络模型反演计算天然气水合物饱和度，相较于现有基于固结成岩的传统油气储层饱和度模型，即可以利用声波和电阻率测井数据相互约束，还可以提高反演计算天然气水合物饱和度的精度及方法的鲁棒性，降低了解释结果的不确定性。
（2）本发明天然气水合物饱和度声电测井联合智能反演方法及系统，提供的声学多尺度岩石物理模型考虑了四种不同的水合物赋存形态，整合了喷射流模型和Biot-Rayleigh双孔理论模型，与实际含天然气水合物沉积物储层赋存形态一致，能够精确描述不同赋存形态的含天然气水合物沉积物的水合物饱和度的声学相应规律，为深度残差网络模型提供声学响应的训练集数据。
（3）本发明提供的天然气水合物饱和度声电测井联合智能反演方法及系统，避免了通过水合物原位岩石物理实验构建饱和度解释模型的困难，且随着水合物的勘探开发深入，测量数据的增加，本发明天然气水合物饱和度声电测井联合智能反演方法及系统构建的深度残差网络模型预测水合物饱和度更加准确。
附图说明
图1为本发明实施例所述天然气水合物饱和度声电测井联合智能反演方法的流程图；
图2为本发明实施例构建的10层层状天然气水合物理论储层模型示意图；
图3为本发明实施例所述天然气水合物饱和度声电测井联合智能反演系统的结构框图；
图4为本发明实施例计算得到的声波速度随天然气水合物饱和度的变化曲线图；
图5为本发明实施例计算得到的声波衰减随天然气水合物饱和度的变化曲线图；
图6为本发明实施例电阻率三相物理模型的理论模拟结果与实验结果对比图；
图7为本发明实施例本发明构建的深度残差网络模型结构示意图；
图8为本发明实施例ReLU函数及其导数图像示意图；
图9为本发明实施例采用本发明构建的深度残差网络模型对孔隙填充型的理论水合物储层的天然气水合物饱和度反演结果示意图；
图10为本发明实施例采用本发明构建的深度残差网络模型对颗粒支撑型的理论水合物储层的天然气水合物饱和度反演结果示意图；
图11为本发明实施例采用本发明构建的深度残差网络模型对接触胶结型的理论水合物储层的天然气水合物饱和度反演结果示意图；
图12为本发明实施例采用本发明构建的深度残差网络模型对颗粒涂层型的理论水合物储层的天然气水合物饱和度反演结果示意图；
图13为本发明实施例采用本发明构建的深度残差网络模型、单一声波和单一电阻率预测A井的天然水合物饱和度结果示意图；
图14为本发明实施例采用本发明构建的深度残差网络模型与现有人工神经网络（ANN）和卷积神经网络（CNN）预测A井的天然水合物饱和度结果示意图。
图中，1、数据获取模块，2、储层模型生成模块，3、参数确定模块，4、声学参数计算模块，5、电阻率计算模块，6、数据集生成模块，7、模型生成及训练模块，8、饱和度计算模块，9、最大池化层，10、全局平均池化层，11、全连接层。
具体实施方式
下面，结合附图通过示例性的实施方式对本发明进行具体描述。然而应当理解，在没有进一步叙述的情况下，一个实施方式中的元件、结构和特征也可以有益地结合到其他实施方式中。
现有天然气水合物饱和度评价解释技术存在两个难题：一、基于传统固结岩石的油气评价模型，与水合物实际地址特征不符，导致估算水合物饱和度与实际情况存在明显偏差；二、利用单一声波或单一电阻率测井数据，容易出现水合物饱和度解释的不确定性。为此，本发明提供了一种天然气水合物饱和度声电测井联合智能反演方法及系统，联合声学与电学两种属性的测井参数构建深度残差网络模型，通过深度残差网络模型反演计算天然气水合物饱和度，能够快速、准确地预测天然气水合物饱和度，降低了解释结果的不确定性。以下结合附图对本发明上述天然气水合物饱和度声电测井联合智能反演方法及系统进行详细说明。
参见图1，本发明第一方面实施例实施例提供了一种天然气水合物饱和度声电测井联合智能反演方法，其具体步骤为：
S1、根据天然气水合物储层的文献调研结果，构建理论储层模型，确定储层物性参数，所述储层物性参数包括孔隙度、渗透率、水合物饱和度。
具体地，通过对重点海域的调研，总结出了储层深度、储层厚度以及储层物性参数的主要变化范围。需要说明的是，调研的天然气水合物重点不限于上述三个海域，可根据实际需求进行选择。
根据文献调研结果，本发明一些实施例中，给定了10层的层状天然气水合物储层模型（参见图2），其中储层深度变化范围设置为200~300 mbsf（中文：海底以下深度，英文：meter below seafloor），厚度变化范围设置为5~15 m，储层中含水合物沉积物的矿物成分为石英，渗透率设置为100 mD，孔隙度变化范围设置为0.3~0.6，天然气水合物饱和度变化范围设置为30%~60%，具体各层的参数取值如表1所示。
表1
| 储层深度(mbsf) | 储层总厚度(m) | 孔隙度(%) | 饱和度(%) | 渗透率(mD) | |
| 01 | 200~210 | 10 | 30 | 40 | 100 |
| 02 | 210~215 | 5 | 40 | 30 | 100 |
| 03 | 215~225 | 10 | 50 | 60 | 100 |
| 04 | 225~240 | 15 | 60 | 50 | 100 |
| 05 | 240~245 | 5 | 40 | 30 | 100 |
| 06 | 245~255 | 10 | 30 | 60 | 100 |
| 07 | 255~270 | 15 | 50 | 40 | 100 |
| 08 | 270~280 | 10 | 40 | 30 | 100 |
| 09 | 280~285 | 5 | 30 | 30 | 100 |
| 10 | 285~300 | 15 | 60 | 50 | 100 |
S2、基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数，所述声学参数包括理论储层模型的声波速度和声波衰减。
具体地，所述声学多尺度岩石物理模型的构建方法为：
S21、将喷射流模型加入到Kuster-Toksöz模型中，计算水合物颗粒的等效体积模量和水合物等效剪切模量。
具体地，所述喷射流模型表示为：
式中，为包裹体的体积模量，/>为包裹体的剪切模量，K1i为水合物包裹体和孔隙中的流体等效体积模量，下标i=1,2,3,4,...,N，N为不同种类包裹体；K H为水合物颗粒的体积模量，τ i为弛豫时间，ω为角频率，μ i为水合物中包裹体的粘度。
将喷射流模型加入到Kuster-Toksöz模型中，计算得到水合物颗粒的等效体积模量和等效剪切模量表示为：
其中：
式中，为水合物颗粒的等效体积模量，/>为水合物颗粒的等效剪切模量，K H为水合物颗粒的体积模量，G H为水合物颗粒的剪切模量，下标H表示水合物颗粒，/>为第i中包裹体的体积分量，/>为包裹体的体积模量，/>为包裹体的剪切模量，ζ H为中间参数，P mi、Q mi为几何因子，α为水和甲烷气包裹体的半径，β m为水合物颗粒与固体颗粒等效为一种颗粒所占体积分数。
S22、根据wood公式计算含水合物沉积物的孔隙流体的等效体积模量和等效剪切模量分别表示为：
式中，K F为孔隙流体的等效体积模量，G F为孔隙流体的等效体积模量，S W为孔隙水饱和度，S G为游离气饱和度，S H为天然气水合物饱和度，K W为孔隙水的体积模量，K G为游离气的体积模量。
根据Hill平均方程计算固体相的等效体积模量和等效剪切模量表示为：
式中，K s为固体相的等效体积模量，G s为固体相的剪切模量；m为固体相中存在的矿物成分的种类数量，K i为固体相中第i种矿物成分的体积模量，G i为固体相中第i种矿物成分的剪切模量，f i为固体相中第i种矿物成分的体积分数。
S23、将步骤S21得到的水合物颗粒的等效体积模量和等效剪切模量代入胶结砂岩模型中，得到接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入胶结砂岩模型中，得到颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将步骤S22得到的固体相的体积模量和剪切模量代入Hashin-Shtrikman模型，得到颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入Hashin-Shtrikman模型中，得到孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量。
具体地，基于胶结砂岩模型得到干岩石骨架的体积模量和干岩石骨架的剪切模量表示为：
式中，K dry为干岩石骨架的体积模量，G dry为干岩石骨架的剪切模量，M c为胶结物或水合物颗粒的体积模量，G c为石英颗粒或水合物颗粒的剪切模量；参数为与胶结的两颗粒组合体的正向成比例，/>与胶结的两颗粒组合体的剪切刚度成比例，参数/>和/>取决于胶结物的含量及胶结物和骨架颗粒的特性；n为颗粒配位数，在含水合物沉积物中取n=8.5；ϕ 0为沉积物的初始孔隙度。
当孔隙度小于临界孔隙度时，所述Hashin-Shtrikman模型表示为：
当孔隙度大于等于临界孔隙度时，所述Hashin-Shtrikman模型表示为：
其中：
式中，ϕ为孔隙度；ϕ c为临界孔隙度，K HM为临界孔隙度下的体积模量，G HM为临界孔隙度下的剪切模量，λ为从K'和G'计算的矿物相的泊松比，K'为固体相或石英颗粒的体积模量，G'为固体相或石英颗粒的剪切模量，P为等效压力。
S24、联合步骤S22中得到的孔隙流体的体积模量和剪切模量以及步骤S23中得到的干岩石骨架的体积模量和剪切模量，根据Biot-Rayleigh双孔理论模型，得到接触胶结水合物、颗粒涂层水合物、颗粒支撑水合物、孔隙填充水合物四种不同赋存形态的含水合物沉积物的声波速度和声波衰减，即声学多尺度岩石物理模型。
具体地，将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到接触胶结水合物沉积物的声波速度和声波衰减。
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒涂层水合物沉积物的声波速度和声波衰减。
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒支撑水合物沉积物的声波速度和声波衰减。
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到孔隙填充水合物沉积物的声波速度和声波衰减。
具体地，所述Biot-Rayleigh双孔理论模型的波动方程为：
式中，N、A、Q i、R i为刚性系数，i=1,2，u为固体平均粒子位移，、/>分别为u的一阶导数和二阶导数，/>、/>、/>为固体、孔1流体与孔2流体的位移场散度；/>,i=1,2,3，j=1,2,3是密度参数，/>是主相流体的密度，/>、/>分别是水饱和孔隙和水合物饱和孔隙的孔隙度，/>、/>是两个区域的局部孔隙度，/>,i=1,2是耗散参数，/>是局部流动过程引起的流体应变增量，/>、/>分别是/>的一阶导数和二阶导数，/>是主相流体的粘度，/>是主相流体的渗透率；若水合物作为孔隙流体，/>是孔隙流体为地层水时的位移，/>是孔隙流体为天然气水合物/游离气的位移；若水合物作为岩石骨架，/>是孔隙流体为在宿主骨架中的流体时的位移，/>是孔隙流体为在水合物包裹体中的流体时的位移。
根据平面波分析法求解上式，带入平面纵波的解析解，得到：
其中，表示接触胶结水合物或颗粒涂层水合物或颗粒支撑水合物或孔隙填充水合物的波数；/>、/>为方程的系数矩阵，i=1,2,3，j=1,2,3；
则声波速度和声波衰减的计算公式表示为：
其中，为含水合物沉积物的声波衰减；v为含水合物沉积物的声波速度；/>和/>分别取k的实部和虚部。
需要说明的是，天然气水合物在沉积物中以不同的赋存形态存在。当水合物是孔隙流体的一部分时，它被添加到流体相中，形成孔隙填充水合物。当水合物是固体颗粒时，它被添加到固体基质中，形成颗粒支撑水合物。另外，水合物被当作胶结物，通过胶结砂岩理论中两种不同的胶结方式来计算干岩石骨架的弹性模量，形成接触胶结水合物或颗粒涂层水合物。
为了验证上述方法构建的声学多尺度岩石物理模型对于含天然气水合物沉积物的水合物不同赋存形态的预测能力，利用本发明实施例提供的上述声学多尺度岩石物理模型，选取了表2所示的声学多尺度岩石物理模型输入参数和表3所示的含水合物沉积物各组分弹性参数，计算得到声波速度随水合物饱和度的变化曲线（参见图4）和声波衰减随水合物饱和度的变化曲线（参见图5）。
表2
表3
| 体积模量/(Pa) | 剪切模量/(Pa) | 密度/(g/cm3) | |
| 石英 | 36.5 × 109 | 45.0 × 109 | 2650.0 |
| 水合物 | 7.9 × 109 | 3.3 × 109 | 910.0 |
| 孔隙水 | 2.17 × 109 | 0 | 1006.0 |
| 甲烷气 | 51.858 | 0 | 0.1334 |
图4、图5中的实验数据使用了“过量水”和“过量气”两种方法生成固体水合物，在砂岩岩样中测量了等效压力为500 KPa、频率为200 Hz下的声波速度和声波衰减。在“过量水”条件下，连续注入的水与一定量的甲烷气发生反应，生成少量不与颗粒接触的水合物，即孔隙填充水合物。在“过量气”条件下，注入大量的气与水反应，在颗粒接触处或颗粒表面形成气体水合物。通过对实验数据和模拟数据的比较，发现孔隙填充型的声波速度与实验数据的变化趋势吻合很好，准确率为97%，相对误差为2.3%，能够大致预测声波衰减峰值，如图4、图5所示。而颗粒涂层型的声波速度与实验数据吻合度很好，准确率为92%，相对误差为4.2%，但声波衰减峰值匹配较差。这是因为“过量气”方法生成的水合物是颗粒涂层、接触胶结和颗粒支撑的混合形态。此外，该声波多尺度岩石物理模型能捕捉测量到的声波速度和声波衰减随水合物饱和度的变化趋势，特别是在水合物饱和度为0%-40%的范围内。总之，基于本发明建立的声学多尺度岩石物理的数值模拟结果与实验数据的对比分析显示了两种结果良好的一致性，也说明本发明建立的声学多尺度岩石物理模型可以表征孔隙填充、颗粒支撑、接触胶结及颗粒涂层四种水合物赋存形态下的声学响应特征，为之后所需声学响应训练集提供了基础。
S3、基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率。
具体地，所述电阻率三相物理模型包括孔隙填充水合物的电阻率三相物理模型和胶结水合物的电阻率三相物理模型，其中：
所述孔隙填充水合物的电阻率三相物理模型表示为：
其中：
式中，R T_p为孔隙填充水合物沉积物的电阻率，R f为孔隙流体的电阻率，ϕ为孔隙度，R m为将水合物颗粒与固体颗粒等效为一种颗粒的等效电阻率，R S为固体颗粒的电阻率，R H为天然气水合物颗粒的电阻率，F为电流偏离比例，G为平均几何因子，β S为固体颗粒的体积分数，β H为天然气水合物颗粒的体积分数，S H为天然气水合物饱和度。
所述胶结水合物的电阻率三相物理模型表示为：
式中，R T_c为胶结水合物沉积物的电阻率，R 2T_c为两相介质的等效电阻率，R f为孔隙流体的电阻率，β为孔隙水或流体相所占的体积分数，β=ϕ(1-S H)。
在一些实施例中，所述孔隙填充水合物的电阻率三相物理模型的构建方法为：
对于孔隙填充水合物，水合物颗粒相当于孔隙流体中的孤立颗粒，将水合物颗粒与固体颗粒等效为一种颗粒，这种颗粒的等效电阻率表示为：
则将三相介质转化为两相介质，得到孔隙填充水合物的电阻率上限R HS+_p和孔隙填充水合物的最大几何等效电阻率R geo_p表示为：
定义电流偏离比例为F，F=(1-ϕ)，根据电流偏离比例F、Hashin-Shtrikman电导率上边界和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
需要说明的是，上述两相介质的真实电导率的表达式中，在流体相占比100%的情况下，真实电导率等于Hashin-Shtrikman电导率上边界/>，随着流体相占比的降低，真实电导率/>向最大几何等效电导率/>靠近。
结合电导率与电阻率的反比关系，将孔隙填充水合物的电阻率上限和最大几何等效电阻率的表达式代入两相介质的真实电导率的表达式，得到孔隙填充水合物沉积物的真实电阻率表示为：
所述孔隙填充水合物沉积物的真实电阻率的表达式即为所述孔隙填充水合物的电阻率三相物理模型。
在一些实施例中，所述胶结水合物的电阻率三相物理模型的构建方法为：
对于胶结水合物，水合物颗粒堵塞相邻孔隙之间的连接，含水合物沉积物中只存在固体颗粒和水合物颗粒两相介质，因水合物颗粒的电阻率远远小于固体颗粒的电阻率，电流会优先流经水合物，即水合物颗粒代替了孔隙流体，则两相介质的电阻率上限和电阻率下限分别表示为：
其中：
式中，R geo_c为两相介质的电阻率上限，R HS+_c为两相介质的电阻率下限，β H1为两相介质中水合物的体积分数；
定义电流偏离比例为F，F=(1-ϕ)，根据电流偏离比例F、Hashin-Shtrikman电导率上边界和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
结合电导率与电阻率的反比关系，将两相介质的电阻率上限和两相介质的电阻率下限的表达式代入两相介质的真实电导率的表达式，得到两相介质的等效电阻率表示为：
将孔隙水随机添加到水合物与固体颗粒组成的两相介质中，电流直接流经孔隙水，而不是绕过，即流体相完全隔离，固体相全部连接，此时，Hashin-Shtrikman电导率下边界σ HS_可表示为：
式中，σ S和σ f分别表示固体相和流体相的电导率，β为孔隙水或流体相所占的体积分数。
结合电导率与电阻率的反比关系及胶结水合物固体相电导率σ S=1/R 2T_c，将两相介质的等效电阻率R 2T_c代入Hashin-Shtrikman电导率下边界σ HS_中，计算得到胶结型水合物沉积物的真实电阻率R T_c，所述胶结型水合物沉积物的真实电阻率表示为：
所述胶结型水合物沉积物的真实电阻率的表达式即为所述胶结型水合物的电阻率三相物理模型。
为了验证本发明提供的电阻率三相岩石物理模型的表现和性能，对比了电阻率三相岩石物理模型的数值模拟结果与含水合物沉积物的电阻率实验数据以及常规的阿尔奇公式的效果（参见图6所示）。在计算中使用表4所示的模型输入参数，模型中孔隙水电阻率设置为0.1Ω·m，沉积物矿物组分为石英颗粒，电阻率为1×106Ω·m，水合物颗粒的电阻率为5000Ω·m，孔隙水的电阻率初始值设置为0.1Ω·m，其在0.02~20.1Ω·m的范围内变化。
表4
| 电阻率/(Ω·m) | |
| 石英固体颗粒 | 1 × 106 |
| 水合物颗粒 | 5000 |
| 孔隙水 | 0.02~2 |
从图6中可以看出，电阻率三相岩石物理模型在孔隙填充模型方面的表现相对较好，与孔隙填充水合物的数字岩心实验数据拟合度较高，准确率达到98%。这表明本发明建立的电阻率三相岩石物理模型能够准确地模拟孔隙填充水合物的电阻率响应规律。此外，胶结模型相对于颗粒包裹水合物的数字岩心实验数据拟合度较好，准确率达到96%，提供了一个更加准确地描述胶结水合物电学特性的模型。然而，阿尔奇公式的经验参数需要进行校正，尽管其在孔隙填充水合物方面的拟合效果较好，但其在胶结水合物方面的变化规律并不符合阿尔齐现象，这也进一步论证了适用于固结纯岩石的阿尔奇公式饱和度模型无法适用于“弱胶结、细粒度、低孔径”未成岩特征的海底含水合物沉积物的含水合物饱和度计算。
S4、构建深度残差网络模型，根据步骤S2得到的声学参数和步骤S3得到的电阻率建立训练集，将训练集输入深度残差网络模型中进行训练得到训练后深度残差网络模型。
在一些实施例中，所述深度残差网络模型采用ResNet-18网络结构（参见图7），所述ResNet-18网络结构包括三个主要部分：输入部分、输出部分和中间卷积部分。整个运算流程分为以下几个步骤：（1）数据进入深度残差网络后先进入输入部分，包括1个卷积层、1个BN层、ReLU激活函数以及最大池化层，实现数据的特征提取和降维；（2）随后进入中间的卷积部分，包括4个卷积层，每个卷积层由2个残差结构构成，通过多次特征提取和深层次学习，实现对数据更加复杂的特征提取和表示；（3）最后数据经过一个全局平均池化层和全连接层输出得到结果。其中，全局平均池化层能够有效地减少过拟合的风险，保证网络的稳定性和泛化性能。深度残差网络和其他卷积神经网络的差异主要在4个卷积层中的残差结构，这种结构通过短路的连接方式，直接把输入传到输出作为初始结果，输出结果为/>：
其中，表示残差。
ResNet-18网络结构将学习目标从一个完整的输出变成了残差，训练目标就是将残差结果逼近于0，使得随着网络加深，准确率不再下降。
在一些实施例中，在卷积层的残差结构中采用ReLU函数作为激活函数。ReLU函数是一种非线性激活函数，能够有效拟合复杂的函数关系，相较于Sigmoid函数，其计算量相对较小。如图8所示，ReLU函数在负半区的导数为0，使得神经网络具有良好的稀疏性，计算速度快，缓解了过拟合问题的发生。
在一些实施例中，在卷积层与激活层之间放置BN层，即“批归一化”，公式如下：
其中，表示输入当前批次的第/>个样本时该层第/>个输入节点的值，/>和/>为对应的均值和标准差，在正向传播时统计得到，/>表示防止除零的极小值，/>和/>表示待学习的参数，分别用于控制方差和均值，在反向传播过程中训练得到。BN层的加入减小了对参数初始化的依赖，更加方便调参；不需要dropout即可减小过拟合，提高了收敛效率。
S5、根据目标工区储层的声学参数和电阻率构建待测试集，将待测试集输入至训练后深度残差网络模型进行迁移学习，得到目标工区储层的天然气水合物饱和度。
需要说明的是，所述迁移学习解决由于训练集数据与目标储层数据（即待测试集数据）的特征空间不同而带来的问题。具体来说，利用迁移学习将训练集中的知识应用于目标储层数据，并尽可能减小两者之间的距离。在此过程中，将训练集称为源域，目标储层数据称为目标域。为了实现这一目标，采用迁移成分分析（简称：TCA）方法，TCA方法是一种最常见的迁移学习方法。通过使用矩阵变换，TCA方法能够在源域和目标域之间进行知识转移，并最小化它们之间的距离，以确保在目标域中获得最佳性能。
本发明上述天然气水合物饱和度声电测井联合智能反演方法，考虑了四种不同的赋存形态，联合了声-电两种储层测井参数构建深度残差网络模型，通过深度残差网络进行反演得到天然气水合物饱和度。相比现有反演方法，声-电联合反演结合了声学数据和电学两种物理性质的测井数据，能够同时反映岩性变化和流体状态，避免了单一属性数据反演计算多解性强的问题，提高了水合物饱和度的计算精度。
本发明第二方面实施例，提供了一种天然气水合物饱和度声电测井联合智能反演系统，用于实现本发明第一方面实施例所述天然气水合物饱和度声电测井联合智能反演方法，包括：
数据获取模块1，用于获取天然气水合物储层的文献调研结果以及目标工区储层的声学参数和电阻率；
储层模型生成模块2，用于根据天然气水合物储层的文献调研结果生成理论储层模型；
参数确定模块3，用于确定储层物性参数；所述储层物性参数包括孔隙度、水合物饱和度、渗透率；
声学参数计算模块4，基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；所述声学参数包括声波速度和声波衰减；
电阻率计算模块5，基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；
数据集生成模块6，根据声学参数计算模块得到的声学参数和电阻率计算模块得到的电阻率生成训练集，并将目标工区储层的声学参数和电阻率进行迁移学习得到待测试集；
模型生成及训练模块7，用于以声学参数和电阻率为输入、以天然气水合物饱和度为输出生成深度残差网络模型，并根据训练集训练深度残差网络模型得到训练后深度残差网络模型；
饱和度计算模块8，通过训练后深度残差网络模型预测待测试集得到目标工区储层的，得到目标工区储层的天然气水合物饱和度。
其中，声学多尺度岩石物理模型、电阻率三相物理模型的构建方法同本发明第一方面实施例所述天然气水合物饱和度声电测井联合智能反演方法，此处不再赘述。
本发明上述天然气水合物饱和度声电测井联合智能反演系统，考虑了四种不同的赋存形态，联合了声-电两种储层测井参数构建深度残差网络模型，通过深度残差网络进行反演得到天然气水合物饱和度。相比现有反演方法，声-电联合反演结合了声学数据和电学两种物理性质的测井数据，能够同时反映岩性变化和流体状态，避免了单一属性数据反演计算多解性强的问题，提高了水合物饱和度的计算精度。
以下结合具体实施例对本发明上述天然气水合物饱和度声电测井联合智能反演的有效性进行说明。
图9至图12显示了孔隙填充型、颗粒支撑型、接触胶结型、颗粒涂层型的理论水合物储层的天然气水合物饱和度反演结果。对比分析不同水合物赋存形态的反演结果，得出以下结论：（1）本发明上述方法和系统（以下简称：声电联合）反演天然气水合物饱和度结合了声学数据和电学数据的优势，能够同时反映岩性变化和流体状态，克服了单一属性数据反演的不足，提高了天然气水合物饱和度的反演精度。（2）不同水合物赋存形态的模拟结果显示声电联合反演天然气水合物饱和度的精确性也受到水合物赋存形态的影响。在孔隙填充型水合物赋存形态中，声电联合反演的天然气水合物饱和度的精度相对较高，能够准确地反演水合物饱和度。但在颗粒支撑型和接触胶结型水合物赋存形态中，声波在颗粒和基质之间发生反射和透射，容易出现波形变形和模拟不准确的情况，因此反演的天然气水合物饱和度精度相对较低。在颗粒涂层型水合物赋存形态中声电联合反演的天然气水合物饱和度的精度与孔隙填充型相当，但要比颗粒支撑型和接触胶结型略高一些。
表5
| 均方误差MSE | 平均绝对误差MAE | 决定系数R2 | |
| 孔隙填充水合物 | 0.000809 | 0.0206 | 0.929 |
| 颗粒支撑水合物 | 0.00107 | 0.0264 | 0.802 |
| 接触胶结水合物 | 0.00215 | 0.0405 | 0.813 |
| 颗粒涂层水合物 | 0.00238 | 0.0357 | 0.792 |
表6
| 均方误差MSE | 平均绝对误差MAE | 决定系数R2 | |
| 孔隙填充水合物 | 0.00193 | 0.0349 | 0.832 |
| 胶结水合物 | 0.00165 | 0.0338 | 0.856 |
表7
| 均方误差MSE | 平均绝对误差MAE | 决定系数R2 | |
| 孔隙填充水合物 | 0.000782 | 0.0235 | 0.932 |
| 颗粒支撑水合物 | 0.0228 | 0.0415 | 0.907 |
| 接触胶结水合物 | 0.00129 | 0.0301 | 0.887 |
| 颗粒涂层水合物 | 0.0357 | 0.0335 | 0.860 |
根据表5可以看出，孔隙填充水合物的单一声学数据反演精度可以达到92.9%，但其他三种赋存形态的反演精度不到81%。表6显示，孔隙填充水合物和胶结水合物的单一电阻率反演精度分别为83.2%和85.6%。在表7中，对于孔隙填充水合物，声-电联合的反演精度相较于单一声学数据反演提高了0.3%，相较于单一电阻率反演提高了10%。对于颗粒支撑水合物，声-电联合的反演精度相较于单一声学数据反演提高了10.5%；对于接触胶结水合物，声-电联合的反演精度相较于单一声学数据反演提高了7.4%；对于颗粒涂层水合物，声-电联合的反演精度相较于单一声学数据反演提高了6.8%。与单一电阻率反演结果相比，声-电联合反演精度也有所提升。
图8展示了A井的天然气水合物饱和度反演效果。可以看出，本发明上述天然气水合物饱和度声电测井联合智能反演方法及系统得到的天然气水合物饱和度与实际岩心测量结果数据相吻合，进一步验证了本发明提供的天然气水合物饱和度声电测井联合智能反演方法及系统可靠且适用于“弱胶结、细粒度、低孔径”的未成岩特征的海底含天然气水合物沉积物饱和度。
同时，为了证明本发明建立的深度残差网络模型（以下称为：ResNet-18）的优越性，分别选用人工神经网络（以下称为：ANN）和卷积神经网络（以下称为：CNN）分别对A井的天然气水合物饱和度进行计算，从图9中可以看出，本发明建立的深度残差网络模型在天然气水合物饱和度峰值的预测效果更好。表8展示了ANN、CNN以及ResNet-18对A井天然气水合物饱和度的反演精度和误差，由表8可以看出，ResNet-18的精度更高，误差更小。
表8
上述实施例用来解释本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。
Claims (8)
1.一种天然气水合物饱和度声电测井联合智能反演方法，其特征在于，其具体步骤为：
S1、根据天然气水合物储层的文献调研结果，构建理论储层模型，确定储层物性参数；
S2、基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；所述声学多尺度岩石物理模型的构建方法为：
S21、将喷射流模型加入到Kuster-Toksöz模型中，计算水合物颗粒的等效体积模量和水合物等效剪切模量；
S22、根据Wood公式计算含水合物沉积物的孔隙流体的等效体积模量和等效剪切模量分别表示为：
式中，K F为孔隙流体的等效体积模量，G F为孔隙流体的等效体积模量，S W为孔隙水饱和度，S G为游离气饱和度，S H为天然气水合物饱和度，K W为孔隙水的体积模量，K G为游离气的体积模量，K H为水合物颗粒的体积模量，G H为水合物颗粒的剪切模量；
根据Hill平均方程计算固体相的等效体积模量和等效剪切模量表示为：
式中，K s为固体相的等效体积模量，G s为固体相的剪切模量；m为固体相中存在的矿物成分的种类数量，K i为固体相中第i种矿物成分的体积模量，G i为固体相中第i种矿物成分的剪切模量，f i为固体相中第i种矿物成分的体积分数；
S23、将步骤S21得到的水合物颗粒的等效体积模量和等效剪切模量代入胶结砂岩模型中，得到接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入胶结砂岩模型中，得到颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将步骤S22得到的固体相的体积模量和剪切模量代入Hashin-Shtrikman模型，得到颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量；将石英颗粒的体积模量和剪切模量代入Hashin-Shtrikman模型中，得到孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量；
S24、联合步骤S22中得到的孔隙流体的体积模量和剪切模量以及步骤S23中得到的干岩石骨架的体积模量和剪切模量，根据Biot-Rayleigh双孔理论模型，得到接触胶结水合物、颗粒涂层水合物、颗粒支撑水合物、孔隙填充水合物四种不同赋存形态的含水合物沉积物的声波速度和声波衰减，即声学多尺度岩石物理模型；
S3、基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；所述电阻率三相物理模型包括孔隙填充水合物的电阻率三相物理模型和胶结水合物的电阻率三相物理模型，其中：
所述孔隙填充水合物的电阻率三相物理模型表示为：
其中：
式中，R T_p为孔隙填充水合物沉积物的电阻率，R f为孔隙流体的电阻率，为孔隙度，R m为将水合物颗粒与固体颗粒等效为一种颗粒的等效电阻率，R S为固体颗粒的电阻率，R H为天然气水合物颗粒的电阻率，F为电流偏离比例，G为平均几何因子，β S为固体颗粒的体积分数，β H为天然气水合物颗粒的体积分数，S H为天然气水合物饱和度；
所述胶结水合物的电阻率三相物理模型表示为：
式中，R T_c为胶结水合物沉积物的电阻率，R 2T_c为两相介质的等效电阻率，R f为孔隙流体的电阻率，β为孔隙水或流体相所占的体积分数，；
S4、以声学参数和电阻率为输入、以天然气水合物饱和度为输出构建深度残差网络模型，根据步骤S2得到的声学参数和步骤S3得到的电阻率建立训练集，将训练集输入深度残差网络模型中进行训练得到训练后的深度残差网络模型；
S5、根据目标工区储层的声学参数和电阻率构建待测试集，将待测试集输入至训练后的深度残差网络模型中进行迁移学习，得到目标工区储层的天然气水合物饱和度。
2.如权利要求1所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，所述孔隙填充水合物的电阻率三相物理模型的构建方法为：
对于孔隙填充水合物，水合物颗粒相当于孔隙流体中的孤立颗粒，将水合物颗粒与固体颗粒等效为一种颗粒，这种颗粒的等效电阻率表示为：
则将三相介质转化为两相介质，得到孔隙填充水合物的电阻率上限R HS+_p和孔隙填充水合物的最大几何等效电阻率R geo_p表示为：
定义电流偏离比例为F，，根据电流偏离比例F、Hashin-Shtrikman电导率上边界/>和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
结合电导率与电阻率的反比关系，将孔隙填充水合物的电阻率上限和最大几何等效电阻率的表达式代入两相介质的真实电导率的表达式，得到孔隙填充水合物沉积物的真实电阻率表示为：
所述孔隙填充水合物沉积物的真实电阻率的表达式即为所述孔隙填充水合物的电阻率三相物理模型。
3.如权利要求1所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，所述胶结水合物的电阻率三相物理模型的构建方法为：
对于胶结水合物，水合物颗粒堵塞相邻孔隙之间的连接，含水合物沉积物中只存在固体颗粒和水合物颗粒两相介质，即水合物颗粒代替了孔隙流体，则两相介质的电阻率上限和电阻率下限分别表示为：
其中：
式中，R geo_c为两相介质的电阻率上限，R HS+_c为两相介质的电阻率下限，β H1为两相介质中水合物的体积分数；
定义电流偏离比例为F，，根据电流偏离比例F、Hashin-Shtrikman电导率上边界/>和最大几何等效电导率/>计算两相介质的真实电导率/>，所述两相介质的真实电导率/>表示为：
结合电导率与电阻率的反比关系，将两相介质的电阻率上限和两相介质的电阻率下限的表达式代入两相介质的真实电导率的表达式，得到两相介质的等效电阻率表示为：
将孔隙水随机添加到水合物与固体颗粒组成的两相介质中，电流直接流经孔隙水，即流体相完全隔离，固体相全部连接，此时，Hashin-Shtrikman电导率下边界σ HS_可表示为：
式中，σ S和σ f分别表示固体相和流体相的电导率，β为孔隙水或流体相所占的体积分数；
结合电导率与电阻率的反比关系及胶结水合物固体相电导率σ S=1/R 2T_c，将两相介质的等效电阻率R 2T_c代入Hashin-Shtrikman电导率下边界σ HS_中，计算得到胶结型水合物沉积物的真实电阻率R T_c，所述胶结型水合物沉积物的真实电阻率表示为：
所述胶结型水合物沉积物的真实电阻率的表达式即为所述胶结型水合物的电阻率三相物理模型。
4.如权利要求1所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，所述步骤S4中，所述深度残差网络模型通过短路的连接方式，直接把输入传到输出作为初始结果，输出结果表示为：
式中，为输出结果，/>为残差，X为输入。
5.如权利要求1所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，所述步骤S4中，所述深度残差网络模型采用ResNet-18网络结构，在卷积层与激活层之间放置BN层，在卷积层的残差结构中采用ReLU函数作为激活函数。
6.如权利要求1所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，在所述步骤S24中：
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的接触胶结水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到接触胶结水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒涂层水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒涂层水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的颗粒支撑水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到颗粒支撑水合物沉积物的声波速度和声波衰减；
将步骤S22得到的孔隙流体的体积模量和剪切模量以及步骤S23得到的孔隙填充水合物的干岩石骨架体积模量和干岩石骨架剪切模量代入Biot-Rayleigh双孔理论模型，得到孔隙填充水合物沉积物的声波速度和声波衰减。
7.如权利要求6所述的天然气水合物饱和度声电测井联合智能反演方法，其特征在于，所述Biot-Rayleigh双孔理论模型的波动方程为：
式中，N、A、Q i、R i为刚性系数，i=1,2，u为固体平均粒子位移，、/>分别为u的一阶导数和二阶导数，/>、/>、/>为固体、孔1流体与孔2流体的位移场散度；/>,i=1,2,3，j=1,2,3是密度参数，/>是主相流体的密度，/>、/>分别是水饱和孔隙和水合物饱和孔隙的孔隙度，/>、/>是两个区域的局部孔隙度，/>,i=1,2是耗散参数，/>是局部流动过程引起的流体应变增量，/>、/>分别是/>的一阶导数和二阶导数，/>是主相流体的粘度，/>是主相流体的渗透率；若水合物作为孔隙流体，/>是孔隙流体为地层水时的位移，/>是孔隙流体为天然气水合物/游离气的位移；若水合物作为岩石骨架，/>是孔隙流体为在宿主骨架中的流体时的位移，/>是孔隙流体为在水合物包裹体中的流体时的位移；
根据平面波分析法求解上式，带入平面纵波的解析解，得到：
其中，表示接触胶结水合物或颗粒涂层水合物或颗粒支撑水合物或孔隙填充水合物的波数；/>、/>为方程的系数矩阵，i=1,2,3，j=1,2,3；
则声波速度和声波衰减的计算公式表示为：
其中，为含水合物沉积物的声波衰减；v为含水合物沉积物的声波速度；/>和分别取k的实部和虚部。
8.一种天然气水合物饱和度声电测井联合智能反演系统，用于实现如权利要求1-7任意一项所述天然气水合物饱和度声电测井联合智能反演方法，其特征在于，包括：
数据获取模块，用于获取天然气水合物储层的文献调研结果以及目标工区储层的声学参数和电阻率；
储层模型生成模块，用于根据天然气水合物储层的文献调研结果生成理论储层模型；
参数确定模块，用于确定储层物性参数；
声学参数计算模块，基于声学多尺度岩石物理模型对理论储层模型进行正演模拟，得到理论储层模型的声学参数；
电阻率计算模块，基于电阻率三相物理模型对理论储层模型进行正演模拟，得到理论储层模型的电阻率；
数据集生成模块，根据声学参数计算模块得到的声学参数和电阻率计算模块得到的电阻率生成训练集，并将目标工区储层的声学参数和电阻率进行迁移学习得到待测试集；
模型生成及训练模块，用于以声学参数和电阻率为输入、以天然气水合物饱和度为输出生成深度残差网络模型，并根据训练集训练深度残差网络模型得到训练后深度残差网络模型；
饱和度计算模块，通过训练后深度残差网络模型预对测试集进行迁移学习，得到目标工区储层的天然气水合物饱和度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311570368.7A CN117316329B (zh) | 2023-11-23 | 2023-11-23 | 天然气水合物饱和度声电测井联合智能反演方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311570368.7A CN117316329B (zh) | 2023-11-23 | 2023-11-23 | 天然气水合物饱和度声电测井联合智能反演方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117316329A CN117316329A (zh) | 2023-12-29 |
| CN117316329B true CN117316329B (zh) | 2024-03-29 |
Family
ID=89281326
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311570368.7A Active CN117316329B (zh) | 2023-11-23 | 2023-11-23 | 天然气水合物饱和度声电测井联合智能反演方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117316329B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111781637A (zh) * | 2019-04-04 | 2020-10-16 | 中国石油天然气集团有限公司 | 近海底水合物储层建模方法及装置 |
| CN112133377A (zh) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种海域天然气水合物赋存状态的判别方法及系统 |
| CN113189674A (zh) * | 2021-05-07 | 2021-07-30 | 中国石油大学(北京) | 一种天然气水合物的饱和度估算方法及系统 |
| CN115308793A (zh) * | 2022-08-05 | 2022-11-08 | 青岛海洋地质研究所 | 一种浅表层天然气水合物富集区智能预测方法 |
| CN115993665A (zh) * | 2023-03-06 | 2023-04-21 | 中国石油大学(华东) | 基于多尺度岩石物理模型的井中水合物饱和度计算方法 |
| CN117031545A (zh) * | 2023-07-28 | 2023-11-10 | 中国海洋大学 | 一种识别砂质储层水合物与游离气共存层的敏感参数方法 |
| CN117055114A (zh) * | 2023-10-09 | 2023-11-14 | 中国石油大学(华东) | 储层沉积物游离气饱和度定量分析方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114109359B (zh) * | 2021-11-16 | 2022-06-17 | 广州海洋地质调查局 | 海底水合物储层垂向含量分布精准评估装置使用方法 |
-
2023
- 2023-11-23 CN CN202311570368.7A patent/CN117316329B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111781637A (zh) * | 2019-04-04 | 2020-10-16 | 中国石油天然气集团有限公司 | 近海底水合物储层建模方法及装置 |
| CN112133377A (zh) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种海域天然气水合物赋存状态的判别方法及系统 |
| CN113189674A (zh) * | 2021-05-07 | 2021-07-30 | 中国石油大学(北京) | 一种天然气水合物的饱和度估算方法及系统 |
| CN115308793A (zh) * | 2022-08-05 | 2022-11-08 | 青岛海洋地质研究所 | 一种浅表层天然气水合物富集区智能预测方法 |
| CN115993665A (zh) * | 2023-03-06 | 2023-04-21 | 中国石油大学(华东) | 基于多尺度岩石物理模型的井中水合物饱和度计算方法 |
| CN117031545A (zh) * | 2023-07-28 | 2023-11-10 | 中国海洋大学 | 一种识别砂质储层水合物与游离气共存层的敏感参数方法 |
| CN117055114A (zh) * | 2023-10-09 | 2023-11-14 | 中国石油大学(华东) | 储层沉积物游离气饱和度定量分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 天然气水合物储层测井评价及其影响因素;宁伏龙;刘力;李实;张可;蒋国盛;吴能友;孙长宇;陈光进;;石油学报(第03期);第591-606页 * |
| 天然气水合物地层岩石物理模型构建;刘欣欣;印兴耀;栾锡武;;中国科学:地球科学(第09期);第1248-1266页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117316329A (zh) | 2023-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Carcione et al. | Computational poroelasticity—A review | |
| Chen et al. | Data assimilation for nonlinear problems by ensemble Kalman filter with reparameterization | |
| WO2010118624A1 (zh) | 地震波速度场构建中井约束横向可变h-v曲线构建方法 | |
| WO2005119276A2 (en) | Method for predicting lithology and porosity from seismic reflection data | |
| CN115220101B (zh) | 一种深层碎屑岩砂体沉积结构的建模方法 | |
| CN105804737A (zh) | 一种基于迭代算法求解地层孔隙度的方法 | |
| Li et al. | Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance | |
| CN115993665A (zh) | 基于多尺度岩石物理模型的井中水合物饱和度计算方法 | |
| Evazi et al. | Dual-porosity coarse-scale modeling and simulation of highly heterogeneous geomodels | |
| Luo et al. | Probabilistic seismic petrophysical inversion with statistical double-porosity Biot-Rayleigh model | |
| Guo et al. | Reservoir modeling, history matching, and characterization with a reservoir graph network model | |
| Pan et al. | A unified effective medium modeling framework for quantitative characterization of hydrate reservoirs | |
| CN108665545B (zh) | 录井参数三维地质模型建立方法 | |
| CN117316329B (zh) | 天然气水合物饱和度声电测井联合智能反演方法及系统 | |
| Yanran et al. | Mathematical modeling for total organic carbon content prediction with logging parameters by neural networks: A case study of shale gas well in South China | |
| Pan et al. | Joint elastic-electrical effective medium modeling for improved hydrate quantification | |
| Yanshu et al. | Muddy interlayer forecasting and an equivalent upscaling method based on tortuous paths: A case study of Mackay River oil sand reservoirs in Canada | |
| Fursov | Quantitative application of 4D seismic data for updating thin-reservoir models | |
| CN114325832B (zh) | 一种裂缝参数和弹性参数同步反演方法及系统 | |
| CN116256797A (zh) | 一种水平分层双重孔隙介质中点源激发地震波场正演方法 | |
| CN113534263B (zh) | 一种不依赖测井资料的含油气饱和度预测方法 | |
| CN113281825B (zh) | 岩石物理模型构建方法及装置 | |
| Walker et al. | Stochastic inversion for facies: A case study on the Schiehallion field | |
| CN104181090A (zh) | 一种基于岩电结构系数的含油气饱和度评价方法 | |
| Pan et al. | Joint Interpretation of elastic and electrical data for petrophysical properties of gas-hydrate-bearing sediments using inverse rock physics modeling method |
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 |