发明名称 一种泥页岩油气饱和度的计算模型
摘要 本发明公开了一种泥页岩油气饱和度的计算模型,包括以下步骤:A.建立等效物理模型;B.给出假设条件;C.确定数学模型公式;D.模型中各参数的确定;E.岩电实验数据拟合模型公式参数a、b、m、n;F.模型验证。本发明能够解决泥页岩储层由于骨架矿物复杂、孔隙结构复杂、孔隙润湿性复杂等问题带来的非阿尔奇现象,准确计算泥页岩储层含油气饱和度,为测井评价储层以及油田开发提供帮助,在计算泥页岩储层含水饱和度方面具备其它饱和度模型无可比拟的优势,实际应用效果显著,因此极具推广价值。在目前公开发表文献和商业应用软件中尚无类似方法的提出与应用。
申请公布号 CN104712329A 申请公布日期 2015.06.17
申请号 CN201510044729.3 申请日期 2015.01.29
申请人 中国石油大学(华东) 发明人 孙建孟;张晋言;赵建鹏;侯庆功;李绍霞;张文姣;闫伟超
分类号 E21B49/00(2006.01)I;G06F19/00(2011.01)I 主分类号 E21B49/00(2006.01)I
代理机构 北京方圆嘉禾知识产权代理有限公司 11385 代理人 董芙蓉
主权项 一种泥页岩油气饱和度的计算模型,首先根据泥页岩骨架矿物组分与孔隙类型给出泥页岩岩石体积物理模型,然后在适当的假设条件下,根据物理模型推导出泥页岩导电数学公式,随后给出公式中各项参数的求取方法,具体步骤如下:A.建立等效物理模型:根据泥页岩骨架矿物类型和孔隙类型给出泥页岩岩石体积物理模型,将页岩分为两大系统和四大孔隙组分;B.给出假设条件:模型的提出是基于以下假设的:(1)在泥页岩储层,骨架矿物除黄铁矿以外,有机质和其余骨架矿物是完全不导电的,在泥页岩储层中主要含有四种孔隙类型,分别为粘土孔隙、基质孔隙、微裂缝和有机孔隙;(2)粘土孔隙、基质孔隙、微裂缝是亲水的水润湿、有机质孔隙是亲油气的油润湿,并且有机质孔隙中是100%含油气的,即油气饱和度为100%;粘土孔隙是粘土矿物之间的微细孔隙,与有机质孔隙不同,粘土孔隙表面表现为亲水性特征,对水分子吸附能力强;在构建游离气体积计算模型时,假设粘土孔隙100%含水,不含油气;(3)泥页岩储层基质孔隙主要为残余原生孔隙和不稳定矿物溶蚀孔;(4)岩石的电阻率主要受四部分影响:①100%含水的粘土;②导电矿物黄铁矿;③有机质;④基质孔隙和微裂隙空间的地层水;假设以上导电成分是并联导电的;C.确定数学模型公式:根据泥页岩岩石体积物理模型,认为泥页岩电阻率主要受粘土、黄铁矿、有机质四部分影响:计算每一部分的电阻率然后并联起来,最终得到泥页岩饱和度导电模型;构建如下页岩并联导电模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mn>1</mn><msub><mi>R</mi><mi>t</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>V</mi><mi>cl</mi></msub><msub><mi>R</mi><mi>cl</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>py</mi><mo>_</mo><mi>cut</mi></mrow></msub></mrow><msub><mi>R</mi><mi>py</mi></msub></mfrac><mo>+</mo><mfrac><msub><mi>V</mi><mi>TOC</mi></msub><msub><mi>R</mi><mi>TOC</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msup><mi>&phi;</mi><mi>m</mi></msup><msubsup><mi>S</mi><mi>w</mi><mi>n</mi></msubsup></mrow><mrow><mi>ab</mi><msub><mi>R</mi><mi>w</mi></msub></mrow></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>cl</mi></msub><mo>-</mo><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000663688890000021.GIF" wi="1131" he="156" /></maths>在实际处理过程中,通常认为有机质不导电,所以上式可以化简为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mn>1</mn><msub><mi>R</mi><mi>t</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>V</mi><mi>cl</mi></msub><msub><mi>R</mi><mi>cl</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>py</mi><mo>_</mo><mi>cut</mi></mrow></msub></mrow><msub><mi>R</mi><mi>py</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msup><mi>&phi;</mi><mi>m</mi></msup><msubsup><mi>S</mi><mi>w</mi><mi>n</mi></msubsup></mrow><mrow><mi>ab</mi><msub><mi>R</mi><mi>w</mi></msub></mrow></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>cl</mi></msub><mo>-</mo><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mi>TOC</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000663688890000022.GIF" wi="1028" he="156" /></maths>变形可得:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mi>w</mi></msub><mo>=</mo><mi>n</mi><msqrt><mrow><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>R</mi><mi>t</mi></msub></mfrac><mo>-</mo><mfrac><msub><mi>V</mi><mi>cl</mi></msub><msub><mi>R</mi><mi>cl</mi></msub></mfrac><mo>-</mo><mfrac><mrow><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>py</mi><mo>_</mo><mi>cut</mi></mrow></msub></mrow><msub><mi>R</mi><mi>py</mi></msub></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><mi>ab</mi><msub><mi>R</mi><mi>w</mi></msub></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>cl</mi></msub><mo>-</mo><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mi>TOC</mi></msub><mo>)</mo></mrow><msup><mi>&phi;</mi><mi>m</mi></msup></mrow></mfrac></mrow></msqrt></mrow>]]></math><img file="FDA0000663688890000023.GIF" wi="1219" he="206" /></maths>则油气饱和度为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>S</mi><mi>h</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mi>n</mi><msqrt><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>R</mi><mi>t</mi></msub></mfrac><mo>-</mo><mfrac><msub><mi>V</mi><mi>cl</mi></msub><msub><mi>R</mi><mi>cl</mi></msub></mfrac><mo>-</mo><mfrac><mrow><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>py</mi><mo>_</mo><mi>cut</mi></mrow></msub></mrow><msub><mi>R</mi><mi>py</mi></msub></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><mi>ab</mi><msub><mi>R</mi><mi>w</mi></msub></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>cl</mi></msub><mo>-</mo><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mi>TOC</mi></msub><mo>)</mo></mrow><msup><mi>&phi;</mi><mi>m</mi></msup></mrow></mfrac></msqrt></mrow>]]></math><img file="FDA0000663688890000024.GIF" wi="1281" he="208" /></maths>式中,R<sub>t</sub>表示页岩的电阻率,实际计算时可以用深探测电阻率代替;R<sub>cl</sub>为100%含水粘土的电阻率,V<sub>cl</sub>为粘土的体积百分数;R<sub>py</sub>为黄铁矿的电阻率,V<sub>py</sub>为黄铁矿的体积百分数;V<sub>py_cut</sub>是黄铁矿体积百分数截止值,只有当当黄铁矿体积百分含量大于截止值时才能构成连通相,具有导电性;V<sub>TOC</sub>为有机质体积百分数,R<sub>TOC</sub>为有机质电阻率;φ为基质孔隙、微裂隙之和,S<sub>w</sub>为以上两类孔隙的含水饱和度,R<sub>w</sub>为地层水电阻率,a、b、m、n为模型参数;D.模型中各参数的确定:①粘土含量、黄铁矿含量、有机质含量的计算:粘土含量V<sub>cl</sub>可以由自然伽马测井或自然伽马能谱测井直接计算也可以由测井曲线与岩心分析数据统计回归计算;黄铁矿含量V<sub>py</sub>可以由测井曲线与岩心分析统计回归直接计算;黄铁矿含量截止值V<sub>py_cut</sub>可以通过数字岩心模拟技术给出;随着黄铁矿含量的增大,岩心的电阻率有降低的趋势;但是在黄铁矿含量低的时候,电阻率几乎没有变化,这是由于黄铁矿分散在岩心中,只有在达到一定含量的时候,与孔隙相相连成为连续相,才对岩心电阻率产生影响,V<sub>py_cut</sub>约为5%到6%之间;有机质V<sub>TOC</sub>由于密度低、声波传播速度第(低)、含氢指数高以及不导电的特性,造成鲜明的测井响应特征,同时富有机质页岩还表现为高自然伽马特征,在自然伽马能谱测井响应上表现为高铀含量特征,可以用这些测井特征定量计算有机质含量,也可以由测井曲线与岩心分析统计回归直接计算;②粘土孔隙度、有机质孔隙度、无机质孔隙度的计算:粘土孔隙是粘土矿物之间的微细孔隙,是束缚水的主要赋存空间;粘土孔隙度的计算公式如下:φ<sub>cl</sub>=φ<sub>tcl</sub>V<sub>cl</sub>式中φ<sub>cl</sub>为粘土孔隙度,φ<sub>tcl</sub>为纯粘土段孔隙度,V<sub>cl</sub>由邻近泥岩测井响应确定;利用扫描电镜测试技术可以直观的确定有机孔隙大小及其分布,估算有机孔隙的面孔率;利用平均面孔率对测井计算的V<sub>TOC</sub>进行刻度,获得有机质孔隙度,其计算公式如下:φ<sub>TOC</sub>=aV<sub>TOC</sub>其中a为有机质面孔率,由扫描电镜分析确定;页岩储层总孔隙度可以通过岩心刻度测井的方法求的;通过实验测量的氦气孔隙度与孔隙度敏感测井曲线,有声波、密度或中子曲线建立关系,可以计算储层总孔隙度φ<sub>t</sub>;因此基质孔隙、微裂隙之和φ可以通过下式计算φ=φ<sub>t</sub>‑φ<sub>TOC</sub>‑φ<sub>cl</sub>;E.岩电实验数据拟合模型公式参数a、b、m、n:利用岩电实验数据回归得出:当没有岩心实验数据时,可以采用默认值,一般认为a=b=1,对于微裂缝发育型,m值约为1.1‑1.5,n值约为1;对于孔隙发育型,m值约为1.5‑2.5,n值约为2;F.模型验证:取a=b=1,m=1.5,n=1.3,φ=0.05,R<sub>cl</sub>=40Ω,R<sub>PY</sub>=0.5Ω,R<sub>W</sub>=0.12,V<sub>cl</sub>=0.2,V<sub>PY</sub>=0.05,TOC=0.05,ρ<sub>TOC</sub>=1.2g/cm<sup>3</sup>,ρ<sub>b</sub>=2.62g/cm<sup>3</sup>;利用如下公式计算了泥页岩导电性与含水饱和度之间的关系,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>S</mi><mi>w</mi></msub><mo>=</mo><mi>n</mi><msqrt><mrow><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>R</mi><mi>t</mi></msub></mfrac><mo>-</mo><mfrac><msub><mi>V</mi><mi>cl</mi></msub><msub><mi>R</mi><mi>cl</mi></msub></mfrac><mo>-</mo><mfrac><mrow><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>py</mi><mo>_</mo><mi>cut</mi></mrow></msub></mrow><msub><mi>R</mi><mi>py</mi></msub></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><mi>ab</mi><msub><mi>R</mi><mi>w</mi></msub></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>cl</mi></msub><mo>-</mo><msub><mi>V</mi><mi>py</mi></msub><mo>-</mo><msub><mi>V</mi><mi>TOC</mi></msub><mo>)</mo></mrow><msup><mi>&phi;</mi><mi>m</mi></msup></mrow></mfrac></mrow></msqrt></mrow>]]></math><img file="FDA0000663688890000041.GIF" wi="1244" he="201" /></maths>并与岩石物理实验结果进行对比,对比结果具有很好的吻合性。
地址 266580 山东省青岛市经济技术开发区长江西路66号