发明名称 水合物热分解引起的喷发灾害的分析与防治方法
摘要 本发明涉及一种水合物热分解引起的喷发灾害的分析与防治方法,从相关的热力学及地理参数入手,通过水合物热分解引起喷发破坏时的静力学平衡条件,建立水合物分解范围与孔压、土层力学强度特性的临界关系式,通过得到的临界关系式,获得实际地层中水合物热分解引起喷发灾害发生的评估条件并提供预防措施。
申请公布号 CN102493787A 申请公布日期 2012.06.13
申请号 CN201110389682.6 申请日期 2011.11.30
申请人 中国海洋石油总公司;中海石油研究中心;中国科学院力学研究所 发明人 李清平;鲁晓兵;姚海元;王淑云;张旭辉;刘乐乐;候金林;王丽勤;王世圣
分类号 E21B43/00(2006.01)I;E21B43/30(2006.01)I 主分类号 E21B43/00(2006.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;关畅
主权项 1.一种水合物热分解引起的喷发灾害的分析与防治方法,其特征在于,1)对水合物热分解引起的喷发灾害的分析包括如下:①水合物热分解区域半径的扩展分析水合物沉积层中热源处首先达到水合物相平衡温度,发生相变,形成分解区域与未分解区域,其界面称为相变阵面,分解区域组成为:水、气体、土骨架,未分解区域组成为:水合物、土骨架;随着温度的升高,分解区域不断扩展,分解区域的温度取决于分解区域的平均密度ρ<sub>1</sub>、平均比热C<sub>1</sub>、平均热传导系数K<sub>1</sub>、加热温度T<sub>h</sub>、初始温度T<sub>0</sub>、相变潜热ΔH、相平衡温度T<sub>e</sub>,ρ<sub>1</sub>=ε<sub>0</sub>ρ<sub>g</sub>+(1-ε<sub>0</sub>)ρ<sub>s</sub>,C<sub>1</sub>=ε<sub>0</sub>C<sub>g</sub>+(1-ε<sub>0</sub>)C<sub>s</sub>,K<sub>1</sub>=ε<sub>0</sub>K<sub>g</sub>+(1-ε<sub>0</sub>)K<sub>s</sub>;未分解区域的温度取决于未分解区域的平均密度ρ<sub>2</sub>、平均比热C<sub>2</sub>、平均热传导系数K<sub>2</sub>,ρ<sub>2</sub>=ε<sub>0</sub>ρ<sub>h</sub>+(1-ε<sub>0</sub>)ρ<sub>s</sub>,C<sub>2</sub>=ε<sub>0</sub>C<sub>h</sub>+(1-ε<sub>0</sub>)C<sub>s</sub>,K<sub>2</sub>=ε<sub>0</sub>K<sub>h</sub>+(1-ε<sub>0</sub>)K<sub>s</sub>,其中:ε<sub>0</sub>为水合物的含量;ρ<sub>h</sub>为水合物密度,ρ<sub>s</sub>为土骨架密度,ρ<sub>g</sub>气体密度,C<sub>h</sub>为水合物比热,C<sub>s</sub>为土骨架比热,C<sub>g</sub>为气体比热,K<sub>h</sub>为水合物热传导系数,K<sub>s</sub>为土骨架热传导系数,K<sub>g</sub>为气体热传导系数;假定各个区域热传导物理参数平均,且相变潜热为常数,写出两个区域内的热学表达式为:分解区域:<img file="FDA0000114260040000011.GIF" wi="1316" he="436" />未分解区域:<img file="FDA0000114260040000012.GIF" wi="1421" he="569" />求解方程(1)得到<img file="FDA0000114260040000013.GIF" wi="263" he="69" />其中,ξ<sub>e</sub>是取决于水合物沉积层各组分热学参数,由关系式<maths num="0001"><![CDATA[<math><mrow><mfrac><mrow><msub><mrow><mn>4</mn><mi>&theta;</mi></mrow><mi>h</mi></msub><msub><mi>K</mi><mn>2</mn></msub></mrow><mrow><msub><mi>&rho;</mi><mi>h</mi></msub><mi>&Delta;</mi><msub><mi>H&epsiv;</mi><mn>0</mn></msub><msub><mi>&kappa;</mi><mn>1</mn></msub></mrow></mfrac><mrow><mo>(</mo><mfrac><msub><mi>&kappa;</mi><mn>1</mn></msub><msub><mi>&kappa;</mi><mn>2</mn></msub></mfrac><mfrac><mrow><msub><mi>&upsi;</mi><mi>e</mi></msub><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>4</mn><mfrac><msub><mi>&kappa;</mi><mn>1</mn></msub><msub><mi>&kappa;</mi><mn>2</mn></msub></mfrac><msub><mi>&xi;</mi><mi>e</mi></msub><mo>)</mo></mrow></mrow><mrow><msqrt><msubsup><mi>&pi;&xi;</mi><mi>e</mi><mo>&prime;</mo></msubsup></msqrt><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>erf</mi><mrow><mo>(</mo><mn>2</mn><msqrt><mfrac><msub><mi>&kappa;</mi><mn>1</mn></msub><msub><mi>&kappa;</mi><mn>2</mn></msub></mfrac><msub><mi>&xi;</mi><mi>e</mi></msub></msqrt><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mrow><mn>4</mn><mi>&theta;</mi></mrow><mi>h</mi></msub><msub><mi>K</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&rho;</mi><mi>h</mi></msub><msub><mi>&Delta;H&epsiv;</mi><mn>0</mn></msub><msub><mi>&kappa;</mi><mn>1</mn></msub></mrow></mfrac><mrow><mo>(</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>&upsi;</mi><mi>e</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>4</mn><msub><mi>&xi;</mi><mi>e</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msqrt><msub><mi>&pi;&xi;</mi><mi>e</mi></msub></msqrt><mo>&CenterDot;</mo><mi>erf</mi><mrow><mo>(</mo><mn>2</mn><msqrt><msub><mi>&xi;</mi><mi>e</mi></msub></msqrt><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>=</mo><mn>1</mn></mrow></math>]]></maths>确定,其中,υ<sub>e</sub>=T<sub>e</sub>-T<sub>0</sub>,θ<sub>h</sub>=T<sub>h</sub>-T<sub>0</sub>,T<sub>1</sub>为分解区域内平均温度,T<sub>2</sub>为未分解区域内平均温度,X<sub>e</sub>为相变阵面的位置,t是时间,κ为热扩散系数,有<img file="FDA0000114260040000022.GIF" wi="215" he="119" /><img file="FDA0000114260040000023.GIF" wi="227" he="119" />分解区域半径R与水合物相变阵面的位置X<sub>e</sub>对应,因此<img file="FDA0000114260040000024.GIF" wi="241" he="69" />可记做:<maths num="0002"><![CDATA[<math><mrow><mi>R</mi><mo>=</mo><msqrt><msub><mi>&xi;</mi><mi>e</mi></msub><msub><mi>&kappa;</mi><mn>1</mn></msub><mi>t</mi></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>②水合物热分解临界喷发分析当水合物分解区域半径达到一临界值时,记为<img file="FDA0000114260040000026.GIF" wi="310" he="69" />气体压力将突破上覆层而发生喷发,喷发所对应的水合物热分解区域半径R<sub>cri</sub>和破坏的特征尺度r由如下参数决定:孔隙气体压力p<sub>g</sub>,上方大气压力p<sub>0</sub>,上覆土体的平均比重<img file="FDA0000114260040000027.GIF" wi="85" he="46" />分解区域上方的土体强度τ<sub>f</sub>,热源至上覆层表面的高度h;因此,分解区域临界半径R<sub>cri</sub>和破坏的特征尺度为r可分别写成:<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>R</mi><mi>cir</mi></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>g</mi><mo>,</mo><msub><mi>p</mi><mi>g</mi></msub><mo>,</mo><msub><mi>&tau;</mi><mi>f</mi></msub><mo>,</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>r</mi><mo>=</mo><mi>s</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>g</mi><mo>,</mo><msub><mi>p</mi><mi>g</mi></msub><mo>,</mo><msub><mi>&tau;</mi><mi>f</mi></msub><mo>,</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>选取<img file="FDA0000114260040000029.GIF" wi="60" he="46" />和h作为单位,式(3)可化为无量纲关系:<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>R</mi><mi>cir</mi></msub><mi>h</mi></mfrac><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>p</mi><mi>g</mi></msub><mo>-</mo><msub><mi>p</mi><mn>0</mn></msub></mrow><mrow><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>gh</mi></mrow></mfrac><mo>,</mo><mfrac><msub><mi>&tau;</mi><mi>f</mi></msub><mrow><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>gh</mi></mrow></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mfrac><mi>r</mi><mi>h</mi></mfrac><mo>=</mo><mi>s</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>p</mi><mi>g</mi></msub><mo>-</mo><msub><mi>p</mi><mn>0</mn></msub></mrow><mrow><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>gh</mi></mrow></mfrac><mo>,</mo><mfrac><msub><mi>&tau;</mi><mi>f</mi></msub><mrow><mover><mi>&rho;</mi><mo>&OverBar;</mo></mover><mi>gh</mi></mrow></mfrac><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>已知<img file="FDA00001142600400000211.GIF" wi="163" he="115" />和<img file="FDA00001142600400000212.GIF" wi="115" he="120" />进行喷发临界条件的静力学分析,则可以判断是否会发生喷发,若发生喷发,则可确定水合物分解区域的特征半径<img file="FDA00001142600400000213.GIF" wi="80" he="108" />及破坏的特征尺度<img file="FDA00001142600400000214.GIF" wi="59" he="98" />2)水合物热分解引起的喷发灾害的防治包括如下:根据步骤1)的分析,从如下方面进行防治:①在水合物开发工程中控制热源加热温度、加热效率和加热时间,并预测分解区域半径的扩展,确定分解区域的孔隙气体压力的发展,使得孔隙压力不足以对上覆层造成破坏;②在获得的现场沉积层和上覆层的力学参数、孔隙和裂隙结构、水合物的横向、纵向合成与分布基础上,选取合适的开采方式和开采井布置方法,并对沉积层或上覆层薄弱区域进行加固处理。
地址 100010 北京市东城区朝阳门北大街25号