CN105003258B - 一种高温高压气层甲烷流体密度骨架参数的获取方法 - Google Patents
一种高温高压气层甲烷流体密度骨架参数的获取方法 Download PDFInfo
- Publication number
- CN105003258B CN105003258B CN201510484768.5A CN201510484768A CN105003258B CN 105003258 B CN105003258 B CN 105003258B CN 201510484768 A CN201510484768 A CN 201510484768A CN 105003258 B CN105003258 B CN 105003258B
- Authority
- CN
- China
- Prior art keywords
- density
- methane stream
- parameter
- matrix parameter
- bearing formation
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种高温高压气层甲烷流体密度骨架参数的获取方法,包括如下步骤:甲烷流体采样;获取多组测试实验数据;构建初级计算模型:利用数据拟合分析方法,对上述多组数据进行参数拟合分析,得到初级计算模型;构建密度骨架参数计算模型:在所述的初级计算模型的基础上,利用密度测井响应原理,将甲烷流体的体积密度转化为视体积密度值,得到密度骨架参数计算模型;将测得的温度值及压力值,代入到密度骨架参数计算模型,得到高温高压气层中甲烷流体的密度骨架参数。利用该方法可以显著提高获取甲烷流体骨架密度参数的准确度,进而提高高温高压气层甲烷气层测井解释评价孔隙度的精度,具有较强的通用性。
Description
技术领域
本发明涉及油气勘探技术领域,具体地说是一种高温高压气层甲烷流体密度骨架参数的获取方法。
背景技术
在高温高压甲烷气层测井解释评价过程中,需要测井技术人员对获得的井下密度测井资料依据地层组分体积模型进行及时的测井解释,获取气层的孔隙度参数,工作的关键在于确定地层条件下甲烷流体密度骨架参数。现有的甲烷流体密度骨架参数的获取方法多为区域经验法,该方法主要是利用井下取到的甲烷气样进行PVT实验,然后得到一定温度压力下的密度值,将其作为区域的骨架参数经验值,然而在长时间的生产实践中技术人员发现上述方法获取的甲烷流体密度骨架参数的准确性较低,究其原因主要高温高压储层其压力系数变化范围较大,以及井间高温高压甲烷气层在井上埋深范围也不一样,导致气层的地层压力、温度值不一致;另外该方法的主观性较大,用到了一些估算的过程;从而使得甲烷流体密度骨架参数经验值与实际地层中的值相差较大。
发明内容
针对以上存在的问题,本发明提出了一种高温高压气层甲烷流体密度骨架参数的获取方法,可以显著提高获取甲烷流体骨架密度参数的准确度,进而提高高温高压甲烷气层测井解释评价孔隙度的精度,具有较强的通用性。
为解决上述技术问题,本申请提供了一种高温高压气层甲烷流体密度骨架参数的获取方法,其特征在于,包括如下步骤:
S1、甲烷流体采样:利用取样仪器从高温高压气田井场选取甲烷流体样品;
S2、获取多组测试实验数据:在不同温度及压力条件下,对甲烷流体样品进行体积密度的测定试验,得到由体积密度参数、温度参数和压力参数构成的多组实验数据;
S3、构建初级计算模型:利用数学统计范畴中的数据拟合分析方法,对步骤S2的多组数据进行参数拟合分析,得到甲烷流体的体积密度随温度及压力变化的初级计算模型;
S4、构建密度骨架参数计算模型:在所述的初级计算模型的基础上,利用密度测井响应原理,将甲烷流体的体积密度转化为视体积密度值,得到高温高压气层的甲烷流体的密度骨架参数随温度及压力变化的密度骨架参数计算模型;
S5、获取密度骨架参数:对高温高压气层中的温度及压力进行实际检测,将测得的温度值及压力值,代入到密度骨架参数计算模型中,得到高温高压储层中甲烷流体的密度骨架参数。
所述步骤S2的测定试验中所用的设备为流体体积密度测量仪。
所述步骤S2的实验数据为20组,其中,温度的选取值为20℃、50℃、100℃及150℃,压力的选取值为20MPa、30MPa、40MPa、50MPa和60MPa。
所述的步骤S3中的初级计算模型具体如下:Y=(0.081487ln(p)-0.102241)e(-0.001366ln(p)-0.001565)T,
其中,Y为实验测量得到的甲烷流体体积密度,单位为(g/cm3),P为地层压力值,单位为(MPa),T为地层温度,单位为(℃)。
所述的S4中的密度骨架参数计算模型具体如下:ρa(ch4)=1.3366((0.081487ln(p)-0.102241)e(-0.001366ln(p)-0.001565)T)-0.1883其中,ρa(ch4)为高温高压气层甲烷流体密度骨架参数,单位为(g/cm3),P为地层压力值,单位为(MPa),T为地层温度,单位为(℃)。
有益效果
本申请具有如下技术效果或优点:
利用该方法可以显著提高获取甲烷流体骨架密度参数的准确度,进而提高高温高压甲烷气层测井解释评价孔隙度的精度,具有较强的通用性;本发明为高温高压气田测井解释评价中甲烷流体骨架密度参数选取提供了更好更快捷的获取方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例的方法流程图;
图2为本申请实施例的多组测试实验数据分布图。
具体实施方式
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
本发明提供了一种高温高压气层甲烷流体密度骨架参数的获取方法,可以显著提高获取甲烷流体骨架密度参数的准确度,进而提高高温高压甲烷气层测井解释评价孔隙度的精度,具有较强的通用性。
如图1所示,本发明所述的一种高温高压气层甲烷流体密度骨架的获取方法,包括如下步骤:
步骤1,甲烷流体采样:甲烷流体样品选取:利用取样仪器从高温高压气田井场采集得到甲烷流体样品;
步骤2,获取多组测试实验数据:对甲烷流体样品进行实验室流体体积密度试验分析,分别在温度为20℃、50℃、100℃及150℃,压力为20MPa、30MPa、40MPa、50MPa与60MPa条件下测定甲烷流体体积密度,共测量得到20组甲烷流体体积密度实验数据,将这20组甲烷流体体积密度实验数据以图的方式进行体现,得到甲烷流体体积密度在不同温度、压力条件下的数据散点图,如图2所示。经实验所得到的数据散点图可知,甲烷流体体积密度与温度、压力呈现的是二元相关关系,在同一温度条件下,随着压力的增大,甲烷流体体积密度增大,而在同一压力条件下,随着温度的升高,甲烷流体体积密度降低,也就是说甲烷流体体积密度在不同温度压力条件下具有不同的数值。
步骤3,构建初级计算模型:利用数学统计范畴中的数据拟合分析方法,对上述多组数据进行参数拟合分析,得到甲烷流体的体积密度随温度及压力变化的初级计算模型。
步骤4,构建密度骨架参数计算模型:在所述的初级计算模型的基础上,利用密度测井响应原理,将甲烷流体的体积密度转化为视体积密度值,得到高温高压储层的甲烷流体的密度骨架参数随温度及压力变化的密度骨架参数计算模型。
其具体实施步骤如下:由密度测井响应原理可知,密度测井测量的是地层的电子密度ρe,用含水石灰岩刻度表示的刻度方程为:
ρa=1.07(ρe)i-0.1883 (1)
其中(ρe)i为电子密度指数,定义为:
式中,N为阿伏伽德罗常数;
经过刻度方程刻度后,密度测井仪器记录的是地层的视体积密度ρa,地层的体积密度与电子密度的关系:
式中,N为阿伏伽德罗常数,Z为原子序数,A为原子量;
将(3)和(2)代入式(1),对于多原子组成的化合物,其刻度后的视体积密度为:
二氧化碳的分子式为CH4,将C、H的原子系数及原子量代入式(4),且将步骤3所得的甲烷流体体积密度与温度、压力参数计算模型代入式(4)可得到地层条件下甲烷的视体积密度数学表达式为:
ρa(ch4)=1.3366((0.081487ln(p)-0.102241)e(-0.001366ln(p)-0.001565)T)-0.1883
(5)
其中,ρa(ch4)为高温高压气层甲烷流体密度骨架参数,单位为(g/cm3)。
步骤5,获取密度骨架参数:对高温高压储层中的温度及压力进行实际检测,将测得的温度值及压力值,代入到密度骨架参数计算模型中,得到高温高压气层中甲烷流体的密度骨架参数。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例演示如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (1)
1.一种高温高压气层甲烷流体密度骨架参数的获取方法,其特征在于,包括如下步骤:
S1、甲烷流体采样:利用取样仪器从高温高压气田井场选取甲烷流体样品;
S2、获取多组测试实验数据:在不同温度及压力条件下,采用流体体积密度测量仪对甲烷流体样品进行体积密度的测定试验,得到由体积密度参数、温度参数和压力参数构成的实验数据20组,其中,温度的选取值为20℃、50℃、100℃及150℃,压力的选取值为20MPa、30MPa、40MPa、50MPa和60MPa;
S3、构建初级计算模型:利用数学统计范畴中的数据拟合分析方法,对步骤S2的多组数据进行参数拟合分析,得到甲烷流体的体积密度随温度及压力变化的初级计算模型,具体如下:Y=(0.081487ln(p)-0.102241)e(-0.001366ln(p)-0.001565)T,
其中,Y为实验测量得到的甲烷流体体积密度,单位为(g/cm3),P为地层压力值,单位为(MPa),T为地层温度,单位为(℃);
S4、构建密度骨架参数计算模型:在所述的初级计算模型的基础上,利用密度测井响应原理,将甲烷流体的体积密度转化为视体积密度值,得到高温高压气层的甲烷流体的密度骨架参数随温度及压力变化的密度骨架参数计算模型,具体如下:ρa(ch4)=1.3366((0.081487ln(p)-0.102241)e(-0.001366ln(p)-0.001565)T)-0.1883其中,ρa(ch4)为高温高压气层甲烷流体密度骨架参数,单位为(g/cm3),P为地层压力值,单位为(MPa),T为地层温度,单位为(℃);
S5、获取密度骨架参数:对高温高压气层中的温度及压力进行实际检测,将测得的温度值及压力值,代入到密度骨架参数计算模型中,得到高温高压储层中甲烷流体的密度骨架参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510484768.5A CN105003258B (zh) | 2015-08-07 | 2015-08-07 | 一种高温高压气层甲烷流体密度骨架参数的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510484768.5A CN105003258B (zh) | 2015-08-07 | 2015-08-07 | 一种高温高压气层甲烷流体密度骨架参数的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105003258A CN105003258A (zh) | 2015-10-28 |
CN105003258B true CN105003258B (zh) | 2017-03-15 |
Family
ID=54376089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510484768.5A Active CN105003258B (zh) | 2015-08-07 | 2015-08-07 | 一种高温高压气层甲烷流体密度骨架参数的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105003258B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105952446B (zh) * | 2016-04-26 | 2019-03-01 | 中国海洋石油集团有限公司 | 一种石油天然气气层成分含量测定方法 |
CN107605471B (zh) * | 2017-07-21 | 2020-06-09 | 中国石油天然气股份有限公司 | 一种湖相碳酸盐岩岩性测井定量识别方法 |
CN109269962B (zh) * | 2018-12-04 | 2024-02-13 | 郑州大学 | 一种超临界二氧化碳高温高压pvt测试及驱替甲烷一体化的实验装置及方法 |
CN117310812A (zh) * | 2023-09-27 | 2023-12-29 | 广东海洋大学 | 一种甲烷流体纵波时差骨架参数获取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1193488A1 (de) * | 2000-09-29 | 2002-04-03 | Ruhrgas Aktiengesellschaft | Verfahren und Vorrichtung zum Ermitteln der Gasbeschaffenheit eines Erdgases |
CN1621805A (zh) * | 2003-11-27 | 2005-06-01 | 赵万山 | 数显式液化石油气密度计装置 |
CN202110114U (zh) * | 2011-05-11 | 2012-01-11 | 江汉大学 | 一种测量空气密度的装置 |
JP5294110B2 (ja) * | 2008-07-07 | 2013-09-18 | 清水建設株式会社 | メタンハイドレートからのメタンガス生産方法及びメタンハイドレートからのメタンガス生産装置 |
CN103867184A (zh) * | 2014-02-10 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种气井临界携液流量确定方法及装置 |
-
2015
- 2015-08-07 CN CN201510484768.5A patent/CN105003258B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1193488A1 (de) * | 2000-09-29 | 2002-04-03 | Ruhrgas Aktiengesellschaft | Verfahren und Vorrichtung zum Ermitteln der Gasbeschaffenheit eines Erdgases |
CN1621805A (zh) * | 2003-11-27 | 2005-06-01 | 赵万山 | 数显式液化石油气密度计装置 |
JP5294110B2 (ja) * | 2008-07-07 | 2013-09-18 | 清水建設株式会社 | メタンハイドレートからのメタンガス生産方法及びメタンハイドレートからのメタンガス生産装置 |
CN202110114U (zh) * | 2011-05-11 | 2012-01-11 | 江汉大学 | 一种测量空气密度的装置 |
CN103867184A (zh) * | 2014-02-10 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种气井临界携液流量确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105003258A (zh) | 2015-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN103983746B (zh) | 一种富火山物质烃源岩有机碳恢复方法 | |
CN103454198B (zh) | 一种泥页岩有机孔隙度检测方法 | |
CN105003258B (zh) | 一种高温高压气层甲烷流体密度骨架参数的获取方法 | |
CN102220865B (zh) | 一种检测灰岩地层孔隙压力的方法 | |
CN104298883A (zh) | 一种油气资源评价中的烃源岩产烃率图版的建立方法 | |
Wiese et al. | Well-based hydraulic and geochemical monitoring of the above zone of the CO 2 reservoir at Ketzin, Germany | |
CN105842416B (zh) | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 | |
CN105930604B (zh) | 一种确定页岩气开采程度的方法及装置 | |
CN101082277A (zh) | 石油钻井地质x射线荧光岩屑录井方法 | |
MX344149B (es) | Metodo para determinar el volumen de materia organica en roca deposito. | |
CN108508182B (zh) | 快速确定笔石相热页岩中生物硅含量的测录井方法 | |
CN105952446B (zh) | 一种石油天然气气层成分含量测定方法 | |
CN101581222B (zh) | 一种石油钻井液核磁共振随钻分析方法 | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
Kurylyk et al. | Heat: An overlooked tool in the practicing hydrogeologist's toolbox | |
CN206458460U (zh) | 井下测量装置以及井下测量系统 | |
CN105651912A (zh) | 岩石热解分析仪及热解分析方法 | |
CN206095862U (zh) | 页岩气现场解吸测试装置 | |
CN106248548B (zh) | 包气带钻孔常水头注水试验饱和渗透系数测试方法 | |
CN105931125A (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
CN103590828B (zh) | 一种录井dck指数法评价地层压力的方法 | |
CN109458176A (zh) | 碳酸盐岩储层压力的预测方法及其应用 | |
WO2012088732A1 (zh) | 梯度地球化学勘探方法 | |
CN112145165A (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee after: CNOOC (China) Limited Zhanjiang Branch Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 China oil tower, 25 Chaoyangmen North Street, Dongcheng District, Beijing Co-patentee before: CNOOC (China) Limited Zhanjiang Branch Patentee before: China National Offshore Oil Corporation |
|
CP03 | Change of name, title or address |