CN102434119A - 一种卡森模式钻井液流变参数的控制方法 - Google Patents

一种卡森模式钻井液流变参数的控制方法 Download PDF

Info

Publication number
CN102434119A
CN102434119A CN2011102806324A CN201110280632A CN102434119A CN 102434119 A CN102434119 A CN 102434119A CN 2011102806324 A CN2011102806324 A CN 2011102806324A CN 201110280632 A CN201110280632 A CN 201110280632A CN 102434119 A CN102434119 A CN 102434119A
Authority
CN
China
Prior art keywords
casson
value
drilling fluid
gamma
viscosity
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.)
Granted
Application number
CN2011102806324A
Other languages
English (en)
Other versions
CN102434119B (zh
Inventor
闫吉曾
郑锋辉
狄伟
邓红琳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Engineering Technology Research Institute Of North China Branch china Petroleum & Chemical Corp
China Petroleum and Chemical Corp
Original Assignee
Engineering Technology Research Institute Of North China Branch china Petroleum & Chemical Corp
China Petroleum and Chemical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Engineering Technology Research Institute Of North China Branch china Petroleum & Chemical Corp, China Petroleum and Chemical Corp filed Critical Engineering Technology Research Institute Of North China Branch china Petroleum & Chemical Corp
Priority to CN201110280632.4A priority Critical patent/CN102434119B/zh
Publication of CN102434119A publication Critical patent/CN102434119A/zh
Application granted granted Critical
Publication of CN102434119B publication Critical patent/CN102434119B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种卡森模式钻井液流变参数的控制方法,首先确定流变参数合理取值范围,通过旋转粘度计测量钻井液不同转速下的读数并转换为相应的剪切速率和剪切应力,对流变参数进行估计,依据其估计值是否在设定的取值范围内,做现场小型试验,确定钻井液调整处理措施。其中流变参数估计基于最小二乘法,通过数学变换建立关于a的一元目标函数P(a),利用一维搜索法求解a值,使P(a)取最小值,进而求出卡森屈服值和卡森粘度估计值。本方法考虑多种因素设定流变参数合理取值范围,通过提高流变参数估计精度,从而对钻井液调整做出合理判断,利用现场小型试验,确定钻井液维护处理措施,节约处理剂用量,降低钻井成本。

Description

一种卡森模式钻井液流变参数的控制方法
技术领域
本发明涉及石油天然气钻井领域,涉及一种符合卡森模式的钻井液流变参数的控制方法。
背景技术
在石油天然气钻井领域,钻井液被比喻为钻井的的血液,它与钻井速度、钻井成本及钻井参数紧密相关,钻井液的性能是钻井成败的重要因素之一。一般情况下,钻井液属于非牛顿流体,符合卡森流变模式的称为卡森流体,卡森模式(Casson Model)由卡森于1959年提出,是常用的流变模式之一,可适用于各种类型的钻井液(见文献《钻井液工艺学》第76页,鄢捷年主编,中国石油大学出版社,2000),流变方程为
Figure BDA0000092835170000011
其流变参数τc称为卡森动切力或卡森屈服值,η称为极限高剪切粘度或卡森粘度,常用来近似表示钻井液在钻头喷嘴处紊流状态下的流动阻力,因此也称为水眼粘度。这两个参数的合理控制对于钻井液有效携岩、安全优快钻井有重要意义(见文献《钻井手册(甲方)上册》第773页,《钻井手册(甲方)》编写组,石油工业出版社,1990),也与钻井液的压降、摩阻计算及优化施工参数紧密相关,根据钻井工艺的要求,确定流变参数的合理范围是十分重要的,因此行业标准推荐τc值范围为τc≥4Pa,η值范围为3~10mPa·s(见文献《SY/T 5234-2004优选参数钻井基本方法及应用》)。目前,现场调整钻井液流变参数的方法存在一些不足与缺陷,主要体现在两个方面:(1)流变参数估计精度不高,导致控制不当,增加工人劳动量,浪费钻井液材料,增加钻井成本;(2)通过泥浆大班感官眼看手摸的办法判断钻井液性能,过多的依赖于泥浆大班的经验,缺乏合理的手段和程序。
目前,流变参数的估计,通常是采用旋转粘度计测量入井流体不同转速下的读数,根据这些实测数据选择一定的方法进行估计,得到相应的卡森屈服值估计值
Figure BDA0000092835170000021
和卡森粘度估计值
Figure BDA0000092835170000022
估计方法主要有三种,一是直接估计法,二是线性回归法,另一个是非线性回归法。
直接估计法主要是采用旋转粘度计100转读数、300转读数和600转读数进行估计,主要缺陷是计算点单一,仅仅代表某一剪切速率下的流变参数,而钻井液在井内的循环过程从钻杆、钻铤、环空至返出井口,剪切速率变化范围很大(见文献《钻井液工艺学》第68页,鄢捷年主编,中国石油大学出版社,2000),所以这种估计方法不能准确反映各剪切速率下的流变性。
为克服直接估计法的不足,另一种现场常用方法是线性回归法,对卡森模式流变方程
Figure BDA0000092835170000023
进行线性化处理,转换为线性方程,Y=a+bX(其中:
Figure BDA0000092835170000024
Figure BDA0000092835170000026
Figure BDA0000092835170000027
),根据线性回归法相关公式计算出τc和η,转换后新方程出现方差异性,不再满足Gauss-Markov假定,参数估计不具有方差一致最小性,结果造成流变参数的估计精度大大降低,对钻井液调整产生“误导”。
非线性回归法由于计算精度要高于线性回归法,越来越受到工程技术人员的重视。目前,利用非线性回归法估计流变参数,通常通过LMF算法(Levenberg-Marquardt-Fletcher)求解,属于二元函数的极值问题,须设置初始值,进行判断转移次数较多,涉及大量矩阵运算,见文献《最优化原理与方法》(薛毅著,北京工业大学出版社,2001);文献《钻井液卡森参数非线性最小二乘估计新算法》(鲁港李晓光等,石油学报2008年第3期)提出了一种方法,主要是求解二元非线性方程组,然后判断根的合理性。但是,这两种方法都过于复杂,编程步骤过多,判断转移次数多,计算量大,给实际应用造成较大困难。
发明内容
本发明的目的是提供一种卡森模式钻井液流变参数的控制方法,以解决现有控制方法中计算量大、方法复杂,不适于现场实时操作的问题。
为实现上述目的,本发明的一种卡森模式钻井液流变参数的控制方法包括以下步骤:
(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定钻井液卡森屈服值τc和卡森粘度η的合理取值范围;
(2)从振动筛前泥浆槽采集钻井液样本,测量并记录N种不同转速Φ下的读数θ,其中N≥3;
(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
(4)根据得到的剪切速率γ和剪切应力τ,估计卡森屈服值τc和卡森粘度值η,得到卡森屈服估计值和卡森粘度估计值
Figure BDA0000092835170000032
(5)将卡森屈服估计值
Figure BDA0000092835170000033
和卡森粘度估计值
Figure BDA0000092835170000034
与步骤(1)确定的取值范围进行比对,判定流变参数是否在合理取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护。
进一步的,所述步骤(2)中用旋转粘度计测量钻井液样本在N种不同转速Φ下的读数θ。
进一步的,所述步骤(3)中是依据旋转粘度计仪器常数,将不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ。
进一步的,所述步骤(4)中卡森屈服值τc和卡森粘度值η的估计步骤如下:
步骤41:对卡森模式流变方程
Figure BDA0000092835170000041
进行参数代换,令
Figure BDA0000092835170000043
基于最小二乘法原理,建立关于变量a的目标函数P(a), P ( a ) = Σ i = 1 N ( τ i - a 2 - b 2 ( a ) γ i - 2 ab ( a ) γ i ) 2 , 并用a来表示b值:
b = b ( a ) = - B ( a ) + B 2 ( a ) - 4 A ( a ) C ( a ) 2 A ( a ) ,
其中:τc为卡森屈服值;η为卡森粘度;N为样本容量,N≥3;
Figure BDA0000092835170000046
简记为Sm=S(m);
Figure BDA0000092835170000047
简记为Tm=T(m),m用于幂指数的简化表示可取值为0,0.5,1,1.5或2;
Figure BDA0000092835170000048
B(a)=3a2S0.5S2-3a2S1S1.5-T0.5S2+T1S1.5
C(a)=Na3S2-a3S0.5S1.5-aT0S2+aT0.5S1.5
步骤42:确定目标函数P(a)的搜索区间[a1,a2];
a 1 = min ( τ i γ j - τ j γ i γ j - γ i ) , a 2 = max ( τ i γ j - τ j γ i γ j - γ i ) ; 其中:i和j,下标,且j>i;
步骤43:根据目标函数P(a)的搜索区间[a1,a2],采用一维搜索法求解a值,使目标函数P(a)取最小值;
步骤44:依据求得的P(a)最小值,估计卡森屈服值τc和卡森粘度η,卡森屈服值估计值记为
Figure BDA0000092835170000051
卡森粘度估计值记为
Figure BDA0000092835170000052
计算公式分别为: τ ^ c = a 2 , η ^ ∞ = b 2 .
进一步的,所述步骤43中一维搜索法为黄金分割法、二分法、分数法、进退法或插值法。
进一步的,所述步骤(4)可通过EXCEL平台实现。
所述步骤(5)中若流变参数估计值未在确定的取值范围内,则根据钻井液基本配方,现场做小型试验,确定钻井液材料加量,对钻井液进行处理。
本发明的卡森模式钻井液流变参数的控制方法通过判定流变参数估计值是否在其合理区间内,依据现场小型试验,确定钻井液处理措施。其流变参数估计根据卡森模式流变方程基于最小二乘法,所得卡森模式流变参数是最小二乘意义下的最优估计。本发明的优点在于,将求解二元函数极值问题,通过一系列数学变换转化为求一元函数P(a)在区间[a1,a2]上的最小值问题,计算过程不涉及矩阵及导数运算,不用设初始值,缩短了搜索区间,减少了计算量,加快了收敛速度,算法清晰、简单、稳定,易于实现,为准确判定流变参数调控范围提供了依据,节约了钻井液材料,减少了钻井成本。
附图说明
图1是流变参数的控制方法流程图;
图2是卡森模式流变参数的估计步骤流程图;
图3是目标函数P(a)搜索区间和最小值示意图;
图4是黄金分割法计算流程图;
图5是目标函数曲线、目标函数最小值和搜索区间图;
图6是流变参数在EXCEL 2003中计算过程图。
具体实施方式
卡森模式钻井液流变参数的控制方法,具体步骤为:
步骤1:根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定钻井液卡森屈服值τc和卡森粘度η的合理取值范围;
步骤2:在振动筛前泥浆槽采集钻井液样本,用旋转粘度计测量其N种不同转速Φ下的读数θ,并记录N种不同转速Φ及相对应的读数θ;
步骤3:依据旋转粘度计仪器常数,将步骤2得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
步骤4:根据步骤3数据,估计卡森屈服值τc和卡森粘度值η,得到卡森屈服估计值
Figure BDA0000092835170000061
和卡森粘度估计值
步骤5:根据步骤1数据、步骤4数据,将卡森屈服估计值和卡森粘度估计值与步骤(1)确定的取值范围进行比对,判定流变参数是否在合理取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护。
一、卡森模式钻井液流变参数估计:
已知某钻井液属于卡森型流体,通过旋转粘度计测量数据利用本发明方法求卡森屈服值τc和卡森粘度η
如图2所示,具体步骤如下:
步骤一:通过常用的范氏35型六速旋转粘度计测量钻井液6种不同转速对应的读数,如下表1所示:
表1
  转速Φ(r/min)   3   6   100   200   300   600
  读数数θ(格)   5   6   28   45   58   95
步骤二:把步骤一得到的数据,根据旋转粘度计仪器常数,转换为相应的剪切应力和剪切速率,转换关系式(见非专利文献《SY/T 5480-2007固井设计规范》)如下:
τ = 0.511 θ γ = 1.7023 Φ
转换结果如下:
τ1=2.555,τ2=3.066,τ3=14.308,τ4=22.995,τ5=29.638,τ6=48.545;
γ1=5.1069,γ2=10.2138,γ3=170.23,γ4=340.46,γ5=510.69,γ6=1021.38;
步骤三:对卡森模式流变方程
Figure BDA0000092835170000072
进行参数代换
Figure BDA0000092835170000073
Figure BDA0000092835170000074
根据其流变方程,基于最小二乘法原理,利用步骤二得到的数据,建立变量a的目标函数P(a): P ( a ) = Σ i = 1 6 ( τ i - a 2 - b 2 ( a ) γ i - 2 ab ( a ) γ i ) 2 , 并用a来表示b, b = b ( a ) = - B ( a ) + B 2 ( a ) - 4 A ( a ) C ( a ) 2 A ( a ) ;
其中:τc,卡森屈服值,Pa;
η,卡森粘度,Pa·s;
A(a)=605272275.1383a;
B(a)=72244589.3940a2-170308871.3301;
C(a)=3868797.3730a3-25325334.7335a;
步骤四:依据步骤二得到的数据,计算得到目标函数P(a)的搜索区间[a1,a2]=[1.0919,1.9086],搜索区间的计算公式为:
Figure BDA0000092835170000081
Figure BDA0000092835170000082
其中:i和j为下标,且j>i
步骤五:根据步骤四得到的数据,采用一维搜索法中的黄金分割法(0.618法)求解a值,使目标函数P(a)取最小值。
Figure BDA0000092835170000083
取ε=0.001,这个精度完全满足工程需要,k为循环次数,具体求解步骤如下,见图4:
Step1:计算x1=a1+(1-λ)(a2-a1),x2=a1+λ(a2-a1);
Step2:判断a2-a1<ε是否成立,若成立,则
Figure BDA0000092835170000084
循环结束;若不成立,则转Step3;
Step3:比较P(x1)-P(x2)>0是否成立,若成立,则令
a1=x1,x1=x2,x2=a1+λ(a2-a1);若不成立,则令
a2=x2,x2=x1,x1=a1+(1-λ)(a2-a1);转Step2;
具体计算过程见下表2:
表2
  k   a1   x1   x2   a2   P1(x1)   P2(x2)  |a2-a1|<ε   a   Pmin
  1   1.0919   1.4039   1.5967   1.9086   50.4170   77.2675   0.8167   1.5003   63.8423
  2   1.0919   1.2847   1.4039   1.5967   361.4299   50.4170   0.5047   1.3443   205.9235
  3   1.2847   1.4039   1.4775   1.5967   50.4170   4.3426   0.3119   1.4407   27.3798
  4   1.4039   1.4775   1.5230   1.5967   4.3426   14.2720   0.1928   1.5003   9.3073
  5   1.4039   1.4494   1.4775   1.5230   11.8655   4.3426   0.1192   1.4635   8.1041
  6   1.4494   1.4775   1.4949   1.5230   4.3426   5.1228   0.0736   1.4862   4.7327
  7   1.4494   1.4668   1.4775   1.4949   5.8807   4.3426   0.0455   1.4721   5.1117
  8   1.4668   1.4775   1.4842   1.4949   4.3426   4.1758   0.0281   1.4808   4.2592
  9   1.4775   1.4842   1.4883   1.4949   4.1758   4.3628   0.0174   1.4862   4.2693
  10   1.4775   1.4816   1.4842   1.4883   4.1704   4.1758     0.0107   1.4829   4.1731
  11   1.4775   1.4801   1.4816   1.4842   4.2096   4.1704     0.0066   1.4808   4.1900
  12   1.4801   1.4816   1.4826   1.4842   4.1704   4.1625     0.0041   1.4821   4.1664
  13   1.4816   1.4826   1.4832   1.4842   4.1625   4.1638     0.0025   1.4829   4.1631
  14   1.4816   1.4822   1.4826   1.4832   4.1640   4.1625     0.0016   1.4824   4.1633
  15   1.4822   1.4826   1.4828   1.4832   4.1625   4.1624     0.0010   1.4827   4.1624
从上表看出,循环15次,就达到了设定的计算精度。计算得到a=1.4827时,函数P(a)取最小值,最小值为Pmin=4.1624,参见图5;
步骤六:依据步骤五得到的数据,计算卡森屈服值和卡森粘度,将计算结果保留4位小数,得到,
Figure BDA0000092835170000091
η ^ ∞ = b 2 = ( - B ( a ) + B 2 ( a ) - 4 A ( a ) C ( a ) 2 A ( a ) ) 2 = 0.0300 ;
可以看到,将求解二元函数 P ( a , b ) = Σ i = 1 6 ( τ i - a 2 - b 2 γ i - 2 ab γ i ) 2 极值问题,通过数学变换转化为求一元函数 P ( a ) = Σ i = 1 6 ( τ i - a 2 - b 2 ( a ) γ i - 2 ab ( a ) γ i ) 2 在区间[a1,a2]上的最小值问题,见图3,由于函数P(a)在区间[a1,a2]=[1.0919,1.9086]上是单峰函数,见图5,利用黄金分割法(0.618法)很容易求解出卡森模式流变参数估计值
整个计算过程没有涉及到矩阵和导数运算,也没有设初始值,搜索区间也大大缩短,计算量大大减少,仅仅循环15次就达到了工程需要的计算精度,算法清晰、简单、稳定,非常易于实现,在常用的EXCEL软件就可实现,图6是本估计法在EXCEL 2003中的计算过程图。
当然,在求解函数P(a)在区间[1.0919,1.9086]上的最小值时,也可用二分法、分数法(Fibonacci法)、进退法或插值法等方法。
二、卡森模式钻井液流变参数的控制方法
大牛地气田某口水平井三级井身结构,一开采用Φ17 1/2″钻头下Φ13 3/8″表层套管,二开采用Φ12 1/4″钻头下Φ9 5/8″技术套管,三开采用Φ8 1/2″钻头,先期裸眼完井,钻井液采用天然高分子钻井液,这种钻井液体系属于卡森流体。钻至井深500.16m后下表层套管,钻至井深2876.29m时已进入山西组,钻井液流变参数的控制具体实施步骤见图1。
第一步,根据二开井眼尺寸为311.2mm,井径扩大率大约为2%,地层属砂泥岩地层,固相颗粒较小,当量直径最大为几个毫米,因此确定卡森屈服值范围为0.5≤τc≤3.0Pa和卡森粘度范围为3.0≤η≤10mPa·s;
第二步,在振动筛前泥浆槽内天然高分子钻井液样本,通过六速旋转粘度计测得钻井液数据如下表3所示:
表3
  六速旋转粘度计转速Φ(r/min)   3   6   100   200   300   600
  天然高分子钻井液读数θ(格)   2   6   13   19   26   40
第三步,转换为相应的剪切速率和剪切应力,转换结果如表4所示:
表4
 剪切速率γ(s-1)   5.1069   10.2138   170.23   340.46   510.69   1021.38
 剪切应力τ(Pa)   1.022   3.066   6.643   9.709   13.286   20.44
第四步,通过本发明方法估计卡森屈服值和卡森粘度,结果为τc=1.45Pa,η=10.95mPa·s;
第五步,根据第四步计算结果,判定卡森粘度偏高,不利于优快钻井,需要对天然高分子钻井液进行调整,在泥浆槽取钻井液样本做现场小型试验,通过实验确定需要加入烧碱、NAT20和IND30。
第六步,调整钻井液,加入烧碱是50kg、NAT20是375kg和IND30是100kg,充分循环后,在井深2800.5m测钻井液性能,测量数据如表5:
表5
  六速旋转粘度计转速Φ(r/min)   3   6   100   200   300   600
  天然高分子钻井液读数θ(格)   2   4   10   16   22   34
用本方法再次对流变参数进行估计,结果为
Figure BDA0000092835170000111
Figure BDA0000092835170000112
在设定合理范围内,则只需对钻井液进行正常维护。
最后所应说明的是:以上实施例仅用以说明而非限定本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应当理解;依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。

Claims (7)

1.一种卡森模式钻井液流变参数的控制方法,其特征在于,该方法包括以下步骤:
(1)根据井眼尺寸、井径扩大率、环空间隙、固相颗粒和泥浆泵排量大小,确定钻井液卡森屈服值τc和卡森粘度η的合理取值范围;
(2)从振动筛前泥浆槽采集钻井液样本,测量并记录N种不同转速Φ下的读数θ,其中N≥3;
(3)将得到的不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ;
(4)根据得到的剪切速率γ和剪切应力τ,估计卡森屈服值τc和卡森粘度值η,得到卡森屈服估计值
Figure FDA0000092835160000011
和卡森粘度估计值
Figure FDA0000092835160000012
(5)将卡森屈服估计值
Figure FDA0000092835160000013
和卡森粘度估计值
Figure FDA0000092835160000014
与步骤(1)确定的取值范围进行比对,判定流变参数是否在合理取值范围内,若是未在确定的取值范围内,则对钻井液进行调整,并重复步骤(2)至(5),直至钻井液流变参数在确定的取值范围内;若是在确定的取值范围内,则对钻井液进行正常维护。
2.根据权利要求1所述的一种卡森模式钻井液流变参数的控制方法,其特征在于:所述步骤(2)中用旋转粘度计测量钻井液样本在N种不同转速Φ下的读数θ。
3.根据权利要求2所述的一种卡森模式钻井液流变参数的控制方法,其特征在于:所述步骤(3)中是依据旋转粘度计仪器常数,将不同转速Φ和对应的读数θ转换为相应的剪切速率γ和剪切应力τ。
4.根据权利要求3所述的一种卡森模式钻井液流变参数的控制方法,其特征在于,所述步骤(4)中卡森屈服值τc和卡森粘度值η的估计步骤如下:
步骤41:对卡森模式流变方程
Figure FDA0000092835160000021
进行参数代换,令
Figure FDA0000092835160000022
基于最小二乘法原理,建立关于变量a的目标函数P(a), P ( a ) = Σ i = 1 N ( τ i - a 2 - b 2 ( a ) γ i - 2 ab ( a ) γ i ) 2 , 并用a来表示b值:
b = b ( a ) = - B ( a ) + B 2 ( a ) - 4 A ( a ) C ( a ) 2 A ( a ) ,
其中:τc为卡森屈服值;η为卡森粘度;N为样本容量,N≥3;
Figure FDA0000092835160000026
简记为Sm=S(m);
Figure FDA0000092835160000027
简记为Tm=T(m),m用于幂指数的简化表示可取值为0,0.5,1,1.5或2;
Figure FDA0000092835160000028
B(a)=3a2S0.5S2-3a2S1S1.5-T0.5S2+T1S1.5
C(a)=Na3S2-a3S0.5S1.5-aT0S2+aT0.5S1.5
步骤42:确定目标函数P(a)的搜索区间[a1,a2];
a 1 = min ( τ i γ j - τ j γ i γ j - γ i ) , a 2 = max ( τ i γ j - τ j γ i γ j - γ i ) ; 其中:i和j,下标,且j>i;
步骤43:根据目标函数P(a)的搜索区间[a1,a2],采用一维搜索法求解a值,使目标函数P(a)取最小值;
步骤44:依据求得的P(a)最小值,估计卡森屈服值τc和卡森粘度η,卡森屈服值估计值记为
Figure FDA00000928351600000211
卡森粘度估计值记为
Figure FDA00000928351600000212
计算公式分别为:
τ ^ c = a 2 , η ^ ∞ = b 2 .
5.根据权利要求4所述的一种卡森模式钻井液流变参数的控制方法,其特征在于:所述步骤43中一维搜索法为黄金分割法、二分法、分数法、进退法或插值法。
6.根据权利要求5所述的一种卡森模式钻井液流变参数的控制方法,其特征在于:所述步骤(4)可通过EXCEL平台实现。
7.根据权利要求1-6中任一项所述的一种卡森模式钻井液流变参数的控制方法,其特征在于:所述步骤(5)中若流变参数估计值未在确定的取值范围内,则根据钻井液基本配方,现场做小型试验,确定钻井液材料加量,对钻井液进行处理。
CN201110280632.4A 2011-09-20 2011-09-20 一种卡森模式钻井液流变参数的控制方法 Expired - Fee Related CN102434119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110280632.4A CN102434119B (zh) 2011-09-20 2011-09-20 一种卡森模式钻井液流变参数的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110280632.4A CN102434119B (zh) 2011-09-20 2011-09-20 一种卡森模式钻井液流变参数的控制方法

Publications (2)

Publication Number Publication Date
CN102434119A true CN102434119A (zh) 2012-05-02
CN102434119B CN102434119B (zh) 2014-03-12

Family

ID=45982444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110280632.4A Expired - Fee Related CN102434119B (zh) 2011-09-20 2011-09-20 一种卡森模式钻井液流变参数的控制方法

Country Status (1)

Country Link
CN (1) CN102434119B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104458501A (zh) * 2014-11-27 2015-03-25 河南理工大学 卡森流体流变参数测定方法
CN105352851A (zh) * 2015-12-07 2016-02-24 中国石油化工股份有限公司 一种入井流体类型的识别方法
CN105370229A (zh) * 2014-08-22 2016-03-02 中国石油化工股份有限公司 钻井液的控制方法
CN110388202A (zh) * 2018-04-13 2019-10-29 中国石油化工股份有限公司 基于高温高压粘度计读数预测的井筒流体流变特性方法
CN112949014A (zh) * 2019-11-26 2021-06-11 中国石油天然气股份有限公司 获取钻井液流变参数的方法、装置、存储介质以及处理器
CN114077812A (zh) * 2020-08-14 2022-02-22 中国石油化工股份有限公司 一种用于设计钻井液流变参数的方法
CN114439500A (zh) * 2021-12-16 2022-05-06 山东大学 基于随钻测试的tbm穿越不良地质智能掘进系统与方法
WO2022133761A1 (zh) * 2020-12-23 2022-06-30 四川维泰科创石油设备制造有限公司 一种弯管钻井液流变性的现场测量方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1096100A (zh) * 1993-05-26 1994-12-07 西安石油学院 一种远距离测量液体流变性能装置
WO2001002832A1 (en) * 1999-07-06 2001-01-11 Sofitech N.V. Modelling the rheological behaviour of drilling fluids as a function of pressure and temperature
CN101532962A (zh) * 2009-04-14 2009-09-16 中国地质大学(武汉) 一种钻井液低温性能测试方法及装置
CN102071926A (zh) * 2010-12-02 2011-05-25 中国石油集团钻井工程技术研究院 一种全井段环空压力测量方法、装置及控制方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1096100A (zh) * 1993-05-26 1994-12-07 西安石油学院 一种远距离测量液体流变性能装置
WO2001002832A1 (en) * 1999-07-06 2001-01-11 Sofitech N.V. Modelling the rheological behaviour of drilling fluids as a function of pressure and temperature
CN101532962A (zh) * 2009-04-14 2009-09-16 中国地质大学(武汉) 一种钻井液低温性能测试方法及装置
CN102071926A (zh) * 2010-12-02 2011-05-25 中国石油集团钻井工程技术研究院 一种全井段环空压力测量方法、装置及控制方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张冠军等: "大牛地气田无土相钻井液流变模式优选", 《探矿工程(岩土钻掘工程)》, vol. 37, no. 7, 31 December 2010 (2010-12-31), pages 22 - 24 *
鲁港等: "钻井液卡森模式流变参数非线性最小二乘估计新算法", 《石油学报》, vol. 29, no. 3, 31 May 2008 (2008-05-31), pages 470 - 474 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105370229A (zh) * 2014-08-22 2016-03-02 中国石油化工股份有限公司 钻井液的控制方法
CN105370229B (zh) * 2014-08-22 2019-01-04 中国石油化工股份有限公司 钻井液的控制方法
CN104458501A (zh) * 2014-11-27 2015-03-25 河南理工大学 卡森流体流变参数测定方法
CN105352851A (zh) * 2015-12-07 2016-02-24 中国石油化工股份有限公司 一种入井流体类型的识别方法
CN105352851B (zh) * 2015-12-07 2017-11-10 中国石油化工股份有限公司 一种入井流体类型的识别方法
CN110388202A (zh) * 2018-04-13 2019-10-29 中国石油化工股份有限公司 基于高温高压粘度计读数预测的井筒流体流变特性方法
CN112949014A (zh) * 2019-11-26 2021-06-11 中国石油天然气股份有限公司 获取钻井液流变参数的方法、装置、存储介质以及处理器
CN114077812A (zh) * 2020-08-14 2022-02-22 中国石油化工股份有限公司 一种用于设计钻井液流变参数的方法
CN114077812B (zh) * 2020-08-14 2022-08-12 中国石油化工股份有限公司 一种用于设计钻井液流变参数的方法
WO2022133761A1 (zh) * 2020-12-23 2022-06-30 四川维泰科创石油设备制造有限公司 一种弯管钻井液流变性的现场测量方法
CN114439500A (zh) * 2021-12-16 2022-05-06 山东大学 基于随钻测试的tbm穿越不良地质智能掘进系统与方法
CN114439500B (zh) * 2021-12-16 2023-09-05 山东大学 基于随钻测试的tbm穿越不良地质智能掘进系统与方法

Also Published As

Publication number Publication date
CN102434119B (zh) 2014-03-12

Similar Documents

Publication Publication Date Title
CN102434119B (zh) 一种卡森模式钻井液流变参数的控制方法
CN103226641B (zh) 深水气液两相流循环温度压力耦合计算方法
CN102425386B (zh) 一种符合幂律模式的钻井液流变参数控制方法
CN103590824A (zh) 经过多段压裂改造后的致密气藏水平井的产能计算方法
CN110656915B (zh) 一种页岩气多段压裂水平井多工作制度产能预测方法
Cheng et al. Mechanics of hydraulic fracturing
RU2717019C1 (ru) Способ вывода на режим скважины, пробуренной в естественно трещиноватом пласте
CN110543745A (zh) 一种热采井水泥环封固完整性的评价方法
Ettehadi et al. Functional and practical analytical pressure surges model through Herschel Bulkley fluids
CN111353205A (zh) 用于致密气藏产水气井地层压力和动态产能的计算方法
Zolfaghari et al. Blade-shaped hydraulic fracture driven by a turbulent fluid in an impermeable rock
CN104866681B (zh) 高温高压油气斜井关井过程中温度压力数值模拟方法
CN110593856B (zh) 一种固井安全作业密度窗口测定方法
US11262282B2 (en) Analyzing drilling fluid rheology at a drilling site
CN115293066B (zh) 考虑地层渗流传热效应的气井温度场计算方法
Zhou et al. A novel prediction model for rheological properties of drilling fluids at HTHP conditions and its evaluation
CN115081352A (zh) 深水高温高压疏松砂岩端部脱砂工艺设计优化方法及装置
Olayiwola et al. Optimizing economic number of transverse fractures in horizontal well: A systematic design for maximum tight gas recovery
Darwesh et al. Kicks controlling techniques efficiency in term of time
Chowdhury et al. Mechanical borehole stability analysis using phase2
Chen et al. Practical shale gas decline analysis in Changning shale gas play-a variable pressure drop study
Zhang et al. A newly developed rate analysis method for a single shale gas well
CN111914494B (zh) 一种致密储层水平井压后压裂液返排率预测方法及系统
CN117703355A (zh) 基于在线监测的井筒压力控制智能决策方法、设备及可读存储介质
CN117933108A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140312

Termination date: 20170920