CN102425386B - 一种符合幂律模式的钻井液流变参数控制方法 - Google Patents
一种符合幂律模式的钻井液流变参数控制方法 Download PDFInfo
- Publication number
- CN102425386B CN102425386B CN201110280092.XA CN201110280092A CN102425386B CN 102425386 B CN102425386 B CN 102425386B CN 201110280092 A CN201110280092 A CN 201110280092A CN 102425386 B CN102425386 B CN 102425386B
- Authority
- CN
- China
- Prior art keywords
- drilling fluid
- rheological parameter
- gamma
- power law
- rheological
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Earth Drilling (AREA)
Abstract
本发明涉及一种符合幂律模式的钻井液流变参数控制方法,根据井眼尺寸、井径扩大率、环空间隙、固相颗粒及泥浆泵排量大小确定流变参数合理取值范围,测量钻井液不同转速下的读数并转换为相应的剪切速率和剪切应力,对流变参数进行估计,依据其估计值是否在确定的取值范围内,做现场小型试验,确定钻井液调整处理措施。本发明的方法考虑多种因素设定流变参数合理取值范围,通过提高流变参数估计精度,从而对钻井液调整做出合理判断,利用现场小型试验,确定钻井液维护处理措施,节约处理剂用量,降低钻井成本。
Description
技术领域
本发明属于石油天然气钻井领域,涉及一种符合幂律模式的钻井液流变参数控制方法。
背景技术
在石油天然气钻井领域,钻井液被比喻为钻井的血液,它与钻井速度、钻井成本及钻井参数紧密相关,钻井液的性能是钻井成败的重要因素之一。一般情况下,钻井液属于非牛顿流体,其流变性用流变方程来描述,幂律模式(Power Law Model)是最常用的流变模式之一,流变方程为τ=Kγn,其流性指数n和稠度系数K的合理控制对于钻井液有效携岩、安全优快钻井有重要意义,也与钻井液的压降、摩阻计算及优化施工参数紧密相关,因此,根据钻井工艺的要求,确定流变参数的合理范围是十分重要的,因此行业标准推荐n值范围为0.4~0.7,K值范围为0.01915~0.1915Pa·sn(见非专利文献《SY/T 5234-2004优选参数钻井基本方法及应用》)。目前,钻井现场钻井液流变参数的方法存在一些不足,主要体现在:(1)流变参数估计精度不高,导致控制不当,浪费钻井液材料,增加钻井成本;(2)通过泥浆大班感官眼看手摸的办法判断钻井液性能,过多的依赖于泥浆大班的经验,缺乏合理的手段和程序。
目前,流变参数的估计,通常是采用旋转粘度计测量入井流体不同转速下的读数,根据这些实测数据选择一定的方法进行计算,得到相应的流性指数估计值和稠度系数估计值。现场钻井液技术员主要采用两种方法,一是直接估计法,另一个是线性回归法。
直接估计法主要是采用旋转粘度计300转读数和600转读数进行估计,主要缺陷是计算点单一,仅仅代表某一剪切速率下的流变参数,而钻井液在井内的循环过程从钻杆、钻铤、环空至返出井口,剪切速率变化范围很大(见非专利文献《钻井液工艺学》第68页,鄢捷年主编,中国石油大学出版社,2000),所以这种估计方法不能准确反映各剪切速率下的流变性。
为克服直接估计法的不足,另一种现场普遍采用的方法是线性回归法,主要是对幂律模式流变方程τ=Kγn两边取对数进行线性化处理,转换为线性方程:Y=a+nX(其中:Y=lnτ,a=lnK,X=lnγ),根据线性回归法相关公式估计出流变参数n和K的值,由于对幂律模式流变方程取对数进行线性化转换后,新方程出现方差异性,不再满足Gauss-Markov假定,参数估计不具有方差一致最小性,结果造成流变参数n和K的估计精度大大降低,对钻井液调整产生“误导”。但是,由于其计算简单,这种方法目前应用较为广泛。
目前,学术上还出现了非线性回归法。非线性回归法通常通过LMF算法(Levenberg-Marquardt-Fletcher)求解,或者转换为一元非线性方程求根,但都没有解析公式,必须进行数值计算,这涉及到初始值的选取、导数运算、确定搜索区间等问题,结果造成循环次数多、计算量偏大、算法复杂且不稳定,给现场实际应用造成较大困难,所以这种方法应用很少。
发明内容
本发明的目的是提供一种符合幂律模式的钻井液流变参数现场控制方法,以解决现有的方法利用线性回归法计算存在方差异性及计算精度较低的问题,克服现场仅凭经验进行钻井液流变参数调整的不足。
为实现上述目的,本发明的一种符合幂律模式的钻井液流变参数控制方法步骤如下:
(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定符合幂律模式的钻井液流性指数n和稠度系数K的取值范围;
(2)从振动筛前泥浆槽内采集钻井液样本,测量并记录钻井液样本在N种不同转速Φ下的读数θ,其中N≥3;
(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
(4)根据得到的剪切速率γ和剪切应力τ估计流性指数n和稠度系数K的值;
(5)将估计的流性指数n和稠度系数K的值与确定的取值范围进行比对,判定流变参数是否在取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护。
进一步的,所述步骤(2)中是用旋转粘度计测量钻井液样本在N种不同转速Φ下的读数θ。
进一步的,所述步骤(3)中是依据旋转粘度计仪器常数,将不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ。
进一步的,所述步骤(4)中流性指数n和稠度系数K的值估计步骤如下:
式中:β=(XTGX)-1XTGY,
β表示二维向量;G表示修正矩阵;λ表示修正矩阵指数;
其中:A=[τ1,τ2,…,τN];
进一步的,所述步骤(5)中如果流变参数估计值未在确定的取值范围内,根据钻井液基本配方,现场做小型试验,确定钻井液材料加量,对钻井液进行处理。
附图说明
图1是流变参数的控制方法流程图;
图2是幂律模式流变参数的估计方法流程图;
图3是本发明方法对第一组数据进行流变曲线拟合图;
图4是使用线性回归法对第一组数据进行流变曲线拟合图;
图5是本发明方法对第二组数据进行流变曲线拟合图;
图6是使用线性回归法对第二组数据进行流变曲线拟合图。
具体实施方式
符合幂律模式的钻井液流变参数控制方法如图1所示,步骤如下:
(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定符合幂律模式的钻井液流性指数n和稠度系数K的取值范围;
(2)从振动筛前泥浆槽内采集钻井液样本,测量并记录钻井液样本在N种不同转速Φ下的读数θ,其中N≥3;
(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
(4)根据得到的剪切速率γ和剪切应力τ估计流性指数n和稠度系数K的值;
(6)将估计的流性指数n和稠度系数K的值与确定的取值范围进行比对,判定流变参数是否在取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护。
下面结合实例对本发明进行详细说明。
第一部分主要用于说明一种幂律模式流变参数估计方法的实施步骤,第二部分主要用于分析对比幂律模式流变参数估计方法与线性回归法的估计结果,第三部分主要用于说明现场钻井液调控方法。
一、幂律模式流变参数估计方法
已知某钻井液属于幂律型流体,通过旋转粘度计测量数据,利用本发明方法求流变参数n和K。
如图2所示,具体实施步骤如下:
步骤一:通过常用的范氏35型六速旋转粘度计测量钻井液,6种不同转速对应的读数如表1所示:
表1
转速Φ(r/min) | 3 | 6 | 100 | 200 | 300 | 600 |
读数θ(格) | 2.5 | 3.5 | 15.5 | 24.5 | 31.5 | 49.5 |
步骤二:根据旋转粘度计仪器常数,把步骤一得到的数据转换为相应的剪切应力和剪切速率,转换关系式(见《SY/T 5480-2007固井设计规范》)如下:
转换结果如下:
τ1=1.2775,τ2=1.7885,τ3=7.9205,τ4=12.5195,τ5=16.0965,τ6=25.2945;
γ1=5.1069,γ2=10.2138,γ3=170.23,γ4=340.46,γ5=510.69,γ6=1021.38;
步骤三:把幂律模式流变方程τ=Kγn两边取自然对数,转换为线性方程Y=a+nX,其中,Y=lnτ,a=lnK,X=lnγ;根据步骤二的数据,基于线性回归法,初步估计流变参数n和K的值,并记为和计算公式如下:
所以,修正矩阵
因此,
其中:
根据步骤二和步骤五得到的数据,得到,
A=[τ1,τ2,τ3,τ4,τ5,τ6]=[1.2775 1.7885 7.9205 12.5195 16.0965 25.2945]
根据流变参数估计结果,绘制流变曲线,见图3和图4。其中,图3是幂律模式流变参数估计方法估计结果绘制的流变曲线,图4是根据线性回归法估计结果绘制的流变曲线。从图3和图4可直观看出,图3较图4实测数据均更贴近拟合曲线。
二、幂律模式流变参数估计方法与线性回归法的估计结果比对
用范氏35型六速旋转粘度计测量了五种不同类型的入井流体,6种不同转速对应的读数如表2所示:
表2
转速Φ(r/min) | 3 | 6 | 100 | 200 | 300 | 600 |
样本1读数θ(格) | 3 | 4 | 34 | 60 | 80 | 138 |
样本2读数θ(格) | 3.7 | 3.7 | 14.5 | 22 | 29.4 | 51.8 |
样本3读数θ(格) | 2 | 3 | 13 | 22 | 29 | 47 |
样本4读数θ(格) | 2 | 3 | 40 | 73 | 100 | 176 |
样本5读数θ(格) | 3 | 4 | 58 | 104 | 142 | 243 |
用幂律模式流变参数估计方法和线性回归法分别估计5个样本的流性指数n和稠度系数K,并分别计算评价拟合效果的几个主要指标,包括:残差平方和、残差样本方差和相关系数。计算结果如表3所示:
表3
由表3分别分析评价拟合效果的几个主要标准指标,
从残差平方和来看,用幂律模式流变参数估计方法计算的残差平方和均大大小于用线性回归法的计算值;
从残差样本方差来看,用幂律模式流变参数估计方法计算的残差样本方差均大大小于用线性回归法的计算值;
从相关系数来看,用幂律模式流变参数估计方法计算的相关系数均高于用线性回归法计算值一个数量级。
从评价拟合效果的这几个主要指标可以看出,幂律模式流变参数估计方法的拟合效果大大优于线性回归法,因此,用幂律模式流变参数估计方法估计的幂律模式流变参数值,其精度比线性回归法大幅度提高。
利用流变曲线拟合效果图,也可以直观的看出幂律模式流变参数估计方法优于线性回归法。
图5是根据幂律模式流变参数估计方法的估计结果绘制的流变曲线;图6是根据线性回归法估计结果绘制的流变曲线。从图5和图6可以明显看出,虽然都是利用相同的实测数据进行流变参数的估计,但图5较图6实测数据均更贴近拟合曲线。其中图5中所有6组实测数据均在拟合曲线上,图6中仅有4组实测数据在拟合曲线上。所以,从流变曲线拟合效果图也可直观看出幂律模式流变参数估计方法优于线性回归法。
三、现场钻井液调控方法
镇泾油田某口水平井,三级井身结构,一开采用Φ143/4″钻头下Φ103/4″表层套管,二开采用Φ91/2″钻头下Φ7″技术套管,三开采用Φ6″钻头,先期裸眼完井,钻井液采用钾铵基钻井液,这种钻井液体系属于幂律流体。钻至井深302.5m后下表层套管,钻至井深1446.5m时已进入侏罗系直罗组,钻井液流变参数的控制具体实施步骤见图1。
第一步,根据二开井眼尺寸为241.3mm,井径扩大率大约为3%,因采用PDC钻头,地层属砂泥岩地层,固相颗粒较小,当量直径最大为几个毫米,因此确定流性指数n范围为0.4~0.7之间,稠度系数K范围在0.02~0.19Pa·sn之间;
第二步,在振动筛前泥浆槽内取钾铵基钻井液样本,通过六速旋转粘度计测得钻井液数据如下表4所示:
表4
六速旋转粘度计转速Φ(r/min) | 3 | 6 | 100 | 200 | 300 | 600 |
钾铵基钻井液读数θ(格) | 1 | 2 | 6 | 11 | 17 | 27 |
第三步,转换为相应的剪切速率和剪切应力,转换结果如表5所示:
表5
剪切速率γ(s-1) | 5.1069 | 10.2138 | 170.23 | 340.46 | 510.69 | 1021.38 |
剪切应力τ(Pa) | 0.511 | 1.022 | 3.066 | 5.621 | 8.687 | 13.797 |
第四步,通过本发明方法估计流性指数n和稠度系数K的值,计算结果为n=0.7633,K=0.0699;
第五步,根据第四步计算结果,判定流性指数偏高,未在确定的取值范围内,不利于钻井液携岩和井底清洁,需要对钾铵基钻井液进行调整,在泥浆槽取钻井液样本做现场小型试验,通过实验确定需要加入坂土、钾盐和CMC,提高钻井液切力,降低流性指数。
第六步,调整钻井液,加入坂土2000kg、钾盐500kg和CMC300kg,充分循环后,在井深1600m测钻井液性能,测量数据如表6:
表6
六速旋转粘度计转速Φ(r/min) | 3 | 6 | 100 | 200 | 300 | 600 |
钾铵基钻井液读数θ(格) | 2 | 3 | 7 | 13 | 19 | 29 |
用幂律模式流变参数估计方法再次估计流性指数n=0.6651,稠度系数K=0.1464,均在合理范围内。
最后所应说明的是:以上实施例仅用以说明而非限定本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应当理解;依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种符合幂律模式的钻井液流变参数控制方法,其特征在于,该方法的步骤如下:
(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定符合幂律模式的钻井液流性指数n和稠度系数K的取值范围;
(2)从振动筛前泥浆槽内采集钻井液样本,测量并记录钻井液样本在N种不同转速Φ下的读数θ,其中N≥3;
(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
(4)根据得到的剪切速率γ和剪切应力τ估计流性指数n和稠度系数K的值;
(5)将估计的流性指数n和稠度系数K的值与确定的取值范围进行比对,判定流变参数是否在取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护;
所述步骤(4)中流性指数n和稠度系数K的值估计步骤如下:
式中:β=(XTGX)-1XTGY,
β表示二维向量;G表示修正矩阵;λ表示修正矩阵指数;
其中:A=[τ1,τ2,…,τN];
2.根据权利要求1所述的一种符合幂律模式的钻井液流变参数控制方法,其特征在于,所述步骤(2)中是用旋转粘度计测量钻井液样本在N种不同转速Φ下的读数θ。
3.根据权利要求2所述的一种符合幂律模式的钻井液流变参数控制方法,其特征在于,所述步骤(3)中是依据旋转粘度计仪器常数,将不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ。
5.根据权利要求1-4中任一项所述的一种符合幂律模式的钻井液流变参数控制方法,其特征在于,所述步骤(5)中如果流变参数估计值未在确定的取值范围内,根据钻井液基本配方,现场做小型试验,确定钻井液材料加量,对钻井液进行处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110280092.XA CN102425386B (zh) | 2011-09-20 | 2011-09-20 | 一种符合幂律模式的钻井液流变参数控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110280092.XA CN102425386B (zh) | 2011-09-20 | 2011-09-20 | 一种符合幂律模式的钻井液流变参数控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102425386A CN102425386A (zh) | 2012-04-25 |
CN102425386B true CN102425386B (zh) | 2014-03-26 |
Family
ID=45959411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110280092.XA Expired - Fee Related CN102425386B (zh) | 2011-09-20 | 2011-09-20 | 一种符合幂律模式的钻井液流变参数控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102425386B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105370229B (zh) * | 2014-08-22 | 2019-01-04 | 中国石油化工股份有限公司 | 钻井液的控制方法 |
CN105352851B (zh) * | 2015-12-07 | 2017-11-10 | 中国石油化工股份有限公司 | 一种入井流体类型的识别方法 |
CN110388202A (zh) * | 2018-04-13 | 2019-10-29 | 中国石油化工股份有限公司 | 基于高温高压粘度计读数预测的井筒流体流变特性方法 |
CN109115656A (zh) * | 2018-07-21 | 2019-01-01 | 中铁十八局集团有限公司 | 一种自动测量砂浆稠度的装置及方法 |
CN110773068B (zh) * | 2018-07-25 | 2022-10-11 | 中国石油化工股份有限公司 | 一种确定原油包水乳状液制备过程中机械搅拌条件的方法 |
CN112949014B (zh) * | 2019-11-26 | 2022-11-04 | 中国石油天然气股份有限公司 | 获取钻井液流变参数的方法、装置、存储介质以及处理器 |
CN114077812B (zh) * | 2020-08-14 | 2022-08-12 | 中国石油化工股份有限公司 | 一种用于设计钻井液流变参数的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4726219A (en) * | 1986-02-13 | 1988-02-23 | Atlantic Richfield Company | Method and system for determining fluid pressures in wellbores and tubular conduits |
CN2033048U (zh) * | 1986-08-02 | 1989-02-22 | 中南工业大学 | 一种可测定钻井浆液三种性能的仪器 |
CN102071926A (zh) * | 2010-12-02 | 2011-05-25 | 中国石油集团钻井工程技术研究院 | 一种全井段环空压力测量方法、装置及控制方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100252259A1 (en) * | 2009-04-01 | 2010-10-07 | Horton Robert L | Oil-based hydraulic fracturing fluids and breakers and methods of preparation and use |
-
2011
- 2011-09-20 CN CN201110280092.XA patent/CN102425386B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4726219A (en) * | 1986-02-13 | 1988-02-23 | Atlantic Richfield Company | Method and system for determining fluid pressures in wellbores and tubular conduits |
CN2033048U (zh) * | 1986-08-02 | 1989-02-22 | 中南工业大学 | 一种可测定钻井浆液三种性能的仪器 |
CN102071926A (zh) * | 2010-12-02 | 2011-05-25 | 中国石油集团钻井工程技术研究院 | 一种全井段环空压力测量方法、装置及控制方法和装置 |
Non-Patent Citations (8)
Title |
---|
不同流变模式下钻井液剪切稀释性评价;王平全;《天然气工业》;19971128;第17卷(第06期);43-45 * |
刘东.幂律钻井液流变参数估计的改进算法.《江汉石油职工大学学报》.2008,第21卷(第02期),31-33. |
幂律钻井液流变参数估计的改进算法;刘东;《江汉石油职工大学学报》;20080320;第21卷(第02期);31-33 * |
王平全.不同流变模式下钻井液剪切稀释性评价.《天然气工业》.1997,第17卷(第06期),43-45. |
钻井液幂律模式非线性最小二乘参数估计的新算法;鲁港等;《钻井液与完井液》;20070910;第24卷;65-68 * |
非牛顿流体流变特性的测量和回归分析;韩洪升等;《大庆石油学院学报》;19930910;第27卷(第03期);19-25 * |
韩洪升等.非牛顿流体流变特性的测量和回归分析.《大庆石油学院学报》.1993,第27卷(第03期),19-25. |
鲁港等.钻井液幂律模式非线性最小二乘参数估计的新算法.《钻井液与完井液》.2007,第24卷65-68. |
Also Published As
Publication number | Publication date |
---|---|
CN102425386A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102425386B (zh) | 一种符合幂律模式的钻井液流变参数控制方法 | |
Quinn et al. | Using constant head step tests to determine hydraulic apertures in fractured rock | |
CN102434119B (zh) | 一种卡森模式钻井液流变参数的控制方法 | |
CN106204304B (zh) | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 | |
US10301936B2 (en) | Tight gas formation pressure determination method | |
CN103867174B (zh) | 一种稠油水平井吸汽状况分析方法及系统 | |
CN105672997A (zh) | 钻井液地层漏失量监测方法 | |
CN105178943A (zh) | 一种实时校正井筒压力的方法 | |
Amani et al. | Prediction of rock strength using drilling data and sonic logs | |
AU5279702A (en) | Method for determining real-time downhole mud properties while drilling | |
CN109915128B (zh) | 地层承压能力动态测试方法及固井方法 | |
CN116150695A (zh) | 一种钻井液漏失位置确定及漏失参数计算方法 | |
Nikitin et al. | Influence of capillary pressure on the restoration of the bottomhole zone permeability at the filtrate-oil interfacial phase | |
CN111963149B (zh) | 一种考虑滞地液量增压的压裂后地层压力求取方法 | |
CN105804713A (zh) | 一种快速确定注水井各小层井口注水启动压力的方法 | |
Wang et al. | Non-uniform Temperature Distribution's Impact on Downhole Weight on Bit Measurement (DWOB) and the Novel Compensatory Method | |
RU2651647C1 (ru) | Способ определения параметров ближней зоны пласта | |
US10564083B2 (en) | Analyzing drilling fluid rheology at a drilling site | |
CN113294147B (zh) | 一种考虑重力因素影响的单洞型断溶体储层试井解释方法 | |
Kårstad | Analysis of ballooning effects during drilling of high pressure high temperature wells | |
CN111520132B (zh) | 一种确定地层中洞距离的方法及系统 | |
US11352883B2 (en) | In-situ rheology behavior characterization using data analytics techniques | |
Wang et al. | Non-uniform temperature distribution's impact on downhole weight on bit (DWOB) measurement and the novel compensatory method | |
Wei et al. | Reservoir evaluation technology during underbalanced drilling of horizontal wells in gas reservoirs | |
CN114673494B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140326 Termination date: 20170920 |
|
CF01 | Termination of patent right due to non-payment of annual fee |