CN112558177A - 一种碎屑岩等效孔隙结构反演方法及系统 - Google Patents
一种碎屑岩等效孔隙结构反演方法及系统 Download PDFInfo
- Publication number
- CN112558177A CN112558177A CN201910908097.9A CN201910908097A CN112558177A CN 112558177 A CN112558177 A CN 112558177A CN 201910908097 A CN201910908097 A CN 201910908097A CN 112558177 A CN112558177 A CN 112558177A
- Authority
- CN
- China
- Prior art keywords
- pore structure
- equivalent pore
- unit
- wave velocity
- parameters
- 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.)
- Pending
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 168
- 239000011435 rock Substances 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 65
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000009826 distribution Methods 0.000 claims description 54
- 238000004088 simulation Methods 0.000 claims description 51
- 238000012360 testing method Methods 0.000 claims description 8
- 230000001172 regenerating effect Effects 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 238000011156 evaluation Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 238000009533 lab test Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种碎屑岩等效孔隙结构反演方法及系统，属于岩石物理学及地球物理测井领域。该碎屑岩等效孔隙结构反演方法利用目的层段的多个岩芯样品的等效孔隙结构参数和测井数据中的纵波速度、横波速度、孔隙度、泥质含量、含水饱和度反演出该目的层段的等效孔隙结构参数。本发明所需要的参数少，只需要纵、横波速度、孔隙度、泥质含量和含水饱和度，就可以预测出岩石等效孔隙结构，避免了传统的经验公式中因统计关系而存在的不确定性，应用范围广，具有计算速度快、稳定性好的优点，可以为工区内测井横波速度预测提供重要参数。
Description
技术领域
本发明属于岩石物理学及地球物理测井领域，具体涉及一种碎屑岩等效孔隙结构反演方法及系统。
背景技术
流体是直接赋存在岩石孔隙中的，所以基于地震资料的孔隙结构研究在油气/流体识别方面有时比弹性三参数(纵、横波速度和岩石密度)反演更直接、有效。但目前对于从测井资料或者从地震资料预测岩石孔隙结构的研究较少。
叠前AVO分析技术和叠前反演技术在储层预测中发挥了重要作用，这些技术方法均需要准确的纵波及横波速度测井曲线，然而在实际测井资料中往往缺乏横波速度信息，这给储层预测工作带来了一定的困难。为此，许多地球物理工作者对测井横波速度预测进行了大量的研究工作，发展了经验公式法和理论模型法两大类测井横波速度预测方法。其中经验公式法基于统计规律，其应用效果依赖于统计样品，受地域性影响，不具普遍性。理论模型法基于岩石物理方程，由于具有较高的预测精度，取得了广泛的应用。其中Xu-White速度模型法(或称为计算流程)预测横波速度精度较高，但是由于该模型需要输入的参数较多，有些参数对预测精度影响大却不能从实验室或者测井数据直接获得，从而限制了该方法的广泛应用。该计算流程综合考虑了泥质砂岩中基质性质、泥质含量、孔隙结构以及孔隙饱含流体性质对岩石速度的影响，这些参数中，泥质砂岩中基质性质可由岩石物理实验室实验观测及等效介质理论相关公式得到、泥质含量可由测井数据得到，而孔隙结构无法直接从测井数据得到。数值试验表明，在Xu-White速度模型法中，孔隙结构对预测的横波速度有重要影响，如果在较大深度范围内采用固定的孔隙结构参数，那么会对预测的横波引入较大的误差。对于岩石孔隙结构的研究，主要是通过地球物理实验室实验方法和测井资料现场评价法来实现，集中在岩石所具有的孔隙和喉道的几何形状、大小、分布及其相互连通关系等方面。其中室内实验法是目前最主要和最广泛的描述和评价岩石孔隙结构的方法，包括毛管压力曲线法、铸体薄片法、扫描电镜法及CT扫描法。在测井资料现场评价法中，有电阻率法、常规孔渗法、核磁共振法。上述两种孔隙结构特征描述方法均有各自的局限性，实验室实验法研究岩石孔隙结构要求必须有取芯岩样，容易受到岩石样品尺寸大小的限制，在岩芯破碎的层段是就无法描述孔隙结构，不具有纵向上的连续性，不能很好地反映一定地区的储层孔喉结构特征，另外实验室实验法研究岩石孔隙结构的费用昂贵。测井资料现场评价法能够沿着井壁连续的对孔隙结构进行评价，但是这些方法不能定量的给定岩石孔隙结构参数。孔隙结构特征描述方法中的实验法和测井资料现场评价法由于各自的局限性，所描述的孔隙结构参数都无法应用到横波速度预测中来提高横波速度预测精度，使得孔隙结构参数的应用范围受到限制。
综上所述，常规岩石孔隙结构特征描述方法存在受限于钻井岩芯质量、不具有纵向上的连续性，不适于横波速度预测方法，应用受到限制等缺点。
发明内容
本发明的目的在于解决上述现有技术中存在的难题，提供一种碎屑岩等效孔隙结构反演方法及系统，利用随机模拟和岩石物理模型结合，综合利用常规测井资料反演等效孔隙结构，为工区内测井横波速度预测提供重要参数。
本发明是通过以下技术方案实现的：
一种碎屑岩等效孔隙结构反演方法，所述方法利用目的层段的多个岩芯样品的等效孔隙结构参数和测井数据中的纵波速度、横波速度、孔隙度、泥质含量、含水饱和度反演出该目的层段的等效孔隙结构参数。
所述方法包括：
(1)根据岩芯样品测试数据统计等效孔隙结构参数的数值分布范围及分布规律；
(2)利用所述数值分布范围及分布规律获得随机模拟等效孔隙结构参数，并计算得到纵波速度和横波速度；
(3)利用计算得到的纵波速度和横波速度以及测井数据中的纵波速度和横波速度建立误差函数f(αs0)；
(4)判断f(αs0)≤ε是否成立，如果是，则转入步骤(5)，如果否，则返回步骤(2)；
(5)将所述随机模拟等效孔隙结构参数作为等效孔隙结构参数的反演结果进行输出。
所述步骤(1)的操作包括：
输入测井数据，所述测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；
统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[Pl，Pu]，其中，Pl为岩芯样品的等效孔隙结构参数中的最小值，Pu为岩芯样品的等效孔隙结构参数中的最大值。
所述步骤(2)的操作包括：
按照所述分布规律在数值分布范围[Pl，Pu]内随机生成一个随机模拟等效孔隙结构参数αs0；
利用所述随机模拟等效孔隙结构参数αs0和孔隙度、泥质含量、含水饱和度，由岩石物理模型计算得到纵波速度Vp(αs0)、横波速度Vs(αs0)；
所述岩石物理模型采用的是适用于碎屑岩的岩石物理模型。
所述步骤(3)的操作包括：
利用公式(1)建立误差函数f(αs0)：
f(αs0)＝||Vp(αs0)-Vpm||+||Vs(αs0)-Vsm||+|Vp(αs0)/Vs(αs0)-Vpm/Vsm| (1)
公式(1)中，Vpm为测井数据中的纵波速度，Vsm为测井数据中的横波速度；
本发明还提供一种碎屑岩等效孔隙结构反演系统，包括：
测井数据输入单元：用于输入测井数据，测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；
统计单元：用于统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[Pl,Pu]，其中，Pl为岩芯样品的等效孔隙结构参数中的最小值，Pu为岩芯样品的等效孔隙结构参数中的最大值；将数值分布范围及分布规律、开始生成参数的信息发送给随机模拟等效孔隙结构参数生成单元；
随机模拟等效孔隙结构参数生成单元：分别与所述统计单元、判断单元连接，用于在接收到统计单元发送来的开始生成参数的信息后，或者在接收到判断单元发送来的重新生成随机模拟等效孔隙结构参数的消息后，按照统计单元发送来的分布规律在数值分布范围[Pl,Pu]内随机生成一个随机模拟等效孔隙结构参数αs0；
速度计算单元：分别与测井数据输入单元、随机模拟等效孔隙结构参数生成单元连接，用于利用随机模拟等效孔隙结构参数生成单元发送来的随机模拟等效孔隙结构参数αs0和测井数据输入单元发送来的测井数据，由岩石物理模型计算得到纵波速度Vp(αs0)、横波速度Vs(αs0)；
误差函数建立单元：分别与测井数据输入单元、速度计算单元连接，用于建立误差函数f(αs0)；
判断单元：分别与随机模拟等效孔隙结构参数生成单元、误差函数建立单元、反演结果输出单元连接，用于判断f(αs0)≤ε是否成立，如果是，则发送输出结果的消息给反演结果输出单元，如果否，则发送重新生成随机模拟等效孔隙结构参数的消息给随机模拟等效孔隙结构参数生成单元；
反演结果输出单元：分别与随机模拟等效孔隙结构参数生成单元、判断单元连接，用于在接收到判断单元发送来的输出结果的消息后，将随机模拟等效孔隙结构参数生成单元中最新的随机模拟等效孔隙结构参数作为反演结果进行输出。
所述速度计算单元采用的岩石物理模型是适用于碎屑岩的岩石物理模型。
所述误差函数建立单元利用公式(1)建立误差函数f(αs0)：
f(αs0)＝||Vp(αs0)-Vpm||+||Vs(αs0)-Vsm||+|Vp(αs0)/Vs(αs0)-Vpm/Vsm| (1)
公式(1)中，Vpm为测井数据中的纵波速度，Vsm为测井数据中的横波速度；
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行的至少一个程序，所述至少一个程序被所述计算机执行时使所述计算机执行本发明的碎屑岩等效孔隙结构反演方法中的步骤。
与现有技术相比，本发明的有益效果是：本发明是一种碎屑岩等效孔隙结构反演方法，利用随机模拟和岩石物理模型结合从常规测井资料中提取岩石等效孔隙结构参数。本发明方法所需要的参数少，只需要纵、横波速度、孔隙度、泥质含量和含水饱和度，就可以预测出岩石等效孔隙结构，避免了传统的经验公式中因统计关系而存在的不确定性，应用范围广。具有计算速度快、稳定性好的优点，可以为工区内测井横波速度预测提供重要参数。
附图说明
图1是实施例1中的输入测井数据；
图2是实施例1中反演的孔隙结构参数；
图3是实施例1中孔隙结构参数与自然伽马交会图；
图4本发明方法的步骤框图；
图5本发明系统的组成结构图。
具体实施方式
下面结合附图对本发明作进一步详细描述：
本发明能够利用常规测井数据反演岩石等效孔隙结构参数，在等效孔隙结构参数可行解范围内随机模拟生成等效孔隙结构参数，然后利用岩石物理模型计算出当前孔隙结构参数下岩石的纵波速度和横波速度，然后将该计算的纵波速度和横波速度与测井的纵波速度和横波速度进行比较，优选出满足精度要求的孔隙结构参数。具有计算速度快、稳定性好的优点，可直接用于测井数据解释工作。
具体的，本发明方法如图4所示，包括以下步骤：
(1)根据岩芯样品测试数据统计等效孔隙结构分布范围及分布规律，具体如下：输入测井数据，包括纵波速度、横波速度、泥质含量、含水饱和度、孔隙度等测井曲线；统计岩芯样品的等效孔隙结构参数(采用现有的测试化验方法对岩芯样品进行测试化验获得等效孔隙结构参数，在此不再赘述。)的数值分布范围及分布规律，该范围表示为[Pl,Pu]，Pl为统计样品的等效孔隙结构参数的最小值，Pu为统计样品的等效孔隙结构参数的最大值。用于统计等效孔隙结构参数的数值分布范围及分布规律的岩芯样品选自研究的目的层段。
(2)随机模拟等效孔隙结构参数，具体如下：按照统计岩芯样品等效孔隙结构参数数值的分布规律在参数范围[Pl,Pu]内随机生成一个等效孔隙结构参数αs0(采用现有的多种用于生成服从某种分布规律的随机数的函数在[Pl,Pu]内生成一个随机数即可)，利用岩石物理模型计算等效孔隙结构参数为αs0时的纵波速度和横波速度。
所述的岩石物理模型是指适用于碎屑岩的岩石物理模型(利用现有技术中的岩石物理模型，例如采用KT模型，在此不再赘述)；所述的利用岩石物理模型计算纵波速度和横波速度，计算时需要等效孔隙结构参数αs0、孔隙度、泥质含量、含水饱和度；所述的孔隙度、泥质含量、含水饱和度为常规测井参数，由测井解释得到，这里作为已知参数给定；
(3)基于岩石物理方程正演计算纵、横波速度，并建立误差函数，具体如下：利用上述步骤2随机模拟等效孔隙结构参数，由岩石物理模型计算得到纵波速度、横波速度，将计算得到的纵波速度、横波速度与测井实测的纵波速度和横波速度进行比较，利用公式(1)得到误差函数f(αs0)：
f(αs0)＝||Vp(αs0)-Vpm||+||Vs(αs0)-Vsm||+|Vp(αs0)/Vs(αs0)-Vpm/Vsm| (1)
所述公式(1)中，Vpm为实测纵波速度，由测井数据得到。Vsm为实测横波速度，由测井数据得到。Vp(αs0)是由岩石物理模型计算的纵波速度，Vs(αs0)是由岩石物理模型计算的横波速度。
所述公式(1)中，||g||表示二范数，|g|表示一范数，此处用g代替了公式(1)中的符号内的表达式；
f(αs0)为误差函数，表示随机模拟等效孔隙结构参数为αs0时纵横波速度与测井实测纵横波速度的比较，该值越小表示理论计算值与实测值越接近。
(4)判断f(αs0)≤ε是否成立(ε是设定的阈值，根据速度数量级大小确定的，如速度分布为3000-6000m/s，那么该阀值可以设置为速度的1％)，如果是，则找到了使误差函数最小的随机模拟等效孔隙结构参数，转入步骤(5)，如果否，则返回步骤(2)；
(5)将该随机模拟等效孔隙结构参数作为等效孔隙结构参数的反演结果进行输出。
如图5所示，本发明还提供一种碎屑岩等效孔隙结构反演系统，包括：测井数据输入单元10、统计单元20、随机模拟等效孔隙结构参数生成单元30、速度计算单元40、误差函数建立单元50、判断单元60、反演结果输出单元70；
测井数据输入单元10：用于输入测井数据；
统计单元20：用于统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[Pl,Pu]，其中，Pl为岩芯样品的等效孔隙结构参数中的最小值，Pu为岩芯样品的等效孔隙结构参数中的最大值；将数值分布范围及分布规律、开始生成参数的信息发送给随机模拟等效孔隙结构参数生成单元；
随机模拟等效孔隙结构参数生成单元30：分别与所述统计单元20、判断单元60连接，在接收到统计单元发送来的开始生成参数的信息后，或者在接收到判断单元发送来的重新生成随机模拟等效孔隙结构参数的消息后，按照统计单元发送来的分布规律在数值分布范围[Pl,Pu]内随机生成一个随机模拟等效孔隙结构参数αs0；
速度计算单元40：分别与测井数据输入单元10、随机模拟等效孔隙结构参数生成单元30连接，利用随机模拟等效孔隙结构参数生成单元发送来的随机模拟等效孔隙结构参数αs0和测井数据输入单元发送来的测井数据，由岩石物理模型计算得到纵波速度Vp(αs0)、横波速度Vs(αs0)；
误差函数建立单元50：分别与测井数据输入单元10、速度计算单元40连接，建立误差函数f(αs0)；
判断单元60：分别与随机模拟等效孔隙结构参数生成单元30、误差函数建立单元50、反演结果输出单元70连接，判断f(αs0)≤ε是否成立，如果是，则发送输出结果的消息给反演结果输出单元，如果否，则发送重新生成随机模拟等效孔隙结构参数的消息给随机模拟等效孔隙结构参数生成单元；
反演结果输出单元70：分别与随机模拟等效孔隙结构参数生成单元30、判断单元60连接，在接收到判断单元发送来的输出结果的消息后，将随机模拟等效孔隙结构参数生成单元中最新的随机模拟等效孔隙结构参数作为反演结果进行输出。
下面以某地区的实际测井数据为例，利用本发明方法进行岩石孔隙结构参数反演，进而说明本发明的效果。
实施例1：
本实施例是本发明用于某地区的实际测井数据进行岩石孔隙结构参数反演的实施实例，方法的步骤如下：
(1)输入测井数据，包括纵波速度、横波速度、泥质含量、密度、含水饱和度、孔隙度等测井曲线。
(2)利用岩芯样品的等效孔隙结构参数确定等效孔隙结构参数的数值分布范围及分布规律。具体的，在目的层段内采集多个岩芯，获得这些岩芯的等效孔隙结构参数，找到这些参数中的最大值和最小值即得到数值分布范围，统计这些参数的分布规律，例如正态分布或者均匀分布，下面再利用这些岩芯的等效孔隙结构参数反演出该目的层段的等效孔隙结构参数。
(3)按照统计岩芯样品等效孔隙结构参数数值的分布规律在参数范围[Pl,Pu]内随机生成一个等效孔隙结构参数αs0，利用岩石物理模型计算等效孔隙结构参数为αs0时的纵波速度和横波速度。
(4)利用上述步骤(3)随机模拟等效孔隙结构参数，由岩石物理模型计算纵波速度、横波速度，将计算的纵波速度、横波速度与测井实测的纵波速度和横波速度进行比较；
(5)迭代进行上述步骤(3)和(4)，从中选出使得岩石物理模型计算结果与实测数据一致性最好的随机模拟的等效孔隙结构参数(即能够满足f(αs)≤ε的等效孔隙结构参数)，此时该参数就作为反演结果输出。
图1是实施例1中的输入测井数据，从左至右侧依次为：含水饱和度Sw、密度Density、孔隙度Porosity、纵波速度VP、横波速度VS、纵横波速度比、泥质含量Vsh(图1中的VCLGR表示泥质含量Vsh。图1中的JP3-6这些是地质分层号，其中有油的层即为目的层段)。图2是实施例1中反演的孔隙结构参数，图2中从左至右依次为含水饱和度、自然伽马、孔隙结构参数曲线，纵横波速度比，SW为含水饱和度曲线，GR_edit为自然伽马，ARS_x为孔隙结构参数，VPVS为纵横波速度比，其中第三个孔隙结构参数曲线是反演结果，图中a所示部分为含流体储层，具有低速度比高孔隙结构的特征。图3是实施例1中孔隙结构参数与自然伽马交会图，含流体部分具有低自然伽马高孔隙结构的特征。
本发明公开一种碎屑岩等效孔隙结构反演方法，属于岩石物理学及地球物理测井技术领域。本发明方法利用随机模拟和岩石物理模型结合从常规测井资料中提取岩石等效孔隙结构参数。本发明方法所需要的参数少，只需要纵波速度、横波速度、孔隙度、泥质含量、含水饱和度，就可以预测出岩石等效孔隙结构，可以避免传统的经验公式中因统计关系而存在的不确定性，应用范围广。具有计算速度快、稳定性好的优点，可以为工区内测井横波速度预测提供重要参数。
上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。
Claims (10)
1.一种碎屑岩等效孔隙结构反演方法，其特征在于：所述方法利用目的层段的多个岩芯样品的等效孔隙结构参数和测井数据中的纵波速度、横波速度、孔隙度、泥质含量、含水饱和度反演出该目的层段的等效孔隙结构参数。
2.根据权利要求1所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述方法包括：
(1)根据岩芯样品测试数据统计等效孔隙结构参数的数值分布范围及分布规律；
(2)利用所述数值分布范围及分布规律获得随机模拟等效孔隙结构参数，并计算得到纵波速度和横波速度；
(3)利用计算得到的纵波速度和横波速度以及测井数据中的纵波速度和横波速度建立误差函数f(αs0)；
(4)判断f(αs0)≤ε是否成立，如果是，则转入步骤(5)，如果否，则返回步骤(2)；
(5)将所述随机模拟等效孔隙结构参数作为等效孔隙结构参数的反演结果进行输出。
3.根据权利要求2所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(1)的操作包括：
输入测井数据，所述测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；
统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[Pl,Pu]，其中，Pl为岩芯样品的等效孔隙结构参数中的最小值，Pu为岩芯样品的等效孔隙结构参数中的最大值。
4.根据权利要求3所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(2)的操作包括：
按照所述分布规律在数值分布范围[Pl,Pu]内随机生成一个随机模拟等效孔隙结构参数αs0；
利用所述随机模拟等效孔隙结构参数αs0和孔隙度、泥质含量、含水饱和度，由岩石物理模型计算得到纵波速度Vp(αs0)、横波速度Vs(αs0)。
5.根据权利要求4所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述岩石物理模型采用的是适用于碎屑岩的岩石物理模型。
6.根据权利要求2所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(3)的操作包括：
利用公式(1)建立误差函数f(αs0)：
f(αs0)＝||Vp(αs0)-Vpm||+||Vs(αs0)-Vsm||+|Vp(αs0)/Vs(αs0)-Vpm/Vsm| (1)
公式(1)中，Vpm为测井数据中的纵波速度，Vsm为测井数据中的横波速度。
7.一种实现权利要求1-6任一项所述的碎屑岩等效孔隙结构反演方法的系统，其特征在于：所述系统包括：
测井数据输入单元：用于输入测井数据，测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；
统计单元：用于统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[Pl,Pu]，其中，Pl为岩芯样品的等效孔隙结构参数中的最小值，Pu为岩芯样品的等效孔隙结构参数中的最大值；将数值分布范围及分布规律、开始生成参数的信息发送给随机模拟等效孔隙结构参数生成单元；
随机模拟等效孔隙结构参数生成单元：分别与所述统计单元、判断单元连接，用于在接收到统计单元发送来的开始生成参数的信息后，或者在接收到判断单元发送来的重新生成随机模拟等效孔隙结构参数的消息后，按照统计单元发送来的分布规律在数值分布范围[Pl,Pu]内随机生成一个随机模拟等效孔隙结构参数αs0；
速度计算单元：分别与测井数据输入单元、随机模拟等效孔隙结构参数生成单元连接，用于利用随机模拟等效孔隙结构参数生成单元发送来的随机模拟等效孔隙结构参数αs0和测井数据输入单元发送来的测井数据，由岩石物理模型计算得到纵波速度Vp(αs0)、横波速度Vs(αs0)；
误差函数建立单元：分别与测井数据输入单元、速度计算单元连接，用于建立误差函数f(αs0)；
判断单元：分别与随机模拟等效孔隙结构参数生成单元、误差函数建立单元、反演结果输出单元连接，用于判断f(αs0)≤ε是否成立，如果是，则发送输出结果的消息给反演结果输出单元，如果否，则发送重新生成随机模拟等效孔隙结构参数的消息给随机模拟等效孔隙结构参数生成单元；
反演结果输出单元：分别与随机模拟等效孔隙结构参数生成单元、判断单元连接，用于在接收到判断单元发送来的输出结果的消息后，将随机模拟等效孔隙结构参数生成单元中最新的随机模拟等效孔隙结构参数作为反演结果进行输出。
8.根据权利要求7所述的系统，其特征在于：所述速度计算单元采用的岩石物理模型是适用于碎屑岩的岩石物理模型。
9.根据权利要求7所述的系统，其特征在于：所述误差函数建立单元利用公式(1)建立误差函数f(αs0)：
f(αs0)＝||Vp(αs0)-Vpm||+||Vs(αs0)-Vsm||+|Vp(αs0)/Vs(αs0)-Vpm/Vsm| (1)
公式(1)中，Vpm为测井数据中的纵波速度，Vsm为测井数据中的横波速度。
10.一种计算机可读存储介质，其特征在于：所述计算机可读存储介质存储有计算机可执行的至少一个程序，所述至少一个程序被所述计算机执行时使所述计算机执行权利要求1-6任一项所述碎屑岩等效孔隙结构反演方法中的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910908097.9A CN112558177A (zh) | 2019-09-25 | 2019-09-25 | 一种碎屑岩等效孔隙结构反演方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910908097.9A CN112558177A (zh) | 2019-09-25 | 2019-09-25 | 一种碎屑岩等效孔隙结构反演方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112558177A true CN112558177A (zh) | 2021-03-26 |
Family
ID=75029182
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910908097.9A Pending CN112558177A (zh) | 2019-09-25 | 2019-09-25 | 一种碎屑岩等效孔隙结构反演方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112558177A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060273788A1 (en) * | 2005-06-03 | 2006-12-07 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| US20080015782A1 (en) * | 2004-05-27 | 2008-01-17 | Saltzer Rebecca L | Method For Predicting Lithology And Porosity From Seismic Reflection Data |
| US20100326669A1 (en) * | 2008-04-09 | 2010-12-30 | Yaping Zhu | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
| CN104316962A (zh) * | 2014-07-28 | 2015-01-28 | 中国石油化工股份有限公司 | 岩石孔隙结构参数反演方法及其系统 |
| CN105223616A (zh) * | 2015-10-29 | 2016-01-06 | 中国石油大学(北京) | 一种页岩储层的孔隙纵横比反演方法 |
| CN105425280A (zh) * | 2015-11-21 | 2016-03-23 | 西南石油大学 | 一种矿物模量和孔隙结构预测方法 |
| CN108957580A (zh) * | 2017-05-22 | 2018-12-07 | 中国石油化工股份有限公司 | 一种反演页岩地层孔隙结构的方法 |
-
2019
- 2019-09-25 CN CN201910908097.9A patent/CN112558177A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080015782A1 (en) * | 2004-05-27 | 2008-01-17 | Saltzer Rebecca L | Method For Predicting Lithology And Porosity From Seismic Reflection Data |
| US20060273788A1 (en) * | 2005-06-03 | 2006-12-07 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| US20100326669A1 (en) * | 2008-04-09 | 2010-12-30 | Yaping Zhu | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
| CN104316962A (zh) * | 2014-07-28 | 2015-01-28 | 中国石油化工股份有限公司 | 岩石孔隙结构参数反演方法及其系统 |
| CN105223616A (zh) * | 2015-10-29 | 2016-01-06 | 中国石油大学(北京) | 一种页岩储层的孔隙纵横比反演方法 |
| CN105425280A (zh) * | 2015-11-21 | 2016-03-23 | 西南石油大学 | 一种矿物模量和孔隙结构预测方法 |
| CN108957580A (zh) * | 2017-05-22 | 2018-12-07 | 中国石油化工股份有限公司 | 一种反演页岩地层孔隙结构的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 叶庆全，冀宝发，王建新: "油气田开发地质", 31 May 1999, 北京：石油工业出版社, pages: 159 * |
| 周枫;徐鸣洁;马中高;: "等效孔隙结构模型在鄂尔多斯致密含气砂岩中的应用", 高校地质学报, vol. 19, no. 04, pages 588 - 5593 * |
| 李宏兵;张佳佳;姚逢昌;: "岩石的等效孔隙纵横比反演及其应用", 地球物理学报, vol. 56, no. 02, pages 608 - 615 * |
| 蔡涵鹏;贺振华;唐湘蓉;何光明;邹文;: "碳酸盐岩孔隙结构影响分析和等效孔隙结构参数计算", 石油物探, vol. 52, no. 06, pages 566 - 573 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210132246A1 (en) | Method for determining a grid cell size in geomechanical modeling of fractured reservoirs | |
| CN109655903B (zh) | 页岩层横波速度预测方法及系统 | |
| CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
| CN103198363B (zh) | 一种基于ct孔隙分析的储层产气量预测方法及装置 | |
| CN107917865A (zh) | 一种致密砂岩储层多参数渗透率预测方法 | |
| CN104533400A (zh) | 一种重构测井曲线的方法 | |
| Elmo et al. | Guidelines for the quantitative description of discontinuities for use in discrete fracture network modelling | |
| CN105134185A (zh) | 储层流体性质识别方法 | |
| CN105425280A (zh) | 一种矿物模量和孔隙结构预测方法 | |
| CN107808032A (zh) | 页岩气储层含水饱和度计算方法及系统 | |
| CN113138106A (zh) | 基于随钻岩屑录井资料的岩石弹性参数确定方法 | |
| CN110907324A (zh) | 一种碳酸盐岩孔隙成分分析方法及系统 | |
| CN107506556B (zh) | 一种确定新鲜完整岩块声波纵波速度值的方法 | |
| CN104316962B (zh) | 岩石孔隙结构参数反演方法及其系统 | |
| CN110794039B (zh) | 利用岩体波速计算帷幕灌浆岩体裂隙充填率的方法 | |
| CN110118994B (zh) | 一种基于地震反演和机器学习的陆相烃源岩定量预测方法 | |
| CN111751878A (zh) | 横波速度的预测方法和装置 | |
| CN112198551A (zh) | 一种通用的定量化评价储层裂缝强度的新方法 | |
| CN112558177A (zh) | 一种碎屑岩等效孔隙结构反演方法及系统 | |
| CN110954942A (zh) | 粘土含量识别的图版制作方法及系统 | |
| Hou et al. | Prediction of microscopic remaining oil distribution using fuzzy comprehensive evaluation | |
| CN114153002B (zh) | 储层天然裂缝三维地质建模方法、装置、电子设备及介质 | |
| CN108415074A (zh) | 一种缝洞储集体地震响应特征分析方法及装置 | |
| CN111830593A (zh) | 一种储层物性预测方法及装置 | |
| CN112558152B (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 |