发明名称 一种卡森模式钻井液流变参数的控制方法
摘要 发明涉及一种卡森模式钻井液流变参数的控制方法,首先确定流变参数合理取值范围,通过旋转粘度计测量钻井液不同转速下的读数并转换为相应的剪切速率和剪切应力,对流变参数进行估计,依据其估计值是否在设定的取值范围内,做现场小型试验,确定钻井液调整处理措施。其中流变参数估计基于最小二乘法,通过数学变换建立关于a的一元目标函数P(a),利用一维搜索法求解a值,使P(a)取最小值,进而求出卡森屈服值和卡森粘度估计值。本方法考虑多种因素设定流变参数合理取值范围,通过提高流变参数估计精度,从而对钻井液调整做出合理判断,利用现场小型试验,确定钻井液维护处理措施,节约处理剂用量,降低钻井成本。
申请公布号 CN102434119B 申请公布日期 2014.03.12
申请号 CN201110280632.4 申请日期 2011.09.20
申请人 中国石油化工股份有限公司;中国石油化工股份有限公司华北分公司工程技术研究院 发明人 闫吉曾;郑锋辉;狄伟;邓红琳
分类号 E21B21/08(2006.01)I 主分类号 E21B21/08(2006.01)I
代理机构 郑州睿信知识产权代理有限公司 41119 代理人 陈浩
主权项 1.一种卡森模式钻井液流变参数的控制方法,其特征在于,该方法包括以下步骤:(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定钻井液卡森屈服值τ<sub>c</sub>和卡森粘度η<sub>∞</sub>的合理取值范围;(2)从振动筛前泥浆槽采集钻井液样本,测量并记录N种不同转速Φ下的读数θ,其中N≥3;(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;(4)根据得到的剪切速率γ和剪切应力τ,估计卡森屈服值τ<sub>c</sub>和卡森粘度值η<sub>∞</sub>;得到卡森屈服估计值<img file="FDA0000400386330000016.GIF" wi="45" he="58" />和卡森粘度估计值<img file="FDA0000400386330000015.GIF" wi="53" he="61" />;(5)将卡森屈服估计值<img file="FDA0000400386330000017.GIF" wi="41" he="50" />和卡森粘度估计值<img file="FDA0000400386330000018.GIF" wi="53" he="49" />与步骤(1)确定的取值范围进行比对,判定流变参数是否在合理取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护;所述步骤(4)中卡森屈服值τ<sub>c</sub>和卡森粘度值η<sub>∞</sub>的估计步骤如下:步骤41:对卡森模式流变方程<img file="FDA0000400386330000011.GIF" wi="308" he="112" />进行参数代换,令<img file="FDA0000400386330000012.GIF" wi="203" he="112" /><img file="FDA0000400386330000013.GIF" wi="170" he="110" />基于最小二乘法原理,建立关于变量a的目标函数P(a),<![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>-</mo><msup><mi>a</mi><mn>2</mn></msup><mo>-</mo><msup><mi>b</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><mn>2</mn><mi>ab</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub></msqrt><msup><mo>)</mo><mn>2</mn></msup></mrow><mo>,</mo></mrow></math>]]></maths>并用a来表示b值:<![CDATA[<math><mrow><mi>b</mi><mo>=</mo><mi>b</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>+</mo><msqrt><msup><mi>B</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>-</mo><mn>4</mn><mi>A</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>C</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msqrt></mrow><mrow><mn>2</mn><mi>A</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中:τ<sub>c</sub>为卡森屈服值;η<sub>∞</sub>为卡森粘度;N为样本容量,N≥3;<img file="FDA0000400386330000022.GIF" wi="343" he="140" />简记为S<sub>m</sub>=S(m);<img file="FDA0000400386330000023.GIF" wi="369" he="140" />简记为T<sub>m</sub>=T(m),m用于幂指数的简化表示可取值为0,0.5,1,1.5或2;<img file="FDA0000400386330000024.GIF" wi="552" he="77" />B(a)=3a<sup>2</sup>S<sub>0.5</sub>S<sub>2</sub>-3a<sup>2</sup>S<sub>1</sub>S<sub>1.5</sub>-T<sub>0.5</sub>S<sub>2</sub>+T<sub>1</sub>S<sub>1.5</sub>;C(a)=Na<sup>3</sup>S<sub>2</sub>-a<sup>3</sup>S<sub>0.5</sub>S<sub>1.5</sub>-aT<sub>0</sub>S<sub>2</sub>+aT<sub>0.5</sub>S<sub>1.5</sub>;步骤42:确定目标函数P(a)的搜索区间[a<sub>1</sub>,a<sub>2</sub>];<![CDATA[<math><mrow><msub><mi>a</mi><mn>1</mn></msub><mo>=</mo><mi>min</mi><mrow><mo>(</mo><mfrac><mrow><msqrt><msub><mi>&tau;</mi><mi>i</mi></msub><msub><mi>&gamma;</mi><mi>j</mi></msub></msqrt><mo>-</mo><msqrt><msub><mi>&tau;</mi><mi>j</mi></msub><msub><mi>&gamma;</mi><mi>i</mi></msub></msqrt></mrow><mrow><msqrt><msub><mi>&gamma;</mi><mi>j</mi></msub></msqrt><mo>-</mo><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub></msqrt></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mi>a</mi><mn>2</mn></msub><mo>=</mo><mi>min</mi><mrow><mo>(</mo><mfrac><mrow><msqrt><msub><mi>&tau;</mi><mi>i</mi></msub><msub><mi>&gamma;</mi><mi>j</mi></msub></msqrt><mo>-</mo><msqrt><msub><mi>&tau;</mi><mi>j</mi></msub><msub><mi>&gamma;</mi><mi>i</mi></msub></msqrt></mrow><mrow><msqrt><msub><mi>&gamma;</mi><mi>j</mi></msub></msqrt><mo>-</mo><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub></msqrt></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中:i和j,下标,且j&gt;i;步骤43:根据目标函数P(a)的搜索区间[a<sub>1</sub>,a<sub>2</sub>],采用一维搜索法求解a值,使目标函数P(a)取最小值;步骤44:依据求得的P(a)最小值,估计卡森屈服值τ<sub>c</sub>和卡森粘度η<sub>∞</sub>,卡森屈服值估计值记为<img file="FDA0000400386330000026.GIF" wi="78" he="72" />卡森粘度估计值记为<img file="FDA0000400386330000027.GIF" wi="90" he="71" />计算公式分别为:<![CDATA[<math><mrow><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mi>c</mi></msub><mo>=</mo><msup><mi>a</mi><mn>2</mn></msup><mo>,</mo><msub><mover><mi>&eta;</mi><mo>^</mo></mover><mo>&infin;</mo></msub><mo>=</mo><msup><mi>b</mi><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>所述步骤(2)中用旋转粘度计测量钻井液样本在N种不同转速Φ下的读数θ。
地址 100728 北京市朝阳区朝阳门北大街22号