CN105317429A - 一种确定页岩地层干酪根含量的方法 - Google Patents
一种确定页岩地层干酪根含量的方法 Download PDFInfo
- Publication number
- CN105317429A CN105317429A CN201410319217.9A CN201410319217A CN105317429A CN 105317429 A CN105317429 A CN 105317429A CN 201410319217 A CN201410319217 A CN 201410319217A CN 105317429 A CN105317429 A CN 105317429A
- Authority
- CN
- China
- Prior art keywords
- log
- mineral
- error
- curve
- kerogen
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种确定页岩地层干酪根含量的方法。包括:建立由骨架矿物、干酪根和孔隙组成的页岩地层体积模型,其中孔隙包括气体和水;基于所述体积模型构造页岩地层在不同体积含量参数条件下的测井响应方程,确定响应方程误差和测量误差;采集实际测井曲线,以页岩地层矿物质量含量为约束条件,根据所述测井响应方程,响应方程误差和测量误差建立测井解释的目标函数;计算使所述目标函数取得最小值的最优体积含量参数,绘制最优体积含量参数对应的测井响应的正演结果曲线;检验所述测井响应的正演结果曲线是否与实际测井曲线匹配,将匹配的正演结果曲线对应的页岩地层最优体积含量参数中的干酪根体积含量作为确定结果。
Description
技术领域
本发明涉及页岩油气勘探开发技术领域,尤其涉及一种确定页岩地层干酪根含量的方法。
背景技术
页岩油气是存储在泥岩、页岩或者粉砂质较重的细粒沉积岩中的石油或天然气。由于页岩在矿物组成、物性和渗流特性上与砂岩有很大区别,因此,传统测井解释理论和方法不适用于页岩油气层。
干酪根是页岩地层中的重要组成部分,干酪根含量直接影响吸附气含量,干酪根含量越多,页岩地层中的吸附气越高。因此,准确识别干酪根的含量可评价页岩地层的生烃潜力和吸附气量。
目前有两种页岩地层的干酪根的求解方法,一种是由有机质含量利用经验公式计算干酪根含量(Thaimar,2011年);另一种是根据核磁共振测井结合孔隙度测井资料计算(Herron,2000年)。
前一种方法页岩地层中干酪根与有机质含量关系并不是完全一致,由于受到有机质成熟度、干酪根性质等方面的影响,不同页岩地区他们之间的转换关系不一样,确定转换系数需要实验方法和经验方法确定。后一种方法利用核磁共振测井资料获得有效孔隙度大小,利用补偿声波、岩性密度、补偿中子测井曲线计算出总孔隙度大小,总孔隙度与有效孔隙度之差为干酪根孔隙度,这种方法存在很大的局限性。例如,在没有采集核磁共振测井资料的情况下,不能应用后一种方法。
因此,亟需一种利用常规测井资料确定页岩地层干酪根含量的方法。
发明内容
本发明针对现有技术的不足,提出了一种确定页岩地层干酪根含量的方法,包括以下步骤:
建立由骨架矿物、干酪根和孔隙组成的页岩地层体积模型,其中孔隙包括气体和水;
基于所述体积模型构造页岩地层在不同体积含量参数条件下的测井响应方程,确定响应方程误差和测量误差;
采集实际测井曲线,以页岩地层矿物质量含量为约束条件,根据所述测井响应方程,响应方程误差和测量误差建立测井解释的目标函数;
计算使所述目标函数取得最小值的最优体积含量参数,绘制最优体积含量参数对应的测井响应的正演结果曲线;
检验所述测井响应的正演结果曲线是否与实际测井曲线匹配,将匹配的正演结果曲线对应的页岩地层最优体积含量参数中的干酪根体积含量作为确定结果。
根据本发明的一个实施例,所述页岩地层体积模型为:
单位体积的页岩地层包括局部均匀的多种骨架矿物、干酪根和孔隙,其中孔隙包括气体和水,满足
其中Vj为第j种骨架矿物体积含量,VK为干酪根体积含量,φg为气体体积含量,φw为水的体积含量,N为骨架矿物数量。
根据本发明的一个实施例,所述构造页岩地层在不同体积含量参数条件下的测井响应方程包括:
补偿声波测井曲线响应方程
其中MACj为第j种骨架矿物对应的声波骨架值;ACK、ACg、ACw分别是干酪根、气体和水的声波骨架值;
补偿中子测井曲线响应方程
其中MCNLj为第j种骨架矿物对应的中子骨架值;CNLK、CNLg、CNLw分别是干酪根、气体和水的中子骨架值;
补偿密度测井曲线响应方程
其中MDENj为第j种矿物对应的密度骨架值;DENK、DENg、DENw分别是干酪根、气体和水的密度骨架值;
上式中Vj为第j种骨架矿物体积含量,VK为干酪根体积含量,φg为气体体积含量,φw为水的体积含量,N为骨架矿物数量。
根据本发明的一个实施例,所述响应方程误差包括:
声波测井响应方程误差其中,δΔtmf、δΔti分别为泥浆滤液和矿物的补偿声波误差,Vj为第j种骨架矿物体积含量,φ为孔隙度;
补偿中子测井响应方程误差其中,δCNLmf、δCNLi分别为泥浆滤液和矿物的补偿中子误差;
补偿密度测井响应方程误差其中,δρmf、δρi分别为泥浆滤液和矿物的补偿密度误差;
上式中Vi为第i种骨架矿物体积含量,φ为孔隙度,N为骨架矿物数量。
根据本发明的一个实施例,所述测井解释的目标函数为:
其中,ai为第i种测井曲线实际测量值,x=(V1,…Vj,VK,φg,φw)为页岩地层体积含量参数,z为当前测井深度,σi为第i种测井曲线的测量误差,τi为构造的第i种测井曲线响应方程误差,fi(x,z)为在z深度时构造的第i种测井曲线响应,m为测井曲线的种类;
gj(x)为x的第j种约束条件,τj为第j种约束误差,p为约束个数。
根据本发明的一个实施例,所述测井解释的目标函数中的约束条件包括:
其中,Vj为第j种页岩地层骨架矿物的体积含量,wj为页岩地层的矿物质量含量,ρb为密度测井值,ρMA为矿物骨架密度。
根据本发明的一个实施例,计算使所述目标函数取得最小值的最优体积含量参数包括:
利用BFGS变尺度法确定搜索方向为其中为方向矩阵,表示目标函数的梯度,在搜索方向上任意一点表示为目标函数F(x,a)转换为单变量函数φ(t);
利用DSPOW抛物线插值法求取单变量函数φ(t)的极小值点t0,根据t0确定页岩地层最优体积含量参数x=(V1,…Vj,VK,φg,φw)。
根据本发明的一个实施例,所述检验所述测井响应的正演结果曲线是否与实际测井曲线匹配包括:根据实际测井曲线、响应方程误差和测量误差绘制正演曲线的置信区间其中ai为第i种测井曲线实际测量值;σi和τi对应为测量误差和测井曲线响应方程误差;
当所述测井响应的正演结果曲线处于正演曲线的置信区间内时,所述检验所述测井响应的正演结果曲线与实际测井曲线匹配。
根据本发明的一个实施例,所述第i种测井曲线的测量误差为实际测量数值大小的5%。
根据本发明的一个实施例,所述页岩地层矿物质量含量由全岩分析法或者ECS元素俘获测井资料解释确定。
本发明带来了以下有益效果:根据矿物质量解释结果,结合常规测井资料中的补偿声波、补偿中子和密度测井资料反演计算出页岩地层矿物体积、干酪根体积和孔隙度大小,可以快速识别页岩地层干酪根含量,减少了岩心实验确定干酪根含量的过程,节省了费用和时间。
本发明的其它特征和优点将在随后的说明书中阐述,并且部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
图1是本发明实施例的技术思路示意图;
图2是本发明实施例一的确定页岩地层干酪根含量的方法的步骤流程图;
图3是本发明实施例一的页岩地层体积模型示意图;
图4是本发明实施例二的实际工程中某井页岩地层干酪根含量反演结果。
具体实施方式
以下将结合附图来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明各实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
另外,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机***中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
如图1所示,本发明实施例的技术思路为,构造响应方程(补偿声波、补偿中子、补偿密度)的理论测井值,根据矿物质量含量结果和常规测井资料(补偿声波、补偿中子、补偿密度)综合优化反演计算页岩地层矿物体积、干酪根体积和孔隙度大小。
实施例一
本实施例提供一种确定页岩地层干酪根含量的方法,以下结合图2对本实施例提供的方法作详细说明。
在步骤S201中,建立由骨架矿物、干酪根和孔隙组成的页岩地层体积模型。参考图3,单位体积的页岩地层包括局部均匀的多种骨架矿物、干酪根和孔隙,其中孔隙包括气体和水,满足
其中Vj为第j种骨架矿物体积含量,VK为干酪根体积含量,φg为气体体积含量,φw为水的体积含量,N为骨架矿物数量。V1,…Vj,VK,φg,φw表示某一地层深度处的体积含量。
具体的,不同地区矿物成份并不相同,常见的骨架矿物包括:伊利石、绿泥石、蒙脱石、石英、长石、云母方解石、白云岩、硬石膏和黄铁矿等。通常情况下,页岩地层中的硅质包括石英、长石和云母,钙质包括方解石和白云岩,泥质包括蒙脱石、伊利石和绿泥石。
在步骤S202中,基于所述体积模型构造页岩地层在不同体积含量参数条件下的测井响应方程,确定响应方程误差和测量误差。
页岩地层不同骨架矿物、干酪根和孔隙条件下的测井响应方程为:
其中fi为第i种测井曲线的响应结果;Vj为第j种页岩地层矿物体积含量;VK为干酪根体积含量;φg为气体体积含量,φw为水的体积含量;MAij为第i种测井曲线对应第j种矿物的骨架值;MAiK、MAig、MAiw分别为干酪根、气体和水对应的第i种曲线骨架值,N为骨架矿物数量。
在本实施例中,构造不同体积含量参数条件下的补偿声波测井、补偿中子测井和补偿密度测井响应方程:
补偿声波测井曲线响应方程
其中MACj为第j种骨架矿物对应的声波骨架值;ACK、ACg、ACw分别是干酪根、气体和水的声波骨架值;
补偿中子测井曲线响应方程
其中MCNLj为第j种骨架矿物对应的中子骨架值;CNLK、CNLg、CNLw分别是干酪根、气体和水的中子骨架值;
补偿密度测井曲线响应方程
其中MDENj为第j种矿物对应的密度骨架值;DENK、DENg、DENw分别是干酪根、气体和水的密度骨架值。
由于测井仪器与数据采集***均存在一定的测量误差,本实施例中某种测井响应的测量误差σi取其实际测量数值大小的5%。
任何测井方程都是根据对实际地层作一系列的数学物理简化之后建立的测井解释模型得出的理论公式,响应方程中存在一定的响应方程误差τi。
在步骤S202中,响应方程误差包括:
声波测井响应方程误差其中,δΔtmf、δΔti分别为泥浆滤液和骨架矿物的补偿声波误差,Vj为第j种骨架矿物体积含量,φ为孔隙度,φ=φg+φw,如上文所述,φg为气体体积含量,φw为水的体积含量;
补偿中子测井响应方程误差
其中,δCNLmf、δCNLi分别为泥浆滤液和骨架矿物的补偿中子误差;
补偿密度测井响应方程误差
其中,δρmf、δρi分别为泥浆滤液和骨架矿物的密度误差。
在步骤S203中,采集实际测井曲线,以页岩地层矿物质量含量为约束条件,根据所述测井响应方程,响应方程误差和测量误差建立测井解释的目标函数:
其中,
ai为第i种测井曲线实际测量值,%;
x=(V1,…Vj,VK,φg,φw)为页岩地层体积含量参数,%;
z为当前测井深度,m;
σi为第i种测井曲线的测量误差;
τi为构造的第i种测井曲线响应方程误差;
fi(x,z)为在z深度时构造的第i种测井曲线响应;
m为测井曲线的种类;
gj(x)为x的第j种约束条件;
τj为第j种约束误差;
p为约束个数。
在本实施例中,fi(x,z)分别为补偿声波测井曲线响应AC,补偿中子测井曲线响应CNL,补偿密度测井曲线响应DEN。
优选的,本实施例中将页岩地层的矿物质量含量作为约束条件:
约束条件gj(x)表示为其中,
Vj为第j种页岩地层骨架矿物的体积含量,wj为页岩地层的矿物质量含量,ρb为密度测井值,ρMA为矿物骨架密度。页岩地层矿物质量含量wj由全岩分析法或者ECS元素俘获测井资料确定。
相应的质量约束误差τj取值为1。
不限于此,还包括以下约束条件:
Vj≥0,j为矿物种类;φg≥0,φw≥0。
在步骤S204中,计算使所述目标函数取得最小值的最优体积含量参数,绘制最优体积含量参数对应的测井响应的正演结果曲线。
由目标参数F(x,a)可构建最优化测井解释的数学模型:
计算使所述目标函数F(x,a)取得最小值的最优体积含量参数x=(V1,…Vj,VK,φg,φw),即可确定与实际体积含量最接近的体积含量参数。
优选的,在本实施例中利用BFGS变尺度法确定搜索方向为在搜索方向上可寻找目标函数的极小值。
其中正定对称矩阵为尺度矩阵,表示目标函数F(x,a)在第k次迭代点处的梯度,在第k次搜索方向上任意一点表示为目标函数F(x,a)转换为单变量函数φ(t)。
利用DSPOW抛物线插值法求取单变量函数φ(t)的极小值点t0,根据t0确定页岩地层最优体积含量参数x=(V1,…Vj,VK,φg,φw)。
上文中BFGS变尺度法和DSPOW抛物线插值法是本领域技术人员公知的测井解释最优化方法,可参考《高等学校教材:测井数据处理与综合解释》(雍世和著,2007年8月1日出版)。
在步骤S205中,检验所述测井响应的正演结果曲线是否与实际测井曲线匹配,将匹配的正演结果曲线对应的体积含量参数中的干酪根体积含量作为确定结果。
具体的,根据实际测井曲线、响应方程误差和测量误差绘制正演曲线的的置信区间其中ai为第i种测井曲线实际测量值;σi和τi对应为测量误差和测井曲线响应方程误差;
当所述测井响应的正演结果曲线处于正演曲线的置信区间内时,所述检验所述测井响应的正演结果曲线与实际测井曲线匹配,即说明反演得到的最优体积含量参数x=(V1,…Vj,VK,φg,φw)充分反应了页岩地层中的实际体积含量。
此外,可根据目标函数最优值的大小检验测井响应的正演结果曲线是否与实际测井曲线匹配。对于利用最优化方法计算的最优体积含量参数x*,目标函数最优值F*(x*,a)≈0。F*(x*,a)是衡量正演测井曲线是否逼近实际测井曲线的总标志。F*(x*,a)值较小,说明测井响应的正演结果曲线与实际测井曲线匹配程度较好。
综上所述,本发明实施例提供的方法可根据页岩地层矿物质量解释结果,结合常规测井资料中的补偿声波、补偿中子和密度测井资料反演计算出页岩地层矿物体积、干酪根体积和孔隙度大小,可以快速识别页岩地层干酪根含量,减少了岩心实验确定干酪根含量的过程,节省了费用和时间。
实施例二
图4是本发明实施例二的实际工程中某井页岩地层干酪根含量反演结果。利用中国西南部区块某井ECS元素测井资料进行反演。该区块优质页岩岩性主要以黄灰色页岩、粉砂质页岩夹薄层透镜状灰岩为主,页岩地层矿物干酪根含量确定是利用元素测井资料计算矿物质量结果作为输入曲线,结合常规测井中的补偿声波、补偿密度和补偿中子资料进行反演,反演结果如图4所示。
图4中第一道为利用ECS元素测井资料计算的矿物质量含量组合,第二道至第四道为补偿声波、补偿密度、补偿密度的正演响应值和正演曲线的置信区间,第六道为目标函数最优值,第七道为反演的矿物体积含量、干酪根和孔隙度大小,其中孔隙度大小包括气体和水的体积含量。
图4中第一道为根据元素测井资料计算的矿物质量含量组合,分别表示出以下质量含量:黄铁矿WPYR、伊利石WILL、绿泥石WCHL、蒙脱石WMON、石英WQRZ、长石WFLD、方解石WCLC、硬石膏WANH、云母WMIC和白云岩WDOL。
第二道至第四道为补偿声波、补偿密度、岩性密度的正演响应值和正演曲线的置信区间,第五道为深度道,第六道为目标函数最优值。
第二道中,曲线ACL和ACH表示补偿声波正演曲线的置信区间;曲线AC表示实际的补偿声波测井曲线;曲线ACT表示正演的补偿声波测井响应最优值。类似的,曲线CNLL和CNLH表示补偿中子正演曲线的置信区间;曲线CNL表示实际的补偿中子测井曲线;曲线CNLT表示正演的补偿中子测井响应最优值。类似的,曲线DENL和DENH表示补偿密度正演曲线的置信区间;曲线DEN表示实际的补偿密度测井曲线;曲线DENT表示正演的补偿密度测井响应最优值。
第七道为反演的矿物体积含量、干酪根、气和水的体积含量。分别表示出以下体积含量:黄铁矿VPYR、伊利石VILL、绿泥石VCHL、蒙脱石VMON、石英VQRZ、长石VFLD、方解石VCLC、硬石膏VANH、云母VMIC和白云岩VDOL,以及干酪根VKER、气体VPRG、水VPGW和其他的体积含量VSM。从而可以准确得到在不同地层深度处干酪根的体积含量VKER。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种确定页岩地层干酪根含量的方法,其特征在于,包括以下步骤:
建立由骨架矿物、干酪根和孔隙组成的页岩地层体积模型,其中孔隙包括气体和水;
基于所述体积模型构造页岩地层在不同体积含量参数条件下的测井响应方程,确定响应方程误差和测量误差;
采集实际测井曲线,以页岩地层矿物质量含量为约束条件,根据所述测井响应方程,响应方程误差和测量误差建立测井解释的目标函数;
计算使所述目标函数取得最小值的最优体积含量参数,绘制最优体积含量参数对应的测井响应的正演结果曲线;
检验所述测井响应的正演结果曲线是否与实际测井曲线匹配,将匹配的正演结果曲线对应的页岩地层最优体积含量参数中的干酪根体积含量作为确定结果。
2.根据权利要求1所述的方法,其特征在于,所述页岩地层体积模型为:
单位体积的页岩地层包括局部均匀的多种骨架矿物、干酪根和孔隙,其中孔隙包括气体和水,满足
其中Vj为第j种骨架矿物体积含量,VK为干酪根体积含量,φg为气体体积含量,φw为水的体积含量,N为骨架矿物数量。
3.根据权利要求2所述的方法,其特征在于,所述构造页岩地层在不同体积含量参数条件下的测井响应方程包括:
补偿声波测井曲线响应方程
其中MACj为第j种骨架矿物对应的声波骨架值;ACK、ACg、ACw分别是干酪根、气体和水的声波骨架值;
补偿中子测井曲线响应方程
其中MCNLj为第j种骨架矿物对应的中子骨架值;CNLK、CNLg、CNLw分别是干酪根、气体和水的中子骨架值;
补偿密度测井曲线响应方程
其中MDENj为第j种矿物对应的密度骨架值;DENK、DENg、DENw分别是干酪根、气体和水的密度骨架值;
上式中Vj为第j种骨架矿物体积含量,VK为干酪根体积含量,φg为气体体积含量,φw为水的体积含量,N为骨架矿物数量。
4.根据权利要求3所述的方法,其特征在于,所述响应方程误差包括:
声波测井响应方程误差其中,δΔtmf、δΔti分别为泥浆滤液和矿物的补偿声波误差;
补偿中子测井响应方程误差 其中,δCNLmf、δCNLi分别为泥浆滤液和矿物的补偿中子误差;
补偿密度测井响应方程误差其中,δρmf、δρi分别为泥浆滤液和矿物的补偿密度误差;
上式中,Vi为第i种骨架矿物体积含量,φ为孔隙度,N为骨架矿物数量。
5.根据权利要求4所述的方法,其特征在于,所述测井解释的目标函数为:
其中,ai为第i种测井曲线实际测量值,x=(V1,…Vj,VK,φg,φw)为页岩地层体积含量参数,z为当前测井深度,σi为第i种测井曲线的测量误差,τi为构造的第i种测井曲线响应方程误差,fi(x,z)为在z深度时构造的第i种测井曲线响应,m为测井曲线的种类;
gj(x)为x的第j种约束条件,τj为第j种约束误差,p为约束个数。
6.根据权利要求5所述的方法,其特征在于,所述测井解释的目标函数中的约束条件包括:
其中,Vj为第j种页岩地层骨架矿物的体积含量,wj为页岩地层的矿物质量含量,ρb为密度测井值,ρMA为矿物骨架密度。
7.根据权利要求6所述的方法,其特征在于,计算使所述目标函数取得最小值的最优体积含量参数包括:
利用BFGS变尺度法确定搜索方向为其中为方向矩阵,表示目标函数的梯度,在搜索方向上任意一点表示为目标函数F(x,a)转换为单变量函数φ(t);
利用DSPOW抛物线插值法求取单变量函数φ(t)的极小值点t0,根据t0确定页岩地层最优体积含量参数x=(V1,…Vj,VK,φg,φw)。
8.根据权利要求6所述的方法,其特征在于,所述检验所述测井响应的正演结果曲线是否与实际测井曲线匹配包括:
根据实际测井曲线、响应方程误差和测量误差绘制正演曲线的置信区间其中ai为第i种测井曲线实际测量值;σi和τi对应为测量误差和测井曲线响应方程误差;
当所述测井响应的正演结果曲线处于正演曲线的置信区间内时,所述检验所述测井响应的正演结果曲线与实际测井曲线匹配。
9.根据权利要求7或8所述的方法,其特征在于,所述第i种测井曲线的测量误差为实际测量数值大小的5%。
10.根据权利要求7或8所述的方法,其特征在于,所述页岩地层矿物质量含量由全岩分析法或者ECS元素俘获测井资料解释确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410319217.9A CN105317429A (zh) | 2014-07-04 | 2014-07-04 | 一种确定页岩地层干酪根含量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410319217.9A CN105317429A (zh) | 2014-07-04 | 2014-07-04 | 一种确定页岩地层干酪根含量的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105317429A true CN105317429A (zh) | 2016-02-10 |
Family
ID=55245632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410319217.9A Pending CN105317429A (zh) | 2014-07-04 | 2014-07-04 | 一种确定页岩地层干酪根含量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105317429A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106596605A (zh) * | 2016-11-29 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种制备高纯干酪根的方法 |
CN106950347A (zh) * | 2017-03-03 | 2017-07-14 | 中国石油大学(华东) | 一种评价泥页岩各组分体积的方法 |
CN107451671A (zh) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | 用于预测页岩地层压裂后初始产能的方法及*** |
CN108204233A (zh) * | 2017-12-20 | 2018-06-26 | 中国石油大学(北京) | 获取页岩气的类型的方法和装置 |
CN109188557A (zh) * | 2018-09-29 | 2019-01-11 | 中国石油天然气股份有限公司 | 一种应用生烃增压模型计算排烃效率的方法 |
CN112368607A (zh) * | 2018-04-23 | 2021-02-12 | 斯伦贝谢技术有限公司 | 用于改进地下地层中的中子解释的方法 |
CN113123782A (zh) * | 2019-12-31 | 2021-07-16 | 中国石油化工股份有限公司 | 一种利用线性规划评价地层组分的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6484102B1 (en) * | 2000-08-24 | 2002-11-19 | Digital Formation Inc. | System for evaluating fluid distributions of subsurface reservoirs |
US7054749B1 (en) * | 2000-11-13 | 2006-05-30 | O'meara Jr Daniel J | Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models |
GB2421794B (en) * | 2004-12-20 | 2007-05-02 | Schlumberger Holdings | Method of evaluating fluid saturation characteristics in a geological formation |
CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
CN103603659A (zh) * | 2013-12-09 | 2014-02-26 | 中国石油集团川庆钻探工程有限公司 | 利用常规测井资料识别复杂孔隙结构储层流体类型的方法 |
-
2014
- 2014-07-04 CN CN201410319217.9A patent/CN105317429A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6484102B1 (en) * | 2000-08-24 | 2002-11-19 | Digital Formation Inc. | System for evaluating fluid distributions of subsurface reservoirs |
US7054749B1 (en) * | 2000-11-13 | 2006-05-30 | O'meara Jr Daniel J | Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models |
GB2421794B (en) * | 2004-12-20 | 2007-05-02 | Schlumberger Holdings | Method of evaluating fluid saturation characteristics in a geological formation |
CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
CN103603659A (zh) * | 2013-12-09 | 2014-02-26 | 中国石油集团川庆钻探工程有限公司 | 利用常规测井资料识别复杂孔隙结构储层流体类型的方法 |
Non-Patent Citations (4)
Title |
---|
孙建孟等: "《地球物理资料综合应用》", 30 September 2001, 石油大学出版社 * |
朱伟林等: "《深海油气的奥秘》", 30 November 2013, 石油工业出版社 * |
潘仁芳等: "页岩气测井技术的应用", 《中国科技信息》 * |
潘保芝等: "《裂缝性火成岩储层测井评价方法与应用》", 31 October 2003, 石油工业出版社 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451671A (zh) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | 用于预测页岩地层压裂后初始产能的方法及*** |
CN106596605A (zh) * | 2016-11-29 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种制备高纯干酪根的方法 |
CN106596605B (zh) * | 2016-11-29 | 2019-11-08 | 中国石油天然气股份有限公司 | 一种制备高纯干酪根的方法 |
CN106950347A (zh) * | 2017-03-03 | 2017-07-14 | 中国石油大学(华东) | 一种评价泥页岩各组分体积的方法 |
CN106950347B (zh) * | 2017-03-03 | 2019-03-05 | 中国石油大学(华东) | 一种评价泥页岩各组分体积的方法 |
CN108204233A (zh) * | 2017-12-20 | 2018-06-26 | 中国石油大学(北京) | 获取页岩气的类型的方法和装置 |
CN112368607A (zh) * | 2018-04-23 | 2021-02-12 | 斯伦贝谢技术有限公司 | 用于改进地下地层中的中子解释的方法 |
CN109188557A (zh) * | 2018-09-29 | 2019-01-11 | 中国石油天然气股份有限公司 | 一种应用生烃增压模型计算排烃效率的方法 |
CN109188557B (zh) * | 2018-09-29 | 2020-07-10 | 中国石油天然气股份有限公司 | 一种应用生烃增压模型计算排烃效率的方法 |
CN113123782A (zh) * | 2019-12-31 | 2021-07-16 | 中国石油化工股份有限公司 | 一种利用线性规划评价地层组分的方法 |
CN113123782B (zh) * | 2019-12-31 | 2024-05-14 | 中国石油化工股份有限公司 | 一种利用线性规划评价地层组分的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105317429A (zh) | 一种确定页岩地层干酪根含量的方法 | |
Zhao et al. | An improved model for estimating the TOC in shale formations | |
CN103617345B (zh) | 用元素含量计算地层岩石矿物成分的方法 | |
Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
CN103954999B (zh) | 一种适用于低孔隙度砂泥岩地层的横波速度预测方法 | |
CN103576195B (zh) | 一种随压力变化的裂隙介质横波速度预测方法 | |
CN104407378B (zh) | 一种各向异性参数反演方法及装置 | |
CN104573150A (zh) | 一种针对富含有机质页岩的各向异性岩石物理建模方法 | |
AU2012385250B2 (en) | Anisotropy parameter estimation | |
CN105370272A (zh) | 页岩气储层的测井评价方法 | |
CN105527653B (zh) | 一种基于地质信息的虚拟井构建方法 | |
CN105089615B (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN104632202B (zh) | 确定干粘土三孔隙度测井参数值的方法及装置 | |
CN103867198B (zh) | 碳酸盐岩天然气层地层密度判别法 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN113901681A (zh) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 | |
CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
CN103485758A (zh) | 一种利用多种资料拟合反演复杂岩性解释方法及装置 | |
CN105221142A (zh) | 一种识别页岩地层矿物质量含量的方法 | |
CN104865614A (zh) | 基于可变骨架参数的复杂储层流体识别方法 | |
CN105298476B (zh) | 一种识别页岩地层矿物体积含量的方法 | |
CN111368245B (zh) | 一种考虑多影响因素的定量计算岩石Biot系数方法 | |
CN106842301A (zh) | 一种凝灰质砂岩有利储层的定量识别与预测方法 | |
CN110501744A (zh) | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 | |
CN105275459A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160210 |
|
RJ01 | Rejection of invention patent application after publication |