发明名称 |
一种基于等效流体声波阻抗的流体识别方法 |
摘要 |
本发明公开了勘探技术领域内的一种基于等效流体声波阻抗的流体识别方法,其先通过测井方法或岩石物理测试方法获得储层的数据集合,再构建饱和流体的岩石密度,然后利用反演得到饱和流体岩石的声波阻抗,将声波阻抗分解成骨架等效声波阻抗与等效流体声波阻抗之和,利用反演的声波阻抗及骨架等效声波阻抗计算等效流体声波阻抗;等效流体声波阻抗的数据集合用石油物探专业的地震显示软件进行图形显示,并依据显示图形进行流体识别。该方法能够反映岩石的孔隙度和所含流体性质,对地震资料要求简单,且计算量小、计算相对简单。能够更好地进行流体识别。 |
申请公布号 |
CN104316959B |
申请公布日期 |
2017.02.15 |
申请号 |
CN201410593709.7 |
申请日期 |
2014.10.29 |
申请人 |
中国石油化工股份有限公司;中国石油化工股份有限公司江苏油田分公司 |
发明人 |
梁兵;鲍祥生;丁建荣 |
分类号 |
G01V1/28(2006.01)I |
主分类号 |
G01V1/28(2006.01)I |
代理机构 |
南京苏科专利代理有限责任公司 32102 |
代理人 |
董旭东 |
主权项 |
一种基于等效流体声波阻抗的流体识别方法,其特征在于包括如下步骤:1)使用测井解释方法或岩石物理测试方法获取储层不同位置油的密度、水的密度、地震波传播速度和孔隙度,形成一个数据集合;2)利用基质密度、孔隙中的流体密度和孔隙度构建饱和流体的岩石密度:ρ=ρ<sub>ma</sub>(1‑φ)+ρ<sub>f</sub>φ (2)其中,ρ为饱和流体的岩石密度,ρ<sub>ma</sub>为基质密度,ρ<sub>f</sub>为孔隙流体的密度,φ为孔隙度;所述基质为岩石中的骨架部分;同时利用声波在骨架中的传播速度及在流体中的传播速度构建地震波在岩石中的传播速度的关系式:v=(1‑φ)<sup>2</sup>v<sub>ma</sub>+φv<sub>f</sub> (4)其中,v<sub>ma</sub>为声波在骨架中的传播速度,v<sub>f</sub>为声波在流体中的传播速度,3)利用反演方法得到饱和流体岩石的声波阻抗AI,AI可以表示如下公式:AI=ρν (1)将(2)和(4)带入(1)得到:<maths num="0001"><math><![CDATA[<mrow><mi>A</mi><mi>I</mi><mo>=</mo><msub><mi>ρ</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mo>[</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>φ</mi><mrow><mo>(</mo><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>ρ</mi><mi>f</mi></msub><msub><mi>v</mi><mi>f</mi></msub><mfrac><mi>φ</mi><mrow><mo>(</mo><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>)</mo></mrow></mfrac><mo>[</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>φ</mi><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000011.GIF" wi="1668" he="223" /></maths>其中,AI为饱和流体岩石的声波阻抗;4)将步骤3)的声波阻抗AI分解成骨架等效声波阻抗<img file="FDA0001124403920000012.GIF" wi="94" he="63" />与等效流体声波阻抗<img file="FDA0001124403920000013.GIF" wi="80" he="71" />之和,可得到:<maths num="0002"><math><![CDATA[<mrow><mi>A</mi><mi>I</mi><mo>=</mo><msubsup><mi>AI</mi><mrow><mi>m</mi><mi>a</mi></mrow><mi>φ</mi></msubsup><mo>+</mo><msubsup><mi>AI</mi><mi>f</mi><mi>φ</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000014.GIF" wi="1726" he="74" /></maths>其中,<img file="FDA0001124403920000021.GIF" wi="91" he="63" />是岩石骨架贡献部分,称之为骨架等效声波阻抗;<img file="FDA0001124403920000022.GIF" wi="85" he="69" />是岩石中流体贡献部分,称之为等效流体声波阻抗,限定:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>AI</mi><mrow><mi>m</mi><mi>a</mi></mrow><mi>φ</mi></msubsup><mo>=</mo><msub><mi>ρ</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mo>[</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>φ</mi><mrow><mo>(</mo><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000023.GIF" wi="1363" he="161" /></maths><maths num="0004"><math><![CDATA[<mrow><msubsup><mi>AI</mi><mi>f</mi><mi>φ</mi></msubsup><mo>=</mo><msub><mi>ρ</mi><mi>f</mi></msub><msub><mi>v</mi><mi>f</mi></msub><mfrac><mi>φ</mi><mrow><mo>(</mo><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>)</mo></mrow></mfrac><mo>[</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>φ</mi><mfrac><msub><mi>v</mi><mi>f</mi></msub><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi></mrow></msub></mfrac><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000024.GIF" wi="1382" he="222" /></maths>令<img file="FDA0001124403920000025.GIF" wi="205" he="125" />R<sub>v</sub>是声波在流体中的传播速度与声波在骨架中的传播速度的比值;R<sub>v</sub>取值为[2.5~3.0],公式(8)和公式(9)又可以表示成:<img file="FDA0001124403920000026.GIF" wi="1406" he="87" /><maths num="0005"><math><![CDATA[<mrow><msubsup><mi>AI</mi><mi>f</mi><mi>φ</mi></msubsup><mo>=</mo><msub><mi>ρ</mi><mi>f</mi></msub><msub><mi>v</mi><mi>f</mi></msub><mfrac><mi>φ</mi><msub><mi>R</mi><mi>v</mi></msub></mfrac><mo>[</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>φ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>φR</mi><mi>v</mi></msub><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000027.GIF" wi="1310" he="127" /></maths>5)基于步骤4)的分解公式,利用反演的声波阻抗及骨架等效声波阻抗计算等效流体声波阻抗;<maths num="0006"><math><![CDATA[<mrow><msubsup><mi>AI</mi><mi>f</mi><mi>φ</mi></msubsup><mo>=</mo><mi>A</mi><mi>I</mi><mo>-</mo><msubsup><mi>AI</mi><mrow><mi>m</mi><mi>a</mi></mrow><mi>φ</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001124403920000028.GIF" wi="1718" he="71" /></maths>6)将步骤5)得到的等效流体声波阻抗的数据集合用石油物探专业的地震显示软件进行图形显示,并利用图形显示进行流体识别。 |
地址 |
225009 江苏省扬州市维扬路188号 |