CN106353813A - 基于阵列声波测井的流体性质识别方法 - Google Patents
基于阵列声波测井的流体性质识别方法 Download PDFInfo
- Publication number
- CN106353813A CN106353813A CN201610941279.2A CN201610941279A CN106353813A CN 106353813 A CN106353813 A CN 106353813A CN 201610941279 A CN201610941279 A CN 201610941279A CN 106353813 A CN106353813 A CN 106353813A
- Authority
- CN
- China
- Prior art keywords
- modulus
- rock
- matrix
- phi
- rock volume
- 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.)
- Granted
Links
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. analysis, for interpretation, for correction
- 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
技术领域
本发明属于石油勘探开发中储层流体性质测井评价技术领域,涉及常规及复杂储层流体性质定量判别方法,特别涉及利用阵列声波测井资料识别流体性质的方法。
背景技术
储层流体性质判别是测井解释评价中的主要内容之一,它直接影响到储量规模、地质认识,是勘探开发中储层评价的一个关键步骤。目前流体性质识别主要依赖基于电法测井的常规测井系列,利用声波测井识别流体性质主要有时差交会法、力学参数法。
基于电法测井的常规测井系列识别流体性质主要依赖储层流体导电能力大小,根据测量岩石电阻率的高低准确判别油气层与水层,利用孔隙度测井判别油、气层。通过测试、取样等手段获取油层电阻率下限值,当储层电阻率大于油层电阻率下限值时为油气层,否则含水。对于气层,中子孔隙度测井存在“挖掘效应”,导致中子、密度曲线出现“镜像”特征,区分油层与气层。
时差交会法根据阵列声波测井纵波时差与纵横波速比交会识别流体性质,时差交会法通过确定泥岩区与砂岩区对岩性进行分类,对于储层当饱含气时,导致纵波时差增大,而横波时差基本不受影响,纵横波速比降低,气层时差交会点远离砂岩区,进而识别气层。
力学参数法根据纵、横波时差及密度计算岩石力学参数,通过纵横波速度比与体积压缩系数进行流体性质识别,体积压缩系数公式如下:
Cb为体积压缩系数,GPa-1;
DTC、DTS分别为纵波时差、横波时差,μs/ft;
ρb为岩石体积密度,g/cm3。
对于气层,体积压缩系数降低、纵横波速比降低,两者交会可以较好的识别气层。
基于电法测井的常规测井系列对于普通砂泥岩储层具有较好效果,但对于低阻油层、复杂岩性储层,电阻率对油气层的响应特征受其他因素影响难以区分,因此对于这类储层常规测井系列识别流体性质面临较大困难。时差交会法与力学参数法对气层识别效果较好,而对于油层与水层,由于纵、横波时差差异较小,难以准确区分,且两种方法受岩性等其他因素影响较大,通常为定性识别,难以达到定量识别精度。
发明内容
为了解决现有技术在流体性质识别中的不足,本发明提供了一种基于阵列声波测井的流体性质识别方法,利用阵列声波测井资料,对提取的纵、横波时差进行岩性校正,提取变骨架参数值,通过计算模量因子消除岩性与孔隙特征影响,精确识别流体性质,达到复杂储层流体性质定量识别目的。
本发明的基于阵列声波测井的流体性质识别方法按照以下步骤进行操作:
步骤一:针对区域内已经能够确定储层流体性质的油井,获取该油井阵列声波测井的纵波时差值和横波时差值,计算岩石体积模量;
步骤二:对岩性校正,计算岩石变骨架时差值;
步骤三:根据岩石变骨架时差值,计算干岩石体积模量和饱含水岩石体积模量;
步骤四:根据岩石体积模量、干岩石体积模量和饱含水岩石体积模量,计算模量因子;
步骤五:通过统计该油井内不同性质的储层流体对应的模量因子,确定各储层流体的模量因子界限值,建立该区域内的储层流体模量因子交会图;
步骤六:对于该区域内待识别储层流体性质的油井,根据步骤二、步骤三和步骤四计算该油井中的待识别储层流体的模量因子;
步骤七:将步骤六中计算出的该油井中的待识别储层流体的模量因子与步骤五中建立的该区域内的储层流体模量因子交会图进行比较,以确定该油井中的待识别储层流体的性质。
其中,步骤一中计算岩石体积模量的公式为:
其中,K为岩石体积模量,GPa;
DTC、DTS分别为纵波时差值、横波时差值,μs/ft;
ρb为岩石体积密度,g/cm3。
步骤二中计算岩石变骨架时差值的方法为:对岩石骨架中两种主要矿物成分进行岩性校正,得到任意组合下的等效岩石骨架时差值,实现测井不同深度点不同矿物含量组合的变骨架时差值,公式为:
其中,DTCma、DTSma分别为纵波变骨架时差值、横波变骨架时差值,μs/ft;
DTC1、DTS1分别为矿物1的纵波时差值、横波时差值,μs/ft;
DTC2、DTS2分别为矿物2的纵波时差值、横波时差值,μs/ft;
v1、v2分别为矿物1、矿物2在岩石骨架中所占体积百分含量,%;
φ为岩石骨架中的孔隙度,%。
步骤三中计算干岩石体积模量和饱含水岩石体积模量的方法为:采用微分等效介质模型进行计算,计算过程采用四阶龙格-库塔算法迭代求解;
干岩石体积模量计算公式为:
饱含水岩石体积模量计算公式为:
其中,Kdry、Gdry分别为干岩石等效体积模量与干岩石等效剪切模量,GPa;Ksatw、Gsatw分别为饱含水岩石等效体积模量与饱含水岩石等效剪切模量,GPa;
Km、Gm分别为岩石骨架体积模量与岩石骨架剪切模量,GPa;
αi为岩石骨架中的孔隙纵横比,无量纲;
vi为岩石骨架中的孔隙纵横比对应的孔隙百分含量,无量纲;
βm为骨架参数,表达式为βm=Gm(3Km+Gm)/(3Km+4Gm);
φ为岩石骨架中的孔隙度,%。
步骤四中计算模量因子的公式为:
其中,M为模量因子,无量纲;
K为岩石体积模量,GPa;
Kdry为干岩石体积模量,GPa;
Ksatw为饱含水岩石体积模量,GPa。
相对于现有技术,本发明的有益效果为:本发明提供了一种基于阵列声波测井的流体性质识别方法,该方法依据阵列声波测井资料,对纵、横波时差进行岩性校正,通过变骨架参数值计算干岩石、饱含水岩石模量,使得模量因子只反映流体性质影响,进而实现流体性质定量识别。本发明所述方法消除岩性影响,模量因子不受储层岩性、孔隙特征限制,对常规储层、复杂储层均适用,为流体性质识别提供一种新的技术方法,具有推广价值。
下面结合附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
为了更清楚的说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图做简单介绍,显而易见,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些图获得其它的附图。
图1为本发明的基于阵列声波测井的流体性质识别方法的流程图;
图2为本发明实施例提供的某油田岩性校正前的时差交会图;
图3为本发明实施例提供的某油田岩性校正后的时差交会图;
图4为本发明实施例提供的某区块三模量参数交会图;
图5为本发明实施例提供的模量因子与纵波时差交会图;
图6为本发明实施例提供的模量因子与纵横波速比交会图;
图7为本发明实施例提供的模量因子与流体模量交会图;
图8为本发明实施例提供的某井模量因子判别流体性质效果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明实施方式做进一步详细描述。
参见图1,本发明的基于阵列声波测井的流体性质识别方法,按照如下步骤进行操作:
步骤一:获取区域内阵列声波测井纵、横波时差值,计算岩石体积模量。
根据区域内经测试、取样等手段证实流体性质的储层,或根据测井资料能够判别流体性质的不同类型储层,对阵列声波测井进行处理,提取纵、横波时差值,选取具有代表性的油、气、水层声波时差,计算岩石体积模量,公式如下:
K为地层体积模量,GPa;
DTC、DTS分别为地层纵、横波时差值,μs/ft;
ρb为地层岩石体积密度,g/cm3。
步骤二:岩性校正,计算岩石变骨架时差值。
理论上纵波时差对流体性质具有敏感性,通常采用纵波时差与纵横波速比交会识别流体性质,由于储层骨架矿物组分不同,不同矿物成分组合使得储层骨架声波时差值差异较大,因此受矿物组分影响传统的时差交会法识别流体性质存在不确定性,如图2所示,部分油层、水层存在重叠现象,导致流体性质不易区分。
岩性校正针对岩石骨架中两种主要矿物成分进行校正,得到任意矿物成分组合下的等效岩石骨架时差值,实现测井不同深度点不同矿物含量组合时的变骨架时差值,公式如下:
DTCma、DTSma分别为纵、横波变骨架时差值,μs/ft;
DTC1、DTS1分别为矿物1的纵、横波时差值,μs/ft;
DTC2、DTS2分别为矿物2的纵、横波时差值,μs/ft;
v1、v2分别为矿物1、矿物2所占体积百分含量,%;
φ为孔隙度,%。
通过岩性校正,消除岩性组分差异对纵、横波时差的影响,将纵、横波时差值校正到单一骨架背景下,突显流体性质影响效果,如图3所示,通过岩性校正,相比校正前油层、气层、水层区分明显。
步骤三:计算干岩石体积模量和饱含水岩石体积模量。
理论与实际处理发现孔隙对阵列声波时差值影响较大,因此需要消除孔隙影响,通常采用计算干岩石的方法消除孔隙影响。干岩石体积模量采用微分等效介质模型进行计算,计算过程采用四阶龙格-库塔(Runge-Kutta)算法迭代求解,其计算公式如下:
Kdry、Gdry分别为干岩石等效体积模量与剪切模量,GPa;
Km、Gm分别为岩石骨架体积模量与剪切模量,GPa;
αi为孔隙纵横比,无量纲;
vi为孔隙纵横比对应的孔隙百分含量,无量纲;
βm为骨架参数,表达式为βm=Gm(3Km+Gm)/(3Km+4Gm)。
由于干岩石体积模量不包含流体性质的影响,因此干岩石体积模量作为储层体积模量的下限值,当储层饱含水时,岩石体积模量达到上限。饱含水岩石体积模量采用微分等效介质模型进行计算,计算过程采用四阶龙格-库塔算法迭代求解,其计算公式如下:
Ksatw、Gsatw分别为饱含水岩石等效体积模量与干岩石等效剪切模量,GPa。
阵列声波测井计算的岩石体积模量介于干岩石体积模量与饱含水岩石体积模量之间,由于计算过程中消除岩性、孔隙的影响,因此通过干岩石体积模量、饱含水岩石体积模量与储层岩石体积模量进行交会能够较好的反映流体性质的影响,如图4所示,三模量参数交会对于油层、气层、水层区分效果较好。
步骤四:计算模量因子。
由于不同储层基质岩性差异较大,导致储层等效骨架声波时差值差异较大,直接采用模量交会方法计算的绝对模量值在同一个区域不同岩性油田差异较大,无法应用于区域内流体性质定量判别,因此通过定义模量因子消除背景岩性影响,定量识别流体性质,模量因子采用下式计算
M为模量因子,无量纲;
Kdry为干岩石体积模量,GPa;
Ksatw为饱含水岩石体积模量,GPa。
步骤五:建立该区域内的储层流体模量因子交会图。
模量因子采用相对值消除了不同岩性骨架值的影响,适用于整个区域流体性质判别,表1为渤海油田某区块不同类型储层油、气、水层模量因子计算结果,通过统计不同流体性质模量因子范围确定该区块气层模量因子上限为0.56,水层模量因子下限为0.85,油层模量因子范围为0.56~0.85。
表1
根据模量因子确定的储层流体界限值,建立区域内模量因子交会图,如图5、图6、图7所示。图5为模量因子与纵波时差交会图,图6为模量因子与纵横波速比交会图,图7为模量因子与流体模量交会图,从图中可以看出,尽管纵波时差、纵横波速比、流体模量在油、气、水层存在一定不确定性,但模量因子区分明显,因此通过交会图便能够有效判别流体性质。
步骤六:对于该区域内待识别储层流体性质的油井,根据步骤二、步骤三和步骤四计算该油井中的待识别储层流体的模量因子。
步骤七:将计算出的该油井中的待识别储层流体的模量因子与建立的该区域内的储层流体模量因子交会图进行比较,以确定该油井中的待识别储层流体的性质。
根据步骤五中确定的不同流体性质对应的模量因子界限,对该区块某井进行处理,处理结果如图8所示。图中第一道为CAL(井径曲线)、GR(自然伽马曲线)和SP(自然电位曲线);第二道为深度;第三道为RD(深电阻率曲线)、RS(浅电阻率曲线)、RMLL(微侧向电阻率曲线);第四道为(ZDEN密度曲线)、CN(中子曲线);第五道中M为利用本发明计算的模量因子曲线;第六道为解释结论。根据步骤五中确定的气层模量因子界限(0.56)与水层模量因子界限(0.85),将模量因子划分为三个区域,如图第五道所示,根据计算模量因子的大小,对其进行流体性质识别,其中1、2、6、7、8号层模量因子小于0.56,解释为气层,9号层解释为油层,3、4、5、10号层解释为水层。结合常规测井可以看出,1、2、7号层常规测井GR值较低,SP正异常且较低,显示储层,ZDEN、CN出现“镜像”交会特征,电阻率较高,显示气层特征,与模量因子解释一致。9号层RD较高,SP正异常且较低,ZDEN、CN无“镜像”交会特征,显示油层特征,与模量因子解释一致。3、4、5、10号层RD较低,SP正异常且较高,ZDEN、CN存在交会,显示水层特征,与模量因子解释一致。6、8号层RD较低,GR较高,SP正开常且较低,ZDEN、CN交会不明显,为低阻储层,常规测井显示气层特征不明显,与模量因子存在一定差异,对8、9号层进行测试发现日产油49.6方,日产气66441方,证实8号层为气层,从而证实模量因子的可靠性,适用于常规及复杂储层流体性质识别当中。
以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改、等同替换、改进,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于阵列声波测井的流体性质识别方法,其特征在于,包括以下步骤:
步骤一:针对区域内已经能够确定储层流体性质的油井,获取该油井阵列声波测井的纵波时差值和横波时差值,计算岩石体积模量;
步骤二:对岩性校正,计算岩石变骨架时差值;
步骤三:根据岩石变骨架时差值,计算干岩石体积模量和饱含水岩石体积模量;
步骤四:根据岩石体积模量、干岩石体积模量和饱含水岩石体积模量,计算模量因子;
步骤五:通过统计该油井内不同性质的储层流体对应的模量因子,确定各储层流体的模量因子界限值,建立该区域内的储层流体模量因子交会图;
步骤六:对于该区域内待识别储层流体性质的油井,根据步骤二、步骤三和步骤四计算该油井中的待识别储层流体的模量因子;
步骤七:将步骤六中计算出的该油井中的待识别储层流体的模量因子与步骤五中建立的该区域内的储层流体模量因子交会图进行比较,以确定该油井中的待识别储层流体的性质。
2.根据权利要求1所述的基于阵列声波测井的流体性质识别方法,其特征在于,步骤一中计算岩石体积模量的公式为:
其中,K为岩石体积模量,GPa;
DTC、DTS分别为纵波时差值、横波时差值,μs/ft;
ρb为岩石体积密度,g/cm3。
3.根据权利要求2所述的基于阵列声波测井的流体性质识别方法,其特征在于,步骤二中计算岩石变骨架时差值的方法为:对岩石骨架中两种主要矿物成分进行岩性校正,得到任意组合下的等效岩石骨架时差值,实现测井不同深度点不同矿物含量组合的变骨架时差值,公式为:
其中,DTCma、DTSma分别为纵波变骨架时差值、横波变骨架时差值,μs/ft;
DTC1、DTS1分别为矿物1的纵波时差值、横波时差值,μs/ft;
DTC2、DTS2分别为矿物2的纵波时差值、横波时差值,μs/ft;
v1、v2分别为矿物1、矿物2在岩石骨架中所占体积百分含量,%;
φ为岩石骨架中的孔隙度,%。
4.根据权利要求3所述的基于阵列声波测井的流体性质识别方法,其特征在于,步骤三中计算干岩石体积模量和饱含水岩石体积模量的方法为:采用微分等效介质模型进行计算,计算过程采用四阶龙格-库塔算法迭代求解;
干岩石体积模量计算公式为:
饱含水岩石体积模量计算公式为:
其中,Kdry、Gdry分别为干岩石等效体积模量与干岩石等效剪切模量,GPa;
Ksatw、Gsatw分别为饱含水岩石等效体积模量与饱含水岩石等效剪切模量,GPa;
Km、Gm分别为岩石骨架体积模量与岩石骨架剪切模量,GPa;
αi为岩石骨架中的孔隙纵横比,无量纲;
vi为岩石骨架中的孔隙纵横比对应的孔隙百分含量,无量纲;
βm为骨架参数,表达式为βm=Gm(3Km+Gm)/(3Km+4Gm);
φ为岩石骨架中的孔隙度,%。
5.根据权利要求4所述的基于阵列声波测井的流体性质识别方法,其特征在于,步骤四中计算模量因子的公式为:
其中,M为模量因子,无量纲;
K为岩石体积模量,GPa;
Kdry为干岩石体积模量,GPa;
Ksatw为饱含水岩石体积模量,GPa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610941279.2A CN106353813B (zh) | 2016-10-25 | 2016-10-25 | 基于阵列声波测井的流体性质识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610941279.2A CN106353813B (zh) | 2016-10-25 | 2016-10-25 | 基于阵列声波测井的流体性质识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106353813A true CN106353813A (zh) | 2017-01-25 |
CN106353813B CN106353813B (zh) | 2018-12-25 |
Family
ID=57864225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610941279.2A Active CN106353813B (zh) | 2016-10-25 | 2016-10-25 | 基于阵列声波测井的流体性质识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106353813B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107907910A (zh) * | 2017-11-16 | 2018-04-13 | 山东胜软科技股份有限公司 | 一种不同岩性油藏横波测井确定方法 |
CN110005396A (zh) * | 2018-12-11 | 2019-07-12 | 中国石油天然气集团有限公司 | 一种测井波列数据深度弹性校正方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887132A (zh) * | 2009-05-15 | 2010-11-17 | 中国石油天然气股份有限公司 | 一种井震结合定量预测砂岩储层流体饱和度的方法 |
CN102096107A (zh) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 |
CN102426390A (zh) * | 2011-10-21 | 2012-04-25 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN103513270A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | 一种基于岩石声学特性的气层识别评价方法及装置 |
CN103792575A (zh) * | 2012-11-01 | 2014-05-14 | 中国石油天然气集团公司 | 一种去骨架影响的致密砂岩气层识别评价方法 |
CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
CN103867197A (zh) * | 2014-04-04 | 2014-06-18 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性天然气层声波时差判别法 |
CN104345339A (zh) * | 2013-07-25 | 2015-02-11 | 中国石油天然气集团公司 | 一种用阵列声波测井资料确定岩石脆性系数的方法 |
CN104360389A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种致密砂岩储层岩石弹性模量计算方法 |
CN104570103A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种低信噪比地震资料的井约束速度谱拾取方法 |
US9176243B2 (en) * | 2006-02-24 | 2015-11-03 | Hanner Zueroher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
CN105589110A (zh) * | 2015-12-10 | 2016-05-18 | 中国石油集团川庆钻探工程有限公司 | 孔隙流体压缩系数判别致密砂岩气层的方法 |
CN105988136A (zh) * | 2015-01-30 | 2016-10-05 | 中石化石油工程技术服务有限公司 | 一种利用纵横波速度信息进行含气性分析的方法 |
-
2016
- 2016-10-25 CN CN201610941279.2A patent/CN106353813B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9176243B2 (en) * | 2006-02-24 | 2015-11-03 | Hanner Zueroher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
CN101887132A (zh) * | 2009-05-15 | 2010-11-17 | 中国石油天然气股份有限公司 | 一种井震结合定量预测砂岩储层流体饱和度的方法 |
CN102096107A (zh) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 |
CN102426390A (zh) * | 2011-10-21 | 2012-04-25 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN103513270A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | 一种基于岩石声学特性的气层识别评价方法及装置 |
CN103792575A (zh) * | 2012-11-01 | 2014-05-14 | 中国石油天然气集团公司 | 一种去骨架影响的致密砂岩气层识别评价方法 |
CN104345339A (zh) * | 2013-07-25 | 2015-02-11 | 中国石油天然气集团公司 | 一种用阵列声波测井资料确定岩石脆性系数的方法 |
CN104570103A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种低信噪比地震资料的井约束速度谱拾取方法 |
CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
CN103867197A (zh) * | 2014-04-04 | 2014-06-18 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性天然气层声波时差判别法 |
CN104360389A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种致密砂岩储层岩石弹性模量计算方法 |
CN105988136A (zh) * | 2015-01-30 | 2016-10-05 | 中石化石油工程技术服务有限公司 | 一种利用纵横波速度信息进行含气性分析的方法 |
CN105589110A (zh) * | 2015-12-10 | 2016-05-18 | 中国石油集团川庆钻探工程有限公司 | 孔隙流体压缩系数判别致密砂岩气层的方法 |
Non-Patent Citations (2)
Title |
---|
尹帅 等: ""基于阵列声波测井的海陆过渡相碎屑岩地层裂缝识别方法"", 《石油钻探技术》 * |
王贵清 等: ""基于变骨架时差的纵横波速度比识别轻质油气层的方法研究"", 《测井技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107907910A (zh) * | 2017-11-16 | 2018-04-13 | 山东胜软科技股份有限公司 | 一种不同岩性油藏横波测井确定方法 |
CN107907910B (zh) * | 2017-11-16 | 2019-12-13 | 山东胜软科技股份有限公司 | 一种不同岩性油藏横波测井确定方法 |
CN110005396A (zh) * | 2018-12-11 | 2019-07-12 | 中国石油天然气集团有限公司 | 一种测井波列数据深度弹性校正方法 |
CN110005396B (zh) * | 2018-12-11 | 2023-02-10 | 中国石油天然气集团有限公司 | 一种测井波列数据深度弹性校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106353813B (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106468172B (zh) | 一种超低渗砂岩油藏低阻储层测井解释方法 | |
Shipton et al. | Structural heterogeneity and permeability in faulted eolian sandstone: Implications for subsurface modeling of faults | |
Worthington | Net pay—what is it? What does it do? How do we quantify it? How do we use it? | |
Nabawy et al. | Reservoir assessment and quality discrimination of Kareem Formation using integrated petrophysical data, Southern Gulf of Suez, Egypt | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN104278991A (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN110058323A (zh) | 一种致密砂岩地层脆性指数计算方法 | |
Infante-Paez et al. | Seismic analysis of paleotopography and stratigraphic controls on total organic carbon: Rich sweet spot distribution in the Woodford Shale, Oklahoma, USA | |
CN104714252A (zh) | 分析流体因子敏感性的方法 | |
CN107817260A (zh) | 一种泥页岩高频层序识别方法 | |
CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
CN106405639A (zh) | 一种叠前地震储层岩性参数的反演方法 | |
CN107229076B (zh) | 一种基于测井资料进行温度响应特征分析的方法 | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
Richardson | Well correlation and Petrophysical analysis, a case study of “Rickie” field onshore Niger Delta | |
CN106033127B (zh) | 基于横波速度变化率的地应力方位地震预测方法 | |
CN115857047B (zh) | 一种地震储层综合预测方法 | |
CN106223940A (zh) | 一种多层砂岩油藏低阻油层综合识别方法及装置 | |
CN111650644B (zh) | 盐胶结砂岩定量预测方法 | |
Slatt | Characterization of Unconventional Resource Shales (Mudstones) The Necessity of Multiscale Scientific Integration | |
Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
Davis | Velocity variations around Leduc reefs, Alberta | |
Ekine et al. | Delineation of hydrocarbon bearing reservoirs from surface seismic and well log data (Nembe Creek) in Niger Delta oil field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |