CN104895560B - 一种深水测试井筒压力、温度场模拟及水合物预测方法 - Google Patents
一种深水测试井筒压力、温度场模拟及水合物预测方法 Download PDFInfo
- Publication number
- CN104895560B CN104895560B CN201510333445.6A CN201510333445A CN104895560B CN 104895560 B CN104895560 B CN 104895560B CN 201510333445 A CN201510333445 A CN 201510333445A CN 104895560 B CN104895560 B CN 104895560B
- Authority
- CN
- China
- Prior art keywords
- temperature
- wellbore
- pressure
- flow
- parameter
- 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.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种深水测试井筒压力、温度场模拟及水合物预测方法,其步骤:准备参数;对井筒温度场和压力场进行初始化;利用气藏渗流数学模型求解地层流量和井底压力分布信息;将地层流量和井底压力分布信息输入井筒流动模型内进行迭代,求解出井筒压力剖面并判断井筒压力是否收敛;计算热物性参数,并计算地层段和海水段的总传热系数;将井筒压力和总传热系数代入井筒温度场模型,计算各节点井筒温度剖面;判断井筒温度剖面是否收敛;判断tn是否小于tn+1,若小于则输出井温‑井深曲线以及井筒压力对应的水合物生成温度‑井深曲线,并根据曲线判断是否有水合生成;反之则返回预测下一时刻的井筒压力、温度场模拟及水合物情况。本发明可以广泛在石油勘探过程中应用。
Description
技术领域
本发明涉及一种深水测试井筒压力、温度场及水合物预测方法,特别是关于一种在海水温度影响下井筒稳态传热和非稳态传热下的深水测试井筒压力、温度场模拟及水合物预测方法。
背景技术
目前我国能源开发形式日益严峻,海洋深水区域逐渐成为勘探开发的焦点。但对于深水油气田来说,普遍具有埋藏深,储层高温、高压、渗透性好等特征,即正常情况下,深海油气田中的生产井往往是高产井,但是伴随着高产的条件下,高流量气体结合海水的降温作用以及针对高压力的多级节流,很容易在井筒中出现水合物,堵塞测试管柱和生产管柱。但是目前还没有一套完整的深水测试井筒压力、温度场模拟及水合物预测的技术方法。
井筒压力剖面计算极大地依赖于井筒温度剖面的预测,目前常规方法是假设井筒温度剖面线性化,井筒流体温度与时间无关,这种处理方式对处于稳定生产状态的油气井压力计算是合适的,但是在油气井不稳定测试过程中流量、压力、温度可能都处于不稳态过程,在此情况下沿用常规方法估计井筒温度误差较大,导致压力计算偏差。针对井筒动态压力和温度场模拟技术及水合物预测的技术,国外已经形成有较为成熟商业软件-PIPESIM。该软件中的黑油模型可以对油、气、水三相,气液两相,以及单相流体进行计算井筒温度场、压力剖面模拟。但是通过PIPESIM软件进行井筒温度场、压力剖面模拟预测,其存在以下缺点:1、井筒压力剖面预测是基于假设井筒温度剖面为线性变化的计算上进行的,对于油气井不稳定测试过程中,计算结果偏差大。2、没有考虑试油测试短期过程中井筒与地层的非稳态传热。3、没有考虑海水温度变化。
发明内容
针对上述问题,本发明的目的是提供一种深水测试井筒压力、温度场模拟及水合物预测方法,该方法对试采阶段的长期稳定生产过程,采用解析方法建立井筒稳态传热模型,预测非线性井温剖面;针对试油测试的短期过程,建立井筒非稳态传热模型,预测不同测试制度和时间下的井温剖面;并能计算井筒温度,压力结果可靠,准确度高。
本发明的另一目的是在温度剖面预测基础上,进行压力剖面和水合物预测。
为实现上述目的,本发明采取以下技术方案:一种深水测试井筒压力、温度场模拟及水合物预测方法,其包括以下步骤:1)准备参数:气藏参数、流体物性参数、井身结构参数和海水温度剖面等;2)对井筒进行离散化,并对井筒温度场和压力场进行初始化,初始化井筒温度等于环境温度;3)利用气藏渗流数学模型求解地层流量和井底压力分布信息;4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代,进而求解出井筒压力剖面,并判断井筒压力是否收敛,若收敛则进入步骤5),反之继续求解井筒压力剖面,再判断井筒压力是否收敛;其中,井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井筒压力计算管流模型;5)计算测试井内流体和环空流体在当前温度压力条件下的热物性参数,并计算地层段和海水段的总传热系数;6)将井筒压力和总传热系数代入井筒温度场模型,从井底往井口递推计算各节点井筒温度剖面,并修正地层段和海水段的环空温度差,更新流体热物性参数;各节点是井筒每离散段的节点;7)判断井筒温度剖面是否收敛,若不收敛,则返回步骤6)进行迭代,直至井筒温度剖面收敛;若收敛,则修正井筒的压力和温度;8)判断tn是否小于tn+1,若小于则输出井温-井深曲线以及井筒压力对应的水合物生成温度-井深曲线,并根据曲线判断是否有水合生成;反之若不小于,则返回步骤3),预测下一时刻的井筒压力、温度场模拟及水合物情况。
所述步骤3)中,所述气藏渗流数学模型为:3.1)无限大均质油气藏中心直井Laplace空间压力相应解为:
式中,为井底Laplace空间压力;u=sf(s),其中s为Laplace空间时间变量,f(s)=1;K0、K1分别为修正的零阶和一阶第二类Bessel函数;S为表皮系数;CD为无因次井筒储集系数;3.2)无限大油藏斜井在Laplace空间压力相应解:
式中,hfD为无因次射孔井段长度;s为Laplace空间时间变量;为无因次半径;n为级数系数;hD为无因次储层厚度;zD为井筒深度z的无因次值、为井段中心井筒深度z的无因次值。
所述步骤4)中,所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井筒压力计算管流模型。
所述稳态井筒温度场模型为:
地层段的流体温度Tf为:
式中,z为井筒深度;ψ为倾斜校正系数;
海水段的流体温度Tf表达式:
式中,Fc为焦耳-汤普森效应与动能作用系数;Tf0为井筒入口温度;Tei0为地层温度;z为井筒深度;LR为松弛距离参数;θ为井筒与水平面的夹角;g为重力加速度;ψs为组合参数;定义组合参数ψs为:
所述非稳态井筒温度场模型为:
式中,a为温度随时间变化的衰减系数;Tei为任意深度处的原始地层温度;Cp为油管流体比热,单位为J/kg-K;LR为松弛距离参数;L为垂直井筒总长度,单位为m;gG为地温梯度,单位为K/m;gsG为海水段温度梯度,单位为K/m;t为时间,单位为s。
所述多相流的井筒压力计算管流模型建立方法如下:多相管流模型的基本方程为:
式中,为位差压力梯度;为摩擦压力梯度;为加速度压力梯度;为总压力梯度;其中,总压力梯度的计算步骤如下:(1)计算弗鲁德数NFr、无滑脱持液率EL和L1,L2,L3,L4:
L1=316El 0.302;
L2=0.0009252EL -2.4682;
L3=0.1EL -1.4516;
L4=0.5EL -6.738;
式中,EL为无滑脱持液率;QL为液体体积流量;(2)根据步骤(1)中计算的基本参数,采用Beggs-Brill法判断流型:分离流、间歇流或分散流:当EL<0.01时NFr<L1,EL≥0.01时NFr<L2,则流型为离流;当EL≥0.01,L2<NFr≤L3时,流型为过渡流;当0.01≤EL<0.4时L3<NFr<L1,EL≥0.4时L3<NFr<L4,流型为间歇流;当EL<0.4时NFr≥L1,EL≥0.4时NFr>L4,流型为分散流;(3)持液率计算:倾角为θ的气液两相流动的持液率HL(θ)为:HL(θ)=HL(0)ψ;其中,HL(0)为同样流动参数下,水平流动时的持液率;
式中,ψ为倾斜校正系数,对于垂直管:ψ=1+0.3C,C为系数,与无滑脱持液率、弗鲁德系数和液相速度有关;b、c均为流型常数;(4)计算液相密度ρL:ρL=ρLHL+ρg(1-HL);(5)确定流动阻力系数λ:
式中:
其中,μL为在P、T下的液体粘度;μg为在P、T下的气体粘度;T为温度;(6)计算井筒压力梯度:将步骤(1)~步骤(5)中计算的参数代入多相管流模型的基本方程中计算得到井筒压力梯度。
所述步骤5)中,所述总传热系数采用热阻串联法计算。
本发明由于采取以上技术方案,其具有以下优点:1、本发明采用针对油气井不稳定测试过程中流量、压力、温度可能都处于的不稳态过程,建立了井筒非稳态传热模型,计算井筒温度、压力结果可靠,准确度高。2、本发明为了克服传统井筒温度、压力场预测方法存在的缺陷,同时解决深水产层高流量气体结合海水的降温作用以及针对高压力的多级节流,很容易在井筒中出现水合物的问题,本发明对试采阶段的长期稳定生产过程,采用解析方法建立井筒稳态传热模型,预测的非线性井温剖面;针对试油测试的短期过程,建立井筒非稳态传热模型,预测不同测试制度和时间下的井温剖面。进一步在温度剖面预测基础上,开展压力剖面和水合物预测。本发明可以广泛在石油勘探过程中应用。
附图说明
图1是本发明的整体流程示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供一种深水测试井筒压力、温度场模拟及水合物预测方法,利用气藏渗流数学模型和井筒流动模型,将两个模型耦合求解地层和井筒的瞬态温度、压力剖面,进而实现对井筒压力、温度场模拟及水合物的预测,其包括以下步骤:
1)准备参数:气藏参数、流体物性参数、井身结构参数和海水温度剖面等;
2)对井筒进行离散化,并对井筒温度场和压力场进行初始化,初始化井筒温度等于环境温度;
3)利用气藏渗流数学模型求解地层流量和井底压力分布信息;其中,气藏渗流数学模型为:
利用渗流力学原理和质量守恒定律建立气藏不同完井方式(裸眼、射孔、筛管等)、边界条件(封闭、定压、无限大)和不同井型(直井、斜井、水平井)下的渗流数学模型。其典型模型如下:
3.1)无限大均质油气藏中心直井Laplace空间压力相应解为:
式中,为井底Laplace空间压力;u=sf(s),其中s为Laplace空间时间变量,f(s)=1;K0、K1分别为修正的零阶和一阶第二类Bessel函数;S为表皮系数;CD为无因次井筒储集系数。
3.2)无限大油藏斜井在Laplace空间压力相应解:
式中,hfD为无因次射孔井段长度;s为Laplace空间时间变量;为无因次半径;n为级数系数;hD为无因次储层厚度;zD为井筒深度z的无因次值、为井段中心井筒深度z的无因次值。
4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代,进而求解出井筒压力剖面,并判断井筒压力是否收敛,若收敛则进入步骤5),反之继续求解井筒压力剖面,再判断井筒压力是否收敛;
其中,井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井筒压力计算管流模型;
5)计算测试井内流体和环空流体在当前温度压力条件下的热物性参数,并计算地层段和海水段的总传热系数;其中,总传热系数采用热阻串联法计算;
6)将井筒压力和总传热系数代入井筒温度场模型,从井底往井口递推计算各节点井筒温度剖面,并修正地层段和海水段的环空温度差,更新流体热物性参数;其中,各节点是井筒每离散段的节点;
7)判断井筒温度剖面是否收敛,若不收敛,则返回步骤6)进行迭代,直至井筒温度剖面收敛,一般迭代3~5次;若收敛,则修正井筒的压力和温度;
8)判断tn是否小于tn+1,若小于则输出井温-井深曲线以及井筒压力对应的水合物生成温度-井深曲线,并根据曲线判断是否有水合生成;反之若不小于,则返回步骤3),预测下一时刻的井筒压力、温度场模拟及水合物情况;
其中,在计算第n+1个时间步长时,将用第n个时间步长的井筒温度场和压力场作为初始条件。
上述步骤4)中,稳态和非稳态下的井筒温度场模型建立方法如下:
由于本发明针对试采阶段的长期稳定生产过程,因此采用解析方法建立稳态井筒温度场模型,预测非线性井温剖面;针对试油测试的短期过程,建立非稳态井筒温度场模型,预测不同测试制度和时间下的井温剖面。根据能量守恒定律、焦耳-汤普森效应、井筒径向传热以及海上温度场影响,得到稳态井筒温度场模型和非稳态井筒温度场模型。
4.1)稳态井筒温度场模型
地层段的流体温度Tf为:
式中,z为井筒深度;ψ为倾斜校正系数。
海水段的流体温度Tf表达式:
式中,Fc为焦耳-汤普森效应与动能作用系数;Tf0为井筒入口温度;Tei0为地层温度;z为井筒深度;LR为松弛距离参数;θ为井筒与水平面的夹角;g为重力加速度;ψs为组合参数;
定义组合参数ψs为:
4.2)非稳态井筒温度场模型为:
式(3)~(7)中,a为温度随时间变化的衰减系数;Tei为任意深度处的原始地层温度;Cp为油管流体比热,单位为J/kg-K;LR为松弛距离参数;L为垂直井筒总长度,单位为m;gG为地温梯度,单位为K/m;gsG为海水段温度梯度,单位为K/m;t为时间,单位为s。
上述步骤4)中,多相流的井筒压力计算管流模型建立方法如下:
由于海上油气田一般是大斜度井或水平井,同时结合多种管流计算模型的适应性条件,本发明优选Beggs和Brill方法作为井筒压力计算模型。本发明中Beggs和Brill根据在长为15m、直径为25.4mm和38.1mm的管子中,用空气和水进行实验的基础上提出了处理所有范围的多相流及任何角度的管流的普适化相关式。其中,Beggs和Brill多相管流模型的基本方程、流型判别条件以及主要参数确定如下:
根据机械能量守恒定律,气液两相气液两相管流的压力降消耗于三个方面:位差、摩擦和加速度,即多相管流模型的基本方程为:
式中各项参数的解释如表1所示。
表1式(8)中各项参数的解释表
其中,总压力梯度的计算步骤如下:
(1)计算弗鲁德数NFr、无滑脱持液率EL和L1,L2,L3,L4:
L1=316El 0.302;
L2=0.0009252EL -2.4682;
L3=0.1EL -1.4516;
L4=0.5EL -6.738;
式中,EL为无滑脱持液率;QL为液体体积流量;
(2)根据步骤(1)中计算的基本参数,采用Beggs-Brill法判断流型:分离流、间歇流或分散流。
当EL<0.01时NFr<L1,EL≥0.01时NFr<L2,则流型为离流;
当EL≥0.01,L2<NFr≤L3时,流型为过渡流;
当0.01≤EL<0.4时L3<NFr<L1,EL≥0.4时L3<NFr<L4,流型为间歇流;
当EL<0.4时NFr≥L1,EL≥0.4时NFr>L4,流型为分散流。
(3)持液率计算
倾角为θ的气液两相流动的持液率HL(θ)为:
HL(θ)=HL(0)ψ (9)
式(9)中,HL(0)为同样流动参数下,水平流动时的持液率。
式(9)中,ψ为倾斜校正系数,对于垂直管:ψ=1+0.3C,C为系数,与无滑脱持液率、弗鲁德系数和液相速度有关;b、c均为流型常数;
(4)计算液相密度ρL:
ρL=ρLHL+ρg(1-HL) (11)
(5)确定流动阻力系数λ:
式中:
其中,μL为在P、T下的液体粘度;μg为在P、T下的气体粘度,T为温度。
(6)计算井筒压力梯度:将步骤(1)~步骤(5)中计算的参数代入式(8)中计算得到井筒压力梯度。
上述各实施例仅用于说明本发明,各步骤都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别步骤进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (7)
1.一种深水测试井筒压力、温度场模拟及水合物预测方法,其包括以下步骤:
1)准备参数:气藏参数、流体物性参数、井身结构参数和海水温度剖面;
2)对井筒进行离散化,并对井筒温度场和压力场进行初始化,初始化井筒温度等于环境温度;
3)利用气藏渗流数学模型求解地层流量和井底压力分布信息;
4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代,进而求解出井筒压力剖面,并判断井筒压力是否收敛,若收敛则进入步骤5),反之继续求解井筒压力剖面,再判断井筒压力是否收敛;其中,井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井筒压力计算管流模型;
5)计算测试井内流体和环空流体在当前温度压力条件下的热物性参数,并计算地层段和海水段的总传热系数;
6)将井筒压力和总传热系数代入井筒温度场模型,从井底往井口递推计算各节点井筒温度剖面,并修正地层段和海水段的环空温度差,更新流体热物性参数;各节点是井筒每离散段的节点;
7)判断井筒温度剖面是否收敛,若不收敛,则返回步骤6)进行迭代,直至井筒温度剖面收敛;若收敛,则修正井筒的压力和温度;
8)判断tn是否小于tn+1,若小于则输出井温-井深曲线以及井筒压力对应的水合物生成温度-井深曲线,并根据曲线判断是否有水合生成;反之若不小于,则返回步骤3),预测下一时刻的井筒压力、温度场模拟及水合物情况。
2.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述步骤3)中,所述气藏渗流数学模型为:
3.1)无限大均质油气藏中心直井Laplace空间压力相应解为:
式中,为井底Laplace空间压力;u=sf(s),其中s为Laplace空间时间变量,f(s)=1;K0、K1分别为修正的零阶和一阶第二类Bessel函数;S为表皮系数;CD为无因次井筒储集系数;
3.2)无限大油藏斜井在Laplace空间压力相应解:
式中,hfD为无因次射孔井段长度;s为Laplace空间时间变量;为无因次半径;n为级数系数;hD为无因次储层厚度;zD为井筒深度z的无因次值、为井段中心井筒深度z的无因次值。
3.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述步骤4)中,所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井筒压力计算管流模型。
4.如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述稳态井筒温度场模型为:
地层段的流体温度Tf为:
式中,z为井筒深度;ψ为倾斜校正系数;
海水段的流体温度Tf表达式:
式中,Fc为焦耳-汤普森效应与动能作用系数;Tf0为井筒入口温度;Tei0为地层温度;z为井筒深度;LR为松弛距离参数;θ为井筒与水平面的夹角;g为重力加速度;ψs为组合参数;
定义组合参数ψs为:
5.如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述非稳态井筒温度场模型为:
式中,a为温度随时间变化的衰减系数;Tei为任意深度处的原始地层温度;Cp为油管流体比热,单位为J/kg-K;LR为松弛距离参数;L为垂直井筒总长度,单位为m;gG为地温梯度,单位为K/m;gsG为海水段温度梯度,单位为K/m;t为时间,单位为s。
6.如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述多相流的井筒压力计算管流模型建立方法如下:
多相管流模型的基本方程为:
式中,为位差压力梯度;为摩擦压力梯度;为加速度压力梯度;为总压力梯度;
其中,总压力梯度的计算步骤如下:
(1)计算弗鲁德数NFr、无滑脱持液率EL和L1,L2,L3,L4:
L1=316El 0.302;
L2=0.0009252EL -2.4682;
L3=0.1EL -1.4516;
L4=0.5EL -6.738;
式中,EL为无滑脱持液率;QL为液体体积流量;
(2)根据步骤(1)中计算的基本参数,采用Beggs-Brill法判断流型:分离流、间歇流或分散流:
当EL<0.01时NFr<L1,EL≥0.01时NFr<L2,则流型为离流;
当EL≥0.01,L2<NFr≤L3时,流型为过渡流;
当0.01≤EL<0.4时L3<NFr<L1,EL≥0.4时L3<NFr<L4,流型为间歇流;
当EL<0.4时NFr≥L1,EL≥0.4时NFr>L4,流型为分散流;
(3)持液率计算
倾角为θ的气液两相流动的持液率HL(θ)为:
HL(θ)=HL(0)ψ,
其中,HL(0)为同样流动参数下,水平流动时的持液率;
式中,ψ为倾斜校正系数,对于垂直管:ψ=1+0.3C,C为系数,与无滑脱持液率、弗鲁德系数和液相速度有关;b、c均为流型常数;
(4)计算液相密度ρL:
ρL=ρLHL+ρg(1-HL);
(5)确定流动阻力系数λ:
式中:
其中,μL为在P、T下的液体粘度;μg为在P、T下的气体粘度;T为温度;
(6)计算井筒压力梯度:将步骤(1)~步骤(5)中计算的参数代入多相管流模型的基本方程中计算得到井筒压力梯度。
7.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法,其特征在于:所述步骤5)中,所述总传热系数采用热阻串联法计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510333445.6A CN104895560B (zh) | 2015-06-16 | 2015-06-16 | 一种深水测试井筒压力、温度场模拟及水合物预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510333445.6A CN104895560B (zh) | 2015-06-16 | 2015-06-16 | 一种深水测试井筒压力、温度场模拟及水合物预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104895560A CN104895560A (zh) | 2015-09-09 |
CN104895560B true CN104895560B (zh) | 2017-06-23 |
Family
ID=54028475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510333445.6A Active CN104895560B (zh) | 2015-06-16 | 2015-06-16 | 一种深水测试井筒压力、温度场模拟及水合物预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104895560B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107130955B (zh) * | 2016-02-26 | 2020-12-11 | 中国石油化工股份有限公司 | 井底流压的确定方法及储集体天然能量的确定方法 |
CN107313759B (zh) * | 2016-04-26 | 2019-08-30 | 中国石油化工股份有限公司 | 低渗稠油油藏直井热采压力分布预测方法及系统 |
CN106968667A (zh) * | 2017-04-18 | 2017-07-21 | 中国海洋石油总公司 | 一种温度场预测方法及装置 |
CN106951666B (zh) * | 2017-05-03 | 2020-03-27 | 西南石油大学 | 一种海洋天然气水合物层钻井井筒温度场计算方法 |
CN109538185B (zh) * | 2018-10-30 | 2022-04-22 | 中国海洋石油集团有限公司 | 一种耦合温度场条件下的多层套管井井筒完整性分析模型 |
CN109723970B (zh) * | 2019-01-04 | 2019-11-08 | 中山大学 | 天然气水合物浆体输运的监测与管理系统及方法 |
CN110057488B (zh) * | 2019-04-08 | 2020-12-08 | 浙江省水利河口研究院 | 涌潮冲击圆柱桩的最大压强及其垂向分布的计算方法 |
CN110321647B (zh) * | 2019-07-10 | 2023-08-11 | 合肥迪斯贝能源科技有限公司 | 一种确定垂直井多层合采光纤测温确定每层油产量方法 |
CN110566164B (zh) * | 2019-10-08 | 2021-09-10 | 长江大学 | 一种海上测试管柱最优产量的确认方法 |
CN113536212B (zh) * | 2020-04-15 | 2024-04-02 | 中国石油化工股份有限公司 | 一种基于气液两相管流模型的井底压力确定方法 |
CN113944458B (zh) * | 2020-06-30 | 2023-06-30 | 中国石油天然气股份有限公司 | 单井试井分析方法 |
CN111648764B (zh) * | 2020-07-20 | 2021-03-19 | 西南石油大学 | 一种多层气藏井下分布式温度监测产出剖面解释评价方法 |
CN112069688B (zh) * | 2020-09-09 | 2024-07-02 | 中国石油化工股份有限公司 | 一种天然气长输管道内腐蚀模拟分析方法 |
CN115688342B (zh) * | 2022-11-22 | 2023-11-17 | 西南石油大学 | 一种变截面油管井筒非临界流压力温度分布预测方法 |
CN117688283A (zh) * | 2023-12-12 | 2024-03-12 | 西南石油大学 | 一种稠油井掺稀比和油管下深的实验确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011044070A2 (en) * | 2009-10-06 | 2011-04-14 | Schlumberger Canada Limited | Formation testing planning and monitoring |
CN102682195A (zh) * | 2012-03-09 | 2012-09-19 | 中国石油大学(华东) | 半潜式平台瞬态钻井井筒温度计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7027968B2 (en) * | 2002-01-18 | 2006-04-11 | Conocophillips Company | Method for simulating subsea mudlift drilling and well control operations |
US7996199B2 (en) * | 2006-08-07 | 2011-08-09 | Schlumberger Technology Corp | Method and system for pore pressure prediction |
-
2015
- 2015-06-16 CN CN201510333445.6A patent/CN104895560B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011044070A2 (en) * | 2009-10-06 | 2011-04-14 | Schlumberger Canada Limited | Formation testing planning and monitoring |
CN102682195A (zh) * | 2012-03-09 | 2012-09-19 | 中国石油大学(华东) | 半潜式平台瞬态钻井井筒温度计算方法 |
Non-Patent Citations (2)
Title |
---|
深水钻井井筒全瞬态传热特征;宋洵成等;《石油学报》;20110731(第4期);第704-708页 * |
试采过程的气井非线形井筒温度剖面预测;唐海雄等;《大庆石油学院学报》;20100630(第3期);第101-103页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104895560A (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104895560B (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
CN107045671A (zh) | 产水气井积液风险预测方法 | |
Xie et al. | The influence of backpressure on severe slugging in multiphase flow pipeline-riser systems | |
Sun et al. | Effects of phase transition on gas kick migration in deepwater horizontal drilling | |
Wang et al. | Numerical simulation of two phase flow in wellbores by means of drift flux model and pressure based method | |
Park et al. | Study on severe slugging in an S-shaped riser: small-scale experiments compared with simulations | |
CN106968667A (zh) | 一种温度场预测方法及装置 | |
Li et al. | Modeling of wellbore dynamics of a CO2 injector during transient well shut-in and start-up operations | |
Mao et al. | Temperature prediction model of gas wells for deep-water production in South China Sea | |
CN102682195B (zh) | 半潜式平台瞬态钻井井筒温度计算方法 | |
CN103726815B (zh) | 一种co2驱采出井井筒流态确定及参数优化方法 | |
EP3707345B1 (en) | Determining wellbore leak crossflow rate between formations in an injection well | |
Zhou et al. | Optimization of methods for liquid loading prediction in deep condensate gas wells | |
Jiayan et al. | Effects of variable thermophysical properties of water on the heat extraction of an enhanced geothermal system: A numerical case study | |
CN117307111A (zh) | 页岩气水平井返排阶段的油嘴制度的确定方法及装置 | |
CN102108851B (zh) | 水平井注汽方法及系统 | |
Wang et al. | Development and application of wellbore heat transfer model considering variable mass flow | |
CN104866681A (zh) | 高温高压油气斜井关井过程中温度压力数值模拟方法 | |
Zhang et al. | Establishment and analysis of temperature field of riserless mud recovery system | |
CN116738113A (zh) | 预测地热井筒中结垢位置的数学模型的构建方法 | |
CN105184061A (zh) | 产气井温度压力分布的数值模拟方法 | |
Marcolini et al. | Modeling of wellbore flow within geothermal reservoir simulations at field scale | |
Xu et al. | Tubular string characterization in high temperature high pressure oil and gas wells | |
Shulyupin et al. | Steam-water flow in geothermal wells | |
Yao et al. | Case study on diagnosis and identify the degree of bottom hole liquid accumulation in double-branch horizontal wells in PCOC |
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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |