混合井网条件下加密井产能计算模型及应用

对于以直井、水平井混合井网开采的油藏,后期加密的水平生产井产能计算应综合考虑其他油水井的干扰.基于镜像反映原理与微元线汇理论,建立任意混合井网条件下加密水平井稳态产能计算半解析模型,考虑混合井网内部其他油水井干扰及井筒内由于井壁摩擦与流体加速度造成的压力降对加密水平井产能的影响,适用于任意混合井网条件下(井型任意性、井数目任意性、井位任意性)加密水平井产能计算.     

水平井体积压裂缝网表征及流动耦合模型

针对岩石脆性系数高且发育天然裂缝的储层,提出表征水平井体积压裂形成裂缝网络的三种基本模式,并将渗流过程划分为油藏流动和缝网内部流动.在此基础上,利用势叠加原理导出油藏流动控制方程,利用有限差分方法建立缝网内部有限导流等式;其次,采用星三角变换法处理人工缝与天然缝的交汇流动;最后,耦合两部分流动矩阵方程得到水平井体积压裂缝网渗流数学模型.该模型表明:当水平井改造段长度一定时,压裂段数与段内分簇数是决定产能的最主要因素,其次是人工裂缝半长和人工缝导流能力,而天然裂缝密度和导流能力对产量影响较小.实例应用表明,实际产油量与模型计算值一致,误差较小.     

使用混合网格计算非达西渗流

针对垂直裂缝井的特殊流动模式,从非达西定律出发,建立二维平面的非达西渗流方程.通过建立一组无量纲量,最终得到无量纲的渗流方程及其定解条件.假定外边界为圆形,用PEBI网格及混合网格对求解区域进行网格划分,用有限差分法对无量纲的方程进行离散,最终得到垂直裂缝井的井底压力数值解.根据此数值解并考虑井筒存储和表皮因子的影响,得到真实垂直裂缝井的井底压力.对计算结果的分析表明,使用混合网格求解非达西渗流井底压力相当准确,该方法也适用于水平井等更复杂井型及复杂边界的问题求解.     

低渗透油藏压裂水平井井筒与油藏耦合的非稳态模型

推导低渗透各向异性油藏压裂水平井井筒与盒式油藏耦合的非稳态模型,并给出求解方法.模型考虑摩阻和加速度的影响,并可以使用不同类型的约束条件.实例计算表明,压裂水平井的流动可分为非稳态阶段和拟稳态阶段.在非稳态阶段,各条裂缝的产量相差不大,总产量随着裂缝条数的增加呈线性增加;在拟稳态阶段,两端裂缝产量高于中部裂缝的产量.受摩阻和加速度压降的影响,空间上位置对称的裂缝在流量上呈现不对称性.井筒内压力损失的存在将使水平井的产量降低,并使井筒内的压力分布不均匀.在最小井底流压的基础上固定流量时,裂缝条数越多,稳产期也越长.     

裂缝性应力敏感油藏流固耦合的数值模拟

以裂缝各向异性全张量渗透率动态耦合渗流模型为基础,建立三维三相裂缝各向异性特低渗透油藏流固耦合数学模型,并运用模块化显示耦合迭代求解建立流固耦合渗流数值模拟方法,编制耦合求解程序.研究表明:对于裂缝性油藏,随着油藏开采过程中压力的变化,岩石变形和储层介质的各向异性应变特性使开发过程中多组裂缝的变形程度不同,引起裂缝渗透率的非线性变化,进而使油藏裂缝渗透率张量主值发生旋转,改变液流方向、影响油藏开发效果.     

分形低渗透缝洞型油藏非线性渗流规律

引入分形理论,建立考虑低速非达西效应的分形三重介质缝洞型油藏数学模型,通过Laplace变换及Stehfest数值反演方法求出井底压力,借助Matlab编程绘制压力动态曲线,划分渗流阶段,分析渗流规律,进行非线性参数敏感性分析.最后结合实际算例,验证模型的正确性.结果表明：分形三重介质油藏渗流过程分为早期纯井储,过渡流,缝洞窜流,拟径向流,基质与溶洞、裂缝窜流及总体径向流6个渗流阶段;低速非达西效应对渗流的影响随时间的推移逐渐增大;启动压力梯度越大,总径向流阶段压力动态曲线上翘幅度越大;分形系数影响整个渗流过程,随着分形系数的增大,裂缝迂曲程度随之增大,致使渗流阻力增加,引起压力动态曲线整体上移幅度增大.     

分形裂缝性页岩气藏多段压裂水平井瞬态压力特征

基于分形理论和连续性假设,考虑页岩气吸附解吸、基质-裂缝窜流等机制,建立分形裂缝性页岩气藏多段压裂水平井试井解释模型,并通过拉氏变换、点源函数及压降叠加原理等方法得到模型的解.绘制无因次压力随时间变化的双对数曲线,研究分形裂缝性页岩气藏多段压裂水平井的压力特征,分析分形指数、分形维数等参数对压力动态的影响.结果表明：分形裂缝性页岩气藏多段压裂水平井的压力动态可划分为7个流动阶段;分形指数越大或分形维数越小,晚期径向流直线段的斜率越大;其它参数对水平井的压力动态也有一定的影响.     

考虑水力压裂缝和天然裂缝动态闭合的三维离散缝网数值模拟

基于嵌入式离散裂缝模型, 提出一种可在三维空间中考虑应力状态影响的裂缝动态闭合表征方法。将任意方向裂缝的开度和渗透率考虑为作用在裂缝平面法向有效应力的函数, 同时用裂缝传导率变化表征支撑剂充填的水力压裂缝与被开启的天然裂缝由于油藏开发过程中地层流体压力下降而发生的动态闭合行为。研究表明: 致密油藏开发以缝控储量为主。对压裂水平井进行产能评价时, 裂缝动态闭合会导致产能的部分损失, 其影响不可忽略; 水力压裂缝的支撑剂材料属性及天然裂缝的刚度是其中的主控因素。因此需要增大支撑剂的浓度、粒径大小并改善支撑剂的性质, 在最大程度上降低裂缝闭合对生产的不利影响。     

有限导流压裂定向井耦合流动模型

考虑人工裂缝、斜井筒和地层中的耦合流动,将压裂定向井的流动分为地层向裂缝渗流、裂缝内流动和斜井筒变质量管流.采用镜像反映和势叠加原理,建立地层中势分布和渗流数学模型.应用边值理论,将人工裂缝面离散为裂缝网格,并与斜井筒变质量管流模型进行耦合,建立有限导流压裂定向井耦合流动数学模型,并形成相应的迭代算法.应用分析表明:该模型与Prats图版法相比,对于垂直缝吻合度较高;裂缝导流能力和井斜角对产量及裂缝面压力分布影响较大.     

各向异性渗透率油藏数值模拟

建立三维三相各向异性油藏渗流数值模拟模型并应用于油藏开发.模型包含完全各向异性渗透率张量,采用十九点有限差分格式和全隐式迭代方法求解.可以描述不同区域不同主方向的各向异性渗透率分布,准确计算各向异性油藏内流体运动过程、剩余油分布和压力分布.     

泡沫混排解堵数学模型及数值模拟

建立泡沫井筒流动模型与泡沫地层渗流模型,得到完整的泡沫混排解堵数学模型;利用数值方法对模型进行耦合求解,得到泡沫压力,质量,密度在井筒中的分布及井口井底压力的变化规律.在固定井口回压情况下,得到井底压差的变化规律.结果表明:随着井深的增加,泡沫压力增加,质量减小,密度增加;固定井口回压,井底压差则逐渐变小.     

注蒸汽井筒沿程热物性参数及热损失新算法

基于干度的定义,对传统蒸汽干度方程进行改进,计算井筒热损失速率时,既按深度分段又按注汽时间分段,并采用Hasan等提出的两相流漂移流动模型计算井筒压力降.结果表明:改进后的模型用新算法的计算结果与现场实测值吻合很好,比传统的Beggs-Brill算法更准确.另外,在相同井深处,随着注汽时间的延长,蒸汽压力、温度和干度几乎不变,但是热损失速率逐渐降低.新算法为准确计算井筒热损失提供了参考.     

纵波速度径向层析成像分析

钻井会造成近井壁附近岩石的声学性质改变,为研究纵波径向探测特性,本文使用径向柱状分层模型,应用二维轴对称有限差分法详细考察了过渡带厚度、声源频率以及接收器源距变化下,单极子测井响应的变化特征,并采用慢度-时间相关法研究了不同条件下所测量到的声速变化。最后基于射线理论,分析了不同条件下的首波到时变化特征。基于地层速度逐渐增加的模型表明,纵波的径向探测特性受源距影响较大,通过不同源距信号的组合,可以实现对近井壁纵波速度进行层析成像。     

Resolution of objects within subwavelength range by using the near field of a dipole

We analyze the far-field resolution of apertures that are illuminated by a point dipole located at subwavelength distances. It is well known that radiation emitted by a localized source can be considered a combination of traveling and evanescent waves, when represented by the angular spectrum method. The evanescent wave part of the source can be converted to propagating waves by diffraction at the aperture; thereby it contributes to the far-field detection. Therefore one can expect an increase in the resolution of objects. We present explicit calculations showing that the resolution at the far zone is improved by decreasing the source-aperture distance. We also utilize the resolution enhancement by the near field of a dipole to resolve two closely located apertures. The results show that without the near field (evanescent field) the apertures are not resolved, whereas with the near field of the dipole the far zone intensity distribution shows improved resolution. This method eliminates the requirements of near-field techniques such as controlling and scanning closely located tip detectors.     

