试井解释模型识别及参数估计
试井解释报告模板
试井解释报告第一部分试井解释的理论基础以均质油藏压降试井为例详细阐述现代试井解释的方法、步骤(包括参数的计算方法和公式);说明双重孔隙介质油藏、均质油藏垂直裂缝井所包含的流动阶段、流动阶段的近似解、以及各流动阶段的诊断曲线、特种识别曲线和导数曲线的特点并画出示意图。
第二部分试井解释报告一、测试目的确定地层参数,掌握油气藏的动态资料,具体包括以下几个方面:1、确定井筒储存系数C;2、确定地下流体在地层内的流动能力,即渗透率和流动系数。
3、评价井底污染情况4、确定原始地层压力;二、基础数据如图2-2-1、2-2-2、2-2-3所示为油井定产量生产时压力降落数据。
油藏和井的基本参数见表2-2-1。
表2-2-1 油藏和井的基本参数图2-2-1图2-2-2图2-2-3三、解释结果1、常规方法①早期纯井筒储存阶段C=99.136;结果如图2-3-1、2-3-2所示,C=1e-1m3;D②径向流动阶段结果如图2-3-2所示,k=0.358mD;kh=15.732mD·m; s=-0.547图2-3-1图2-3-22、典型曲线拟合C D=400.00;k=0.350mD;kh= 323.676 mD·m; s=-0.600图2-3-3图2-3-4图2-3-53、一致性检验由常规分析方法和图版拟合方法计算的参数值见表2-3-1表2-3-1 结果对比四、结论1、常规分析方法主要以均质各向同性介质油藏的渗流理论为基础,方法的优点是理论完善,原理简单,易于应用。
但也存在不可避免的缺点,如要求测试时间较长,从而影响生产,无法准确估计井筒储存的特性等。
而现代试井解释方法在一定程度上克服的常规方法存在的问题,使得结论更加的精确.2、由拟合结果k=0.350mD可知,该地层的渗透性属于中等。
因为s=-0.600,所以该油井属于超完善井,可能采取了酸化、压裂等增产措施。
兰亭序永和九年,岁在癸丑,暮春之初,会于会稽山阴之兰亭,修禊事也。
模型识别与参数估计问题分析与研究
模型识别与参数估计问题分析与研究一、引言自从上个世纪六十年代非线性系统控制理论的建立以来,模型识别和参数估计一直是系统控制理论中最重要的问题之一,它们是获得模型参数优化和自适应控制算法的基础。
模型识别和参数估计不仅适用于线性和非线性系统的建模,而且适用于金融、经济、生物、环境等领域中各种系统的建模分析。
为了满足实际生产和科学研究的需要,我们需要研究和开发有效的模型识别和参数估计算法。
二、模型识别和参数估计的定义和意义1. 模型识别的定义和意义模型识别是指利用实验数据对系统动态模型的结构和参数进行估计的过程。
它是系统控制理论和自动化技术中的基础问题。
模型识别不仅适用于系统的建模和仿真,而且适用于控制器设计和优化、系统状态估计和诊断、预测和决策等问题。
通过模型识别,可以获得系统的动态模型,进而进行控制器的设计和优化,以使系统达到所要求的性能指标。
模型识别的目的是在实际应用过程中获得性能优异的控制器。
2. 参数估计的定义和意义参数估计是指利用模型识别过程中的实验数据,对系统动态模型中的参数进行估计的过程。
它是模型识别的一个重要环节,对于系统应用性能的提高至关重要。
参数估计的目的是获得准确、可靠的系统参数,进一步进行控制算法的设计优化,使控制器具有更好的鲁棒性和适应性,提高系统的控制性能和稳定性。
三、模型识别和参数估计的方法和算法1. 基于频域分析的方法基于频域分析的方法是指根据系统的输入输出数据,在频域上对系统的特性进行分析,并在分析的基础上估计系统的模型参数。
其中,经典的频域系统识别方法包括最小平方频域辨识、系数与误差状态空间模型识别、ARX模型识别等。
2. 基于时域分析的方法基于时域分析的方法是指利用系统输入输出数据,按时间序列的顺序分析系统特性,并在分析的基础上估计系统的模型参数。
其中,常见的时域方法包括最小二乘法(OLS)估计、脉冲响应辨识、参数自适应控制法等。
3. 基于神经网络的方法基于神经网络的方法是指利用神经网络的学习能力,对系统的输入输出数据进行训练,从而获得系统的模型参数。
第四章05测井
• 2. (双水)体积模型
固体
液体
骨架 骨架
粉砂
干粘土
粘土束 缚水
湿粘土
泥岩
毛管束 可动 缚水 水 烃 自由水
有效孔隙度
总孔隙度
• 3.基本关系式
f
b
S wt
t
b
S wb t
f
S wf t
t
f
b
h
自由水孔隙度; 束缚水孔隙度; 含油气孔隙度; 总孔隙度
第四章 储层参数定量计算
本章主要内容
• 测井解释模型 • 孔隙度的计算 • 含水饱和度的计算 • 渗透率的计算 • 泥质含量的确定 • 测井解释参数的确定
第一节 测井解释模型
一. 模型的基本含义
所谓模型:就是客观事物被认识后,经 过抽象,再组合为易于理解的形象,即 形象化的抽象。模型的建立过程是:
如:纯砂岩水层=砂岩+孔隙(只有水)
水
砂岩 骨架
V=1个单位
1.通式。
考虑一般情况:地层含泥质、孔隙中含油气,那么体 积模型可画为:
含油气 泥质砂 岩体积 模型
砂岩骨架
Vma
泥质 油气
VSVshhsn
水 Sw
V=1
设我们研究的物理参数为X,用
X ,X X X , ,
h
w
ma
sh
分别表
Qv
CEC(1t )G t
t 泥质砂岩的总孔隙度,小数;
G 岩石的平均颗粒密度, g / cm3
CEC 岩石的阳离子交换能力, mmol / g
Qv 岩石的阳离子交换容量, mmol / cm3
聚合物驱试井-于海洋
3、双区复合油藏聚合物驱试井解释模型
(2)水-聚合物(低浓度聚合物-高浓度聚合物)复合模型
rm
注水井
外区
Ⅴ Ⅰ
内区
பைடு நூலகம்
Ⅳ Ⅱ
Ⅲ
物理模型
牛顿-非牛顿双区复合油藏聚合物驱试井模型的双对数典型曲线存在五个明显的流动阶
段:第Ⅰ段是纯井筒储集阶段,压力和压力导数重合;第Ⅱ段是过渡段,描述纯井筒储集阶
段到内区径向流阶段的过渡;第Ⅲ段是内区径向流段,压力导数曲线为0.5的水平直线段; 第Ⅳ段是牛顿流体到非牛顿流体的过渡段;第Ⅴ段是总系统作用阶段,描述牛顿流体和非牛 顿流体区域共同作用的阶段,由于受到内区非牛顿流体性质的影响 ,曲线表现出小幅度上翘。
16
3、双区复合油藏聚合物驱试井解释模型
聚合物溶液(外区)、水(内区)粘度的影响
内区水相粘度越小或者外区聚合物粘度越大,内外区过渡段和两区共同作用 段上翘幅度越大,而外区聚合物溶液的非牛顿特性则导致了两区共同作用段出
现了小幅度上翘。
17
3、双区复合油藏聚合物驱试井解释模型
内区半径rm的影响
内区半径rm影响内区径向流段和内外区过渡段,rm越大,内区径向流段 持续时间越长,内外区过渡段出现的越晚,内外区共同作用段出现的也就越 晚。
Ⅰ Ⅱ Ⅲ
p pwf p S r rw r r r rw
p(r , t ) p i 外边界条件: lim r
Ⅰ段:井筒储集 Ⅱ段:过渡段 Ⅲ段:径向流段
6
2、单层均质油藏聚合物驱试井解释模型
聚驱 vs 水驱
聚合物驱油藏典型曲线和水驱油藏典型曲线在井储段重合,从径向流段开始, 聚合物驱油藏出现上翘,这是由于聚合物溶液的粘弹效应造成的,显示出了类似 幂律流体的典型曲线特征。
测井解释模型
测井解释模型
测井解释模型是石油勘探开发活动中的一种重要工具,它是通过将已测定义的测井数据进行建模、数据分析、地震解释等步骤,利用这些结果来评价油藏的属性、研究其开发成藏的潜力。
它是有效地利用测井资料评价油气藏的特征、解释岩性变化、提供生产预测的重要技术方法。
测井解释模型包括井眼层析和绘图分析、岩芯分析和应力测试、地质参数分析和相关反演分析、重要层系拾取和解释、储层岩性描述、相关地球物理技术和测井技术。
井眼层析和绘图分析是测井解释模型的核心,也是最重要的技术步骤,它的引申目的是评价主要层系储层的性状,成藏潜力,从而可以对油气藏做出可靠的评价和分析。
此外,测井解释模型在钻井工程中也起着重要作用,例如,它可以为钻井设计提供定量的参数,例如地层厚度、地层性质等。
它还可以获得地质情况下油气藏开发中存在的定量参数和集合属性,从而有助于确定有利的钻井方案、进行定向钻井和排层测气判度,从而为油气藏的开发提供科学的依据。
总而言之,测井解释模型是理解油气藏的最佳技术方法之一,它可以将已经测定了的测井解释转换为有效的信息,为油气藏的开发提供了基础的技术支持,对勘探开发活动具有重要意义。
04第4章 试井解释模型20131018
类别 内边界模型 储层模型
外边界模型 流体模型 流量模型
学
位
模型
Modern well test
第四章、试井解释模型及现代试井解释方法
4.1 试井解释模型
4.3 现代试井解释方法
学
4.2 从系统分析看试井解释
位
6
Modern well test
4.2 从系统分析看试井解释
位
油藏
输入信号 产量变化
学
位
17
Modern well test
现代试井解释方法
现代试井解释方法的重要手段之 一是解释图版拟合,或称为样板曲线 解释图版拟合 拟合(Type Curve Match)。 通过图版拟合,可以得到关于油 藏及油井类型、流动阶段等多方面的 信息,还可以算出K、S、C等参数。
18
学
位
Modern well test
Modern well test
学
第四章 试井解释模型及现代 试井解释方法
1
位
Modern well test
第四章、试井解释模型及现代试井解释方法
4.1 试井解释模型
4.3 现代试井解释方法
学
4.2 从系统分析看试井解释
位
2
Modern well test
4.1 试井解释模型
油气藏在岩石类型、物理性质、埋藏深浅、 压力大小、流体种类和组分等方面都各不相同。但 在试井过程中,所呈现的性态却是有限的。 这是因为油气藏只不过像一个精度不太高的反 应器,只当输出讯号的差别足够大时,地层的差异 才能显现出来,试井才能探测得到。此外,所有各 种性态都只由若干个基本“部分”或“部件”组成。具 体来说,试井解释模型由基本模型、内边界条件和 外边界条件三大部分组成,每一大部分在测试的不 同时间起着支配作用。显然,要想得心应手地选择 试井解释模型,就得对组成解释模型的所有各个基 本“部分”或“阶段”,以及它们的特征有清楚的了解 。
测井解释参数的选择
FS aR 对含泥质少的地层 R 假设a、Rmf为常数
n mm
xo
mf xo
与Hingle图类似,则有
1 Fm ( Sxo )n
m Rxo
aRmf
Sxo 1上的任意点:Rmf FmRxo a
作Sxo=1的线,选含泥质少的水层,把资料点
点在交会图上。
Rmf FmRxo a
Rxo
使用条件:泥质中粘土含量少
3、泥浆滤液的矿化度 (1)、查图版 已知温度、电阻率,用P206图4-23的图版
(2)、已知RmfN(24°C),则可用下式计算 Pmf(NaCL,mg/l)(矿化度):
Pmf 10x
单位:( ppm)
x [3.562 lg(RmfN 0.0123] / 0.955
4、求rmf T
P
在24°C及101.325KPa下,则有下式:
i
式中Pi与Ki——第i种离子的矿化度与系数
用(Pwe、T)查出Rw
为了便于计算机处理,采用近似计算法 Rwn=0.0123+3647.54/Pwn0.955 PwN和 RwN :24度时的地层水总矿化度,地层水电 阻率
任意温度下的Rw :Rw= RwN [45.5/ T(oc)+21.5] 2、用自然电位求Rw
RweN大于0.12 RwN=-0.58+10(0.69RweN - 0.24) RweN小于0.12 RwN=(77RweN+5)/(146-337RweN)
5)求任意温度下的Rw Rw=45.5RwN/[T℃]+21.5]
用自然电位求RW的流程图P209图4-34
在P泥浆>P地时,因过滤电位存在,不能用SP求Rw
油藏课件-油藏工程3-7现代试井解释方法
CD
C
2 hCt rw2
2
4.83102 0.1516.151.422103 0.08782
290
第四步:图版拟合。已经算得 CD, 2我90们只需在接近这个数值的 C曲D 线组,即那 CD 一 1组00样板曲线中进行拟合。拟合结果:实测曲线与 该组中的样板曲线相重合,拟合点为(图5-4):
pD 0.78 p 1MPa
Kh
1.842 103 qB(
pD p
)M
➢ 地层系数:
Kh
1.842
103
q
B(
pD p
)M
➢ 有效渗透率:
K
1.842103
qB
h
(
pD p
)M
第七节 现代试井分析方法简介
➢ 由时间拟合值计算储能系数:
hCt
3.6Kh 1
rw2
(tD t
)M
➢ 由曲线拟合值计算井筒存储系数C 和表皮系数S :
2
CD
) 0.80907 ln(CDe2s )]
pD'
dpD d ( tD
)
1 2( tD
)
CD
CD
pD' tD 1
CD
2
第七节 现代试井分析方法简介
3、压降分析方法与步骤
第一步:初拟合 ➢ 绘图:在比例尺寸与图版相同的双对数坐标纸上绘实测压
力和压力导数曲线Δp~t ,Δp't~t
➢ 拟合:实测曲线与图版拟合,找出一条与实测曲线相吻 合的样板曲线 (初拟合) ,并读出其 CDe2S 值;
Ce 2 s
2Ct hrw2
tD 7.2 Kht
CD
C
第七节 现代试井分析方法简介
第九章 气井的现代试井解释方法
第九章气井的现代试井解释方法气井与油井、水井之间的不同点:(1)气体是可压缩的,μ、压缩系数都是压力的函数;真实气体存在Z(偏差子数)(2)气体渗流不符合达西定律。
那么,如何将已有的油井解释图版用于气井呢?新概念“拟压力”(pseudo-pressure)的概念:⎰=P Pdp ZPP02)(μψ(152)注:P o为参考压力点,一般取P o=0MPa.由此可导出与油井形式完全相同的气井气体渗流方程:因而,只要能算出ψ(P),油井的一切解释方法就可用于气井。
第一节 拟压力的计算方法(重点)一般采用最简单的数值积分方法——梯形法计算拟压力: )]()2()2[(212)(111--=-+==∑⎰i i i i ni PP P P ZPZPdp ZPP μμμψ(P n =P, P 0=0)μ,Z 随P 的变化通过实验得出。
第二节 试井解释方法气井解释区别油井解释唯一的不同是:(1)ψ(P )代替油井的P ,ψD (P )代替油井的P D ,或者说把P D 重新定义。
(2)求出的结果S a 是把表皮函数(包含非达西流的影响)以Gringarten 和Bourdet 图版为例加以说明。
一、Gringarten 图版拟合分析)(489.78)(027143.0)(P qTKh P P T T q KhP P fscf sc D D ψψψ∆=∆== (153)式中: ψ(P i )-ψ[P wf (t)] 压降 △ψ(P )=ψ[P ws (△t)-ψ(P wf )] 压恢 q ——气井产量,104m 3/d T t ——气层温度,K ;P sc ——标准状态下的压力1atm=0.101325MPa T sc ——标准状态下的温度,20℃=293.15K(开压) K ——气层渗透率,μm 2 h ——气层厚度,m. t D 与C D 定义同油井。
注意:英制单位下,Tsc=520ºR(=60ºF=15.55℃=288.75K)与法定单位制所规定的值不等同。
实用现代试井解释模型
1 pD = [ln t D + 0.80907 + 2 S ] 2
p′ = D
dpD tD 1 CD tD ⋅ = ⋅ = 0.5 d( tD / CD ) CD 2 tD CD
Copyright: luopei ,Chongqing University of Science and Technology
pwD = pwDi + pwDb
存在外边界影响的 井底压力 无限大油藏 井底压力 边界影响产生的 井底压力
Copyright: luopei ,Chongqing University of Science and Technology
厚德 博学 砺志 笃行 三、 几种常见边界的压力特征
1. 一条直线封闭边界
厚德 博学 砺志 笃行 (四)均质无限大油藏压力压力特征
101 100
10-1
m =1
10-3 10-2 10-1 100
p′ = 0.5 D
101 102
10-2 10-4
tD
均质无限大油藏压力和压力导数曲线的基本特征
Copyright: luopei ,Chongqing University of Science and Technology
1
2
3 4 5
k = 50mD
k = 100mD
储层渗透率对半对数曲线特征的影响
Copyright: luopei ,Chongqing University of Science and Technology
厚德 博学 砺志 笃行
25 P (MPa) 20 15 10 5 0 10-4 10-3 10-2 10-1 100 101 102 t (hr)
常规试井解释方法
常规试井解释方法常规试井是一种在钻完井以后进行的测试方法,旨在评估井中地层的性质和井的产能。
常规试井通常包括测井、射孔和产量测试。
本文将详细介绍常规试井的原理、步骤以及数据的解释和分析方法。
常规试井的原理是利用测井工具测量井中各点的物理参数,并根据这些参数来推断地层的性质。
其中,测井工具通过电、声、密度和放射性等物理信号来测量地层中的电阻率、声波速度、密度和放射性等参数。
这些参数与地层的含油气性、渗透率和孔隙度等特征相关联。
常规试井的步骤通常包括以下几个阶段:油管下入、测井、射孔和产能测试。
首先,油管被下入井中,将测井仪器下放到需要测试的地层段。
测井仪器包括电阻率测井仪、声波测井仪、密度测井仪和放射性测井仪等。
这些工具通过钢丝绳连接,可以测量不同参数。
测井数据会通过电缆传送到地面。
其次,根据测井的数据,可以计算电阻率、声波速度、密度和放射性等地层参数。
其中,电阻率可以推断出地层的含油气性,电阻率低的地层通常是含油气的。
声波速度和密度可以用来估计地层的渗透率和孔隙度。
放射性数据可以帮助确定地层的组成和厚度。
接下来,通过射孔器进行射孔操作。
射孔是指用爆炸、冲击或冲击弹射等方式在井中形成孔洞,以便使地层与井筒直接相连。
射孔有助于增加地层与井筒的接触面积,提高地层的产能和采收率。
最后,进行产能测试。
产能测试的目的是确定井的流体产能,即每天可产出的油或气的数量。
产能测试可以通过油水分离器和测试管,以及计量和记录仪器来完成。
产能测试时,可以通过控制井口压力和流体的流量来测量不同压力下的流体产能。
在解释和分析常规试井数据时,需要综合考虑各个参数的变化趋势和互相之间的关系。
例如,电阻率降低、声波速度增加、密度增加和放射性增加可能表明地层中存在含油气的区域。
而电阻率增加、声波速度降低、密度降低和放射性降低则可能表示地层中存在含水区域。
此外,在解释常规试井数据时还需要结合地质模型和其他地质信息进行综合分析。
例如,通过对比试井数据和岩心样品的分析结果,可以验证常规试井数据的准确性,并对地层进行更详细的描述和解释。
控制系统中的模型识别与参数估计
控制系统中的模型识别与参数估计在控制系统中,模型识别和参数估计起着至关重要的作用。
模型识别是指通过实验或者数学建模的方法,从系统输入和输出数据中推断出系统的动态特性,以便设计控制器。
参数估计是指用已知的模型结构和一定的观测数据,通过参数估计算法来估计模型中未知参数的值。
本文将介绍控制系统中模型识别与参数估计的基本概念和方法,并探讨其在实际应用中的重要性和挑战。
一、模型识别的基本概念与方法模型识别是控制系统设计的关键步骤之一。
它可以通过实验或者数学建模的方法,从系统的输入输出数据中推断出系统的动态特性,进而用于控制器的设计。
在模型识别中,我们通常假设系统具有一定的结构形式,并用参数来描述系统动态特性。
常用的模型识别方法包括系统辨识方法、频域分析方法和基于机器学习的方法。
其中,系统辨识方法是应用最广泛的方法之一。
它基于系统辨识理论,利用输入输出数据来估计系统的参数。
常见的系统辨识方法包括最小二乘法、极大似然估计和递推最小二乘法等。
频域分析方法则通过对系统的频率响应进行分析,推断出系统的动态特性。
基于机器学习的方法主要利用大数据和人工智能算法,从系统的历史数据中学习和推断出系统的模型。
二、参数估计的基本概念与方法参数估计是指在已知系统模型结构的情况下,通过实验或者观测数据来估计模型中未知参数的值。
参数估计对于控制系统的设计和性能分析具有重要意义。
常用的参数估计方法包括最小二乘法、极大似然估计、扩展卡尔曼滤波和粒子滤波等。
最小二乘法是应用最广泛的方法之一,它通过最小化实际输出与模型预测输出之间的误差来估计参数的值。
极大似然估计则是基于统计学原理,通过最大化观测数据的似然函数来估计参数的值。
扩展卡尔曼滤波和粒子滤波是一类递推滤波算法,可以用于非线性系统的参数估计。
三、模型识别与参数估计在实际应用中的重要性和挑战模型识别与参数估计在控制系统的设计和性能分析中具有重要的作用。
它们可以帮助工程师从实际系统中推断出系统的动态特性,并根据推断结果设计出合适的控制器。
试井分析步骤与过程
辅助试井分析的其它信息
• • • • • • 地质研究 测井 岩心 地震 盆地经验与模拟 井生产动态资料
试井解释框图
数据有效性评价
解释模型识别 油藏参数计算
解释模型检验
不相符
相符否? 结束
五、分析结果报告
六、试井设计
七、建立可供下次参考与使用 的原始资料和解释成果档案
八、试井资料的综合应用及网 络应用系统
MDH曲线分析 当关井前注入时间tp远大于关井时间△t时,即 tp>>△t,有: pws=pws(△t=1h)+ mlog(△t) 流动系数的计算 :
Kh/μ=2.121×10^-3 q.B/m
表皮因子的计算: S=1.1513{(pws1h-pwf)/m-log
[K/(μ.φ.Ct.rw2)] -0.9077} 平均地层压力的计算
100
±¼ Ê ä
ѹÁ¦µ¼Êý
¹ Á Ñ ¦ Ñ ¹ Á ¦ µ ¼ Ê ý Í ¼
三、模型参数计算
1、特征直线分析法 2、 典型曲线分析法 3、非线性回归分析
非线性回归分析技术 是通过改变未知的储层参数使模型与测试数据 充分地拟合,即使测试压力与模型压力差的平方和最小。 n 2 测试 i 模型 i i 1
实际试井分析的步骤与过程
西安石油大学 林加恩
目 录
一、数据确认与数据预处理 1、井与储层的基础数据的确认 2、压力时间数据的确认 3、流量史的确认 二、模型识别 三、模型参数计算 四、模型检验 五、分析结果报告 六、试井设计 七、建立可供下次参考与使用的原始资料和解释成果档案 八、试井资料的综合应用及网络应用系统
一、数据确认与数据预处理
1、井与储层的基础数据的确认
测井解释参数选择
5)求任意温度下的Rw Rw=45.5RwN/[T℃]+21.5]
用自然电位求RW的流程图P209图4-34
在P泥浆>P地时,因过滤电位存在,不能用SP求Rw
3、泥浆滤液的矿化度 (1)、查图版 已知温度、电阻率,用P206图4-23的图版
(2)、已知RmfN(24°C),则可用下式计算 Pmf(NaCL,mg/l)(矿化度):
Pm1f 0 x
单位:( ppm)
x[3.562lgR(mf N0.01]2/03.955
4、求rmf T
P
在24°C及101.325KPa下,则有下式:
(1)用R(1)与R(T)之间的关系求最好
(2)用Rm(地层温度下的)求Rmf Rmf=C•Rm1.07
(系数C为泥浆密度有关,可用rm查出)
(3)用Archie公式求Rmf 条件:m、n、a、F准确,在纯水层处Sx0=1
RmfaRxoFmSxona
在纯水层处Rxo=RMSFL
Rmf=Rmfa
(4)交会图法
aR 对含泥质少的地层 FS R 假设a、Rmf为常数
n mm
xo
mf xo
与Hingle图类似,则有
1 Fm( Sxo )n
m Rxo
aRmf
Sxo1上的任意R点mf : FmRxoa
作Sxo=1的线,选含泥质少的水层,把资料点
点在交会图上。
Rm fFmRxoa
Rxo
使用条件:泥质中粘土含量少
R(T2)=R(T)[T1+x]/(T2+x)。 式中x=10-0.340396lgR(1)+0.641 R(1):起始温度T(1)(OF)时水溶液的电阻率 (2)用微电极曲线求
注水井试井模型及其数值解
② 流体 模 型 : 水微 可 压 缩 , 略 毛管 力 和重 力 油 忽
引
言
的影 响 。 ③ 初 始 条 件 : 注 人 开 始 时 整 个 地 层 压 力 为 常 在
注 人 井 试 井 应 用 于 防止 水 窜 、 持 压 力 和 提 高 保 采 收率 采 油 操作 等 方 面 , 与 生 产 井 试 井 同样 可 确 除 定 地 层 渗透 率 、 附 近地 带 储 层 的伤 害 程度 、 均地 井 平
虑 注水 井 周 围水 饱 和 度 分 布 的条件:
内边 界 条 件 :
P : =P o l
S f o “ () : =s 3
法 的主要缺点是它要求预先 知道相对渗透率 曲线。 18 99年 Y h A awl 用 油 藏 模 拟 软 件 模 拟 了 大 e 和 gra使 量 的注水井压力 落差数据 , 出了不需 要相对渗透 提
层 压 力外 , 可 以确定 井 周 围水 的饱 和 度分 布 、 还 油水
量 ; 藏 中 的油气 饱 和 度是 常 数 ; 水 井完 全 穿透 油 油 注
藏 , 匀注 人 。 均
④ 内边 界 条 件 : 虑井 储 效 应 和 表皮 效 应 。 考 ⑤ 外边 界 条 件 : 藏外 边 界 封 闭或 外 边 界定 压 。 油 ⑥ 水驱 油 符 合 BclyIvrt模 型 。 uk —ee t e _ e
维普资讯
1 2
油
气
井
测
试
20 O 2年 8 月
度, m;≯ — 地 层 孔 隙 度 ;p、 — — 油 、 — 0 水 密 度 ;| — 表 皮 系 数 ;S — 原 始 水 饱 和 s — w —
当 sf | 时 ≥ s 其 中 式 中: ——含水率;
现代试井分析理论与解释方法
8)半球面流、球面流 油藏由于存在气顶或者底水,为了防止底水锥进或者气顶气窜,只打开油层顶 部或者底部,油层中的流体类似于从半球体的四面方向流向油层顶部的打开部位, 此时的流动称为“半球形流动”。 如果只在油层中某一部位打开,油层流体从射孔孔眼的上下、左右、前后四面 八方流向孔眼,此时的流动称为“球形流动”。 厚油层局部打开时可以在“早期段”出现“半球形”或者“球形”流动。
哪些数据点呈现直线关系
20世纪50年代至今,都在使用这种半对数分析法,被称为“常规试 井解释方法”。在直角坐标纸上绘制出井底流动压力pwf与开井生产时间t 的对数lgt关系曲线,或在半对数坐标纸上绘制出pwf与开井生产时间t的关 系曲线就得到一条“压力降落曲线”。根据该曲线的斜率m就能计算出流 动系数、流度、渗透率和表皮。
7
9)直角坐标图、半对数坐标图、双对数坐标图 试井常用直角坐标图是指纵坐标为时间轴、横坐标为压力轴,在该坐标系中跟实际 测试的数据绘制出压力随时间的变化关系曲线。 半对数坐标图,是相对于直角坐标图而言,纵坐标为时间的对数轴;绘制压力随时 间对数值的变化关系曲线。用于常规试井解释。 双对数坐标图,是相对于直角坐标图而言,纵坐标为时间的对数轴,横坐标为时间 的对数轴;绘制压力对数值与时间对数值的变化关系曲线。用于现代时间图版拟合解释。
15
8
三、试 井 分 析 方 法
简化地质模型
建立数学模型
分离变量 积分变换等
数学模型求解
不同坐标系
寻找直线规律、拟合点 求取参数
直线段的斜率和截距 K、S、d
9
气井的现代试井解释方法
qTf Kh
lg tp t t
pw2s
pw2f
42.42 C0 psc Tsc
qTf Kh
Kt
(lg Ctrw2
0.9077 0.8686Sa )
22
压力平方法
(p)
P2
压降曲线 pw2f lg t
霍纳曲线
pw2s
lg
tp
t t
MDH曲线 pw2s lg t
pwf lg t
pws
pi
pwf (1h) m
K
lg Ctrw2
0.9077
Sa
1.151
pws
(1h) m
pwf
K
lg Байду номын сангаасCtrw2
0.9077
30
图版拟合分析
pD
0.054286Kh Tsc
q
Tf
pi
i Zi
p psc
Kh 157
qTf
pi
iZi
p
格林加坦图版拟合
K
qiZiTf
157 pih
(
P 2P
(P)
dP
Po Z
气体渗流方程
2 1 1
r2
r
r
3.6
t
2
§1 拟压力的计算
可用最简单的“梯形法”计算拟压力:
P 2P
( p)
dP
P0 Z
n 1 2P
2P
j1 2 [( Z ) j ( Z ) j1]( Pj Pj1 )
(P)
P,MPa
3
§2 试井解释方法
43.04 443.75 78.489 17.1
1.70 104
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
、 一 介
丈
.1 统一归一化理论曲故和实测曲钱
.2 分 l 理通过对子特征的等份化处理, 在消除实测曲线时间轴上不完整性的同时, 将 实测曲线移动到 样本曲 线位置, 但是, 等份化处理的同时必然带来曲 线实际形态的变化, 从 而可能改变区别模式类别的特征, 如图3 、图4 0
以 式中 x. 间 数 最 值, 力 数 数 最 值, 丛 是 间 . 时 对 的 小 Y.压 导 对 的 小 N和 分别 时 上 ,f i 是 i 是 和 力 对 周数 处 实 曲 时 要 采 与论 线 同N N 压 的 数 期 。 理 侧 线 , 求 用 理 曲 相 的, , 和 值。
S E 4 5 P 662
试井解释模型识别及参数估计 的人工神经网络方法
邓远忠 王家宏
〔 袖勘探开发科学研究院 ) 石
陈钦雷
{ 石油大学 r ) 北京 )
摘
要
试井解释模型的正确识别是试井解释的墓础,压力导数的线特征是拼识试井解释棋型的主 要对象。试井解释模型的自动识别对实现自 动试井分析具有重要意义,神经网络方法由于其其 有较张的模式识别能力和抗噪声干扰能力,因而被作为实现自动识别试井解释模型的方法而得 到广泛研究。神经网络能否正确识别实测曲线的试井解释棋型类别, 取决于构造神经网络的方 法。由于试井问月的特殊性,对于点向f方法构造的神经网络,传统的归一化方法不能有效的 将同一模式的实测曲线统一到学习样本的模式域中,因此网络对实侧试井曲线的模式识别能力
CiaeBue 版,A Ds 线的 数 合。 模型的 力导数曲 a g n od 图 V e为曲 参 组 不同 n r - rt t Cz 压 线的总 特征不
同, 线的 曲 早期、 中期和晚期特征分别反映 井筒、 油藏及外边界类型。 在双对数坐标中的压 力导数曲线特征具 有一定的共性: 井筒和表皮的影响使得早期压力导数曲线出现蜂值, 且出现时机及形态只与井筒和表
②比 缩 : =.0"; .I 叼=.0" ) 例 小 对 0+.(一 jN、 0十.(一 / 1 8x x ) . . 1 8 漏 凡 y ,
③分 理 可= x ) 处 : (一-A二一. 、 (一m} ̄一o 段 x ; x x 对= v )v Y ) , j Y }C m ,
分段处理过程中, Yx i 则令广 0 如果曲 n 如果 mYr, o n , 线由 个子部分组成, 则先将曲线 划分成月 个部分, 分别利用上述方法进行处理, 然后再将处理后的数据按照原曲 线的 组合顺 序连接起来, 最后, } ,x 将x ` 1, 通7 } , , 换到[ 1 x= / l ' A A 0 1 区间。
3 试井问厄的特殊性 .
对于点序列神经网络输人矢量, 传统的归一化处理方法除不能有效将实侧曲 线统一到 理论样本曲线的模式域中外, 还常伴随曲 线特征的形变, 使得预处理后曲 线的 模式发生了转
移。
图1 、图2 为点序列的二范数处理结果, 可见归一化处理未能将属于同 一模式的理论曲 线和实测曲线统一起 来, 尤其是不当的预处理还改变了实测曲线的 形态 ( ) 图2 。对于这类 点序列输人矢t, 神经网络将作为不同的模式对待。
前
言
试井解释一般包括试井解释模型识别和试井解释参数计算两个过程, 其中 试井解释模 型的正确识别是试井解释的 基础。目 前试井解释模型的 识别主要由 解释人员根据试井压力导 数特征和油田 储层地质情况综合判断确 定,自 动化识别程度较低。 近年来,由于人工神经网 络 ( N 具有较强的抗噪声干扰能力和非线性映射特性, A N) 其被广泛应用于多种领域处理 各种复杂问题。9年代初期, 0 试井分析领 域开始将神经网 络技术引 入试井分析过程, 到目 前 为止, 神经网络技术仍然主 要用于识别试井解释模型。在有关的研究报道〔 ‘ 均采用压 ’〕 一中, 力导数曲线的点序列作为神经网 络的输人来构造神经网络。为了肪止神经网络出现饱和现象 或者统一实测、 理论曲 线到同 一值域, 文献中利用二范数或者比 例变换等方法对点向 量进行 归一化处理, 然后作为神经网 络的 输人。由 于试井问 题的 复杂性, 实测压力导数曲线不可能 作为u i练神经网络的学习 q 样本, 所以均采用无因次压力导数曲 线作为 训练网 络的 样本。 学习 与样本曲 线相比, 实测压力导数曲 线总存在一定的位移, 而且早期和后期常不完整, 这些差 异使得用理论曲 线点序列训练的 神经网络很难正确识别实测曲 线的 模式类别。这是因为, 用 曲线的点序列作为训练神经网 络的 输人矢量。 传统的各种归一化处理以 及比 例变换方法都难
38一 8
1.—
下了 下石一一一-一一 了
以 有效的 将实测曲 线统一到神经网 络的学习样本曲线模式域中。部分研究者【6 51 ,注愈到这一 问题, 并提出了 不同的 解决力案, , 但都不是最有效的方法。 可见, 量方法构造神经网 点向 络 来识别试井解释模型存在很大的 局限性。针对点向 量方法的不足, 本文提出了基于神经网络 的坐标变换方法和二值化处理方法。 首先, 利用神经网 络建立实测曲 线与样本曲线间的映射 关系, 实现实测曲 线与学习样本的 初步拟合, 然后,由 拟合参 数建立实测曲线的坐标平移关 系, 将实测曲 线平移到对应样本曲线位置或附近, 并利用拟合参数对试井参数进行初步估计。 训练神经网 络的学习样本是压力导数曲 线的二值化矢量, 由于实测曲 线已 经平移到学习样本 位置, 故样本曲线与实测曲 线的二值化方法相同。这一方法有效的 解决了 点向 量方法无法解 决的模式转移问题, 保障了实测曲 线特征向 量在其对应样本曲线特征向 量的同一模式域内。 研究表明, 对具有较大噪声的实测曲 线, 点向 量方法构造的 神经网 络不能正确识别其试井解 释模型类别, 标变 而坐 换及二值化处理方法构造的 神经网络方法仍然能够正确识别其试井解 释模型类别。
文献的回顾
1 点系列翰入向且A N . N 方法 由 于利用实测曲 线训练神经网络的学习样本过于庞大,网络将难以收敛,因此现有的 相关研究 报告中均选用理论曲线作为训练神经网 络的学习 样本。 用反映曲 线特征的 点序列作 为训练神经网 络的 输人向 量, 通过点序列的 归一化处理将实测曲线与理论曲线统一到同 一值 域。 设实测或理论曲 线点系列为(. 它们是压力导数和时间的 ( Y, x r ,) 双对数值, 即:
网 络可建立高 h 线参数Ce及蜂值 差A与曲 DS Z 点坐标( ) ( 之间的映 x , Y 射关系e 实测曲 解 1 [ ] 。由 线的 和 顶点坐标洲Y, 9 利用神 络 经网 便能获得最佳 拟合的 样本曲 线参数 DS 拟合点坐 Ce值及 2 标,
实测曲线的平移处理
对于实侧压力导数曲 线, 用理论曲线点序列作为样本训练出 来的神经网络难以正确识 别其 模式 类别, 其主要原因 是与理论曲 线相比, 实测曲 线存在平移、 不完整以及噪声等, 不 完整和平移的综合影响将使得点向量处理方法难以有效的将同一模式的实测曲线统一到训练 网络的学习样本摸式域中。如果先将实测曲线移动到其对应的样本曲线位皿或附近,即先将 实测曲 线移动到其a练 i 样本的 l 模式域中, 则实测曲线的 早期和后期不完整性将被处理为拾人 矢t的噪声点, 从而可保障神经网 络对实测压力导数曲 线模式类别的 正确识别。 许多数学变换如傅立叶变换等可以实现曲 线的平移且保持曲线特征不变, 但由于实测 曲 线并不是样本曲 线的完整平移,因 而难以找到有效的数学变换将其平移到对应样本曲 线位 置。 不管实测压力导数曲 线相对样本曲 线的噪声及平移量多大, 要其反映试井解释模型类 只 别的特征存在, 则其必然与某类样本曲 线的特征一致。 基于实测曲 线与样本曲线的共同 特征, 利用神经网 络的函 数通近和多物人输出映射特性, 可以建立实测曲 线与其样本曲 线的映射关 系。 从而将实测曲线平移到其对应的 样本曲线位置。 1 导数曲线的共同特征 . 无因次压力导数曲线通常都是在某种参数的组合下绘制的, 如用于图版拟合分析的
皮的 组合值 l 有关, 藏和外 模型 别无关; C e 与油 边界 类
无限大径向流阶段, 所有油藏类型均出 现导数0 线, . 非均质油藏在0 线上出现相应的 5 . 5
波动。
不同试井解释模型的导数曲线均具 有上 述共同特征, 传统图版拟合分析首 先比 较的 也是实测曲 线与样本曲 线的 早 期峰值及0 水平线特征, . 5 通过这两个
.....口 甲叨..日日
龚 ・
. . … , . … . , 闷 二
居 。
2, ,、 ’. 、
奋月
. 月
. 月
, 自
.3 预处理前的曲 续特征
圈4 分段处理后的曲公特征
10 9
窍
一 百 万1
曲线早期和晚期特征都是识别试井解释模型的重要组成部分,它们分别代表组成试井 解释模型三大部分中的 井筒和边界两大部分,因 而这种特征变化必然影响神经网 络对试井解 释模型的正确VM。 I吐
x= ()或 x= ( ) ; g; l t o , g . lt o
y=o [t)或 Y= [ tr ) i ( ' 1 pl g ; r ( v ) l P n o g
则具 有代表性的 归一化方法有以 下几种:
① 疙(2, Y }Y = x x x 、{ , ,) : 一 1 ) ; , 二/ , (2
很差。本文从试井理论分析出发,利用压力导数曲线的固 ) 有特征,由 神经网 络建立实测曲线与 其对应理论曲线间的映射关系,采用坐标平移的方式将实测曲线移动到对应理论曲线位里。利 用二位化方法将试井压力导数曲线表示成二值向t形式,可有效维护理论曲线和实测曲线的原 始模式特征,从而使得所构造的神经网络能够正确识别实测曲 线的棋式类别。研究表明,对于 存在早期和后期不完整的实测试井曲线的试井解释模型,点向 t神经网络方法的正确识别率很 低,而坐标平移与二值化处理的神经网络方法可全部正确识别。此外,在应用坐标平移及二值 化神经网络方法识别试井解释模型的同时,还可以对模型参数进行初步佑计。