Numerical investigation of a coupled moving boundary model of radial flow in low-permeable stress-sensitive reservoir with threshold pressure gradient

The threshold pressure gradient and formation stress-sensitive effect as the two prominent physical phenomena in the development of a low-permeable reservoir are both considered here for building a new coupled moving boundary model of radial flow in porous medium. Moreover, the wellbore storage and skin effect are both incorporated into the inner boundary conditions in the model. It is known that the new coupled moving boundary model has strong nonlinearity. A coordinate transformation based fully implicit finite difference method is adopted to obtain its numerical solutions. The involved coordinate transformation can equivalently transform the dynamic flow region for the moving boundary model into a fixed region as a unit circle, which is very convenient for the model computation by the finite difference method on fixed spatial grids. By comparing the numerical solution obtained from other different numerical method in the existing literature, its validity can be verified. Eventually, the effects of permeability modulus, threshold pressure gradient, wellbore storage coefficient, and skin factor on the transient wellbore pressure, the derivative, and the formation pressure distribution are analyzed respectively.     

Complex ray method and its applications to underwater acoustics

In this paper,the principle of complex source point and various analyt-ical methods of complex ray are summarized.Some possible applications ofcomplex ray method to underwater acoustics are proposed,including complexray analysis of directional acoustic fields,complex ray simulation of directivityof acoustic radiator arrays,complex ray transformation of acoustic fields fromnear zone to far zone,as wel1 as theoretical prediction of sonar target strengthby complex ray method.     

敏感换道对下匝道系统交通流的影响

在NaSh模型基础上,充分考虑驾驶员在匝口指示牌的诱导作用下驾驶方式的变化,定义了车辆在匝口上游的换道、直行驾驶规则,提出了敏感换道的元胞自动机下匝道交通流模型. 通过计算机数值模拟,结果表明：敏感换道过程能减少直行车道上的转出车辆比例,对非必要的换道行为有明显的抑制作用,且随敏感换道区长度增加,该作用越明显;匝口提示位置并非越长越好,系统转出车辆比例越小,系统所需的最佳敏感换道区长度越短. 工程设计中根据转出车辆比例选取一个适宜距离安放匝口指示牌,能有效增加系统流量和临界加入概率. 关键词： 交通流 元胞自动机 驾驶方式 计算机数值模拟     

各向异性油藏水平井多角度人工裂缝线性单元计算方法

考虑各向异性油藏渗透率张量的表征,利用Green函数和拉普拉斯变换建立裂缝流动的一维单元,裂缝的流量分布采用节点线性插值,裂缝内的流动处理为线性积分,耦合地层与人工裂缝的流动,建立有限导流裂缝井底压力的求解方法.结果表明:多裂缝压裂水平井存在压裂裂缝线性流、地层线性流、系统径向流3种流动形态,压裂裂缝条数越多,相同的生产时间,无因次井底压降越小;裂缝条数对流动影响明显.随着裂缝条数的增加,压降变化减小;裂缝长度和导流能力有相似的变化.人工裂缝与井筒角度越大,产能越大,当裂缝垂直于井筒时,产量最大;地层最大渗透率方向垂直于人工裂缝时产量最大,平行裂缝时产量最小.当人工裂缝垂直于井筒,并同时垂直于地层最大渗透率方向时,达到最大产量值.     

基于比例积分控制原理预测钻井全过程原始地层温度的新方法研究

应用比例积分控制原理将瞬态传热模型预测结果与出口 温度实测数据逐步进行反馈可准确预测原始静态地层温度. 为此, 本文基于井下各控制组件质量、动量及能量守恒原理, 建立了实际井身结构与钻具组合条件下循环和停止循环期间井筒-地层温度分布全瞬态传热模型, 应用全隐式有限差分法进行求解, 并引入比例积分控制原理对比分析实测温度与预测温度的误差范围进而精确、 快速获取原始地层温度. 结合一口深井基础数据计算表明, 套管下入深度改变了井筒-地层间热交换效率, 进而影响了近井壁地层温度分布状况; 同时, 钻井过程中循环和停止循环作业过程改变了井下各控制组件的初始条件与边界条件, 致使近井壁原始地层温度分布距离产生变化. 建立的数学模型和研究方法可为石油钻井、地热井开采及地球深部原始地层温度信息准确、 经济、快速获取提供理论基础. 关键词： 原始地层温度 循环与停止循环 瞬态传热模型 比例积分控制原理