非饱和流问题的数值模拟研究

利用有限元方法 ,建立了非饱和土壤水流动的数值计算模型 ,用以模拟均质土壤 ,地下水埋藏很深且在不同的初始和边界条件下的水分运动 .对于入渗或蒸发问题 ,用通常的离散方法容易产生振荡非物理解 .我们采用有限元集中质量法 ,将质量矩阵的非对角元集中到主对角元 ,避免了这种非物理解的产生 ,并且较好地处理了这类边界条件 .     

基于阵列探针组合仪的垂直井油水两相数据分析与模型建立

传统垂直井油水两相产出剖面采用单一持率仪和单一涡轮以居中的方式进行测量,只采集井筒中心处的局部持率和局部流速信息,只反映井筒中心处的流体信息,不能直观地反映井筒横截面上的流体流动状态,井筒中含水率从低含水到高含水的过渡时存在测量盲区．为此,本文首次尝试采用阵列探针组合仪（包括：阵列电阻持率仪、阵列电容持率仪和阵列涡轮仪）,与常规单一电容持率仪、全井眼涡轮仪进行组合的方法,在垂直井油水两相中开展模拟实验．首先,通过实验数据和理论流型,绘制模拟实验流型图,并对流型流态进行二维成像．其次,利用三种持率仪进行了局部持率和平均持率的计算,对两种涡轮仪进行了局部流速的计算,并利用体积法计算出平均流速．然后,针对不同的仪器组合分别建立了漂移模型,进行了油水流量的计算及对比分析．研究表明：垂直井油水两相实验中,采用阵列探针组合仪测得的数据具有较高的精度,利用它们建立的漂移模型,能更精确的计算出油相和水相的流量,对垂直井油水两相中的动态监测提供一定的参考．     

水平井生产测井组合仪模拟井测量数据分析与模型建立

由于水平油水两相流局部流速及局部相含率沿管径方向存在复杂分布,致使其流量及相含率测量非常困难.本文利用伞集流涡轮流量计、电导传感器、过流式电容传感器构成的生产测井组合仪,在水平油水两相流模拟井中开展了动态实验,分析了流动特性及管子角度倾斜对涡轮流量计及含水率传感器响应特性的影响.根据电容及电导传感器在不同总流量及含水率范围的响应分辨特性,分别采用变系数漂移模型及统计模型预测分相流量,发现变系数漂移模型能够自适应调整模型中的流型参数,且对油水两相分相流量具有较高预测精度.     

自由场地液化振动台实验结果分析与数值模拟

通过拉格朗日差分程序FLAC3D完全非线性动力分析模块,模拟了自由场地液化振动台实验并对实验结果做出分析,结果表明:该动力分析方法能够较好地反映液化土层地震响应,实验结果也验证了模拟结果的精确性;液化土层滤波作用起始于孔压比上升的时候且较为明显,液化后加速度时程的卓越周期有所增加;本文实例中的液化土层对地震动具有放大作用,液化范围越大,其放大作用越明显,且在长周期段放大作用更为突出。     

5°拐折构造变形物理场的实验研究与数值模拟

以实验方法和数值模拟方法研究了含５°拐折的断层在变形过程中物理场的时空演化，结果表明在一次事件序列中拐折部位的声发射往往发生在序列开始时和结束时，震级虽不大，但却有重要意义，在断层错动中拐折部位两侧的断层交替地作为变形主体构造，拐折部位则起了错动阀的作用。     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

松散承压含水层水位的同震响应实验与数值模拟

为深入理解井水位同震响应机理,本文开展了向完整井-松散含水层系统输入由不同频率和振幅（加速度）组成的正弦波荷载的振动台实验。以实验模型为物理模型,建立了振动作用下松散承压含水层中孔隙水压力响应的流固耦合模型和含水层水流与井流的相互作用模型,并运用多物理场耦合模拟软件COMSOL Multiphysics对实验过程进行了数值模拟。实验中观测到的四种典型水位变化形态与野外场地同震井水位变化形态相似。数值模拟结果显示,本研究建立的数学模型能较好地反映松散承压含水层中孔隙水压力和水位的响应情况。本文研究对解释地下水同震响应机制、岩体渗流稳定性和安全问题具有重要意义。     

不同掘进工艺煤与瓦斯流固耦合数值模拟研究

井下瓦斯事故严重威胁着煤矿的安全,研究煤与瓦斯耦合作用规律并采取相应措施可以有效地防止事故发生。考虑瓦斯气体的可压缩性和吸附、解吸特性,假设瓦斯气体在煤体孔隙中的流动遵循D arcy定律,建立了考虑爆破因素的煤与瓦斯流固耦合数学模型,并利用com so l对其求解。结果表明：炮掘对煤体的破坏能力大于机掘;炮掘工艺下,爆破应力对煤体应力、位移的作用大于应力耦合的作用,机掘工艺应该考虑应力耦合的作用;煤层瓦斯压力随着煤壁暴露时间的增加而发生改变,整个巷道附近的瓦斯压力有降低的趋势,呈现漏斗状分布;煤层上表面煤体竖向位移边界呈以巷道顶端为中心的＂O＂型分布,且以＂O＂型中心处的位移值最大。     

基于多相流模型和探地雷达正演模拟的LNAPLs探测研究

探地雷达(Ground Penetrating Radar:GPR)不仅可以用于估计土壤含水量,还可以用于探测和监测轻非水相液体(Light Non Aqueous Phase Liquids:LNAPLs)在土壤中的运移.建立接近实际情况的模型是利用GPR正演模拟开展LANPLs迁移和分布研究的关键问题.以往的地球物理模型大都属于概念模型,存在物性突变的界面.而在LNAPLs污染区域的各种物性参数通常是渐变的,大多数情况下并不存在突变的物性界面.因此,为了建立更符合实际的地球物理模型,本文基于多相流渗流理论模拟了LNAPLs在土壤中的泄漏过程,得到了渗漏后不同时刻土壤含水饱和度、含油饱和度的变化分布.然后利用混合介质的介电模型将流体饱和度分布转换为介电常数分布,获得了地球物理模型,显示了因LNAPLs迁移引起的土壤介电常数的细节变化.随后,基于时域有限差分开展了GPR正演模拟.正演模拟结果显示了雷达波对LNAPLs污染区域、潜水面的响应,与实验室实测数据具有很好的一致性.由以上分析、对比可知,本文提出的地球物理建模的方法和流程与污染场地的实际情况更为符合.基于多相流渗流理论建立的地球物理正演模型把地下LNAPLs迁移的水文模型与GPR探测相结合,为复杂的实际场地地球物理建模提供了思路,也为GPR更有效地探测LNAPLs在土壤中的渗流提供了分析和解释手段.     

大吨位双出杆粘滞阻尼器开发的数值模拟与实验研究

阐述了大吨位双出杆粘滞阻尼器开发的数值与实验研究结果.基于硅油的幂律特性,通过对硅油介质的流态分析,介绍了非牛顿幂律流体模型对粘滞阻尼器耗能和阻尼出力的影响,与粘滞阻尼器的分析模型和计算公式.在商用流体软件平台上利用动网格技术对样品粘滞阻尼器阻尼出力进行了计算;综合实验研究结果,最终确定了大吨位土木用粘滞阻尼器设计的经验方法.另一方面,实验曲线和数值曲线基本吻和,验证了通过数值手段对粘滞阻尼器的出力特性进行预分析,从而提高开发精度、指导产品设计和缩短产品的开发周期的可行性.     

Electromagnetic holographic sensitivity field of two-phase flow in horizontal wells

Electromagnetic holographic data are characterized by two modes, suggesting that image reconstruction requires a dual-mode sensitivity field as well. We analyze an electromagnetic holographic field based on tomography theory and Radon inverse transform to derive the expression of the electromagnetic holographic sensitivity field (EMHSF). Then, we apply the EMHSF calculated by using finite-element methods to flow simulations and holographic imaging. The results suggest that the EMHSF based on the partial derivative of radius of the complex electric potential φ is closely linked to the Radon inverse transform and encompasses the sensitivities of the amplitude and phase data. The flow images obtained with inversion using EMHSF better agree with the actual flow patterns. The EMHSF overcomes the limitations of traditional single-mode sensitivity fields.     

生产油井井下温度场数值模拟分析

An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature field model in cylindrical coordinates is developed. In the model, we considered general heat conduction as well as the heat convection due to fluid flow from porous formation to the borehole. We also take into account the fluid velocity variation in the wellbore due to multiple pay zones. We present coupled boundary conditions at the interfaces between the wellbore and adjacent formation, the wellbore and pay zone, and the pay zone and adjacent formation. Finally, an alternating direction implicit difference method （ADI） is used to solve the temperature model for the downhole temperature distribution. The comparison of modeled temperature curve with actual temperature log indicates that simulation result is in general quite similar to the actual temperature log. We found that the total production rate, production time, porosity, thickness of pay zones, and geothermal gradient, all have effects on the downhole temperature distribution.     

大斜度井和水平井井地三维电阻率数值模拟和联合反演

井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向.

     

大斜度井和水平井井地三维电阻率数值模拟和联合反演

井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向.     

Study on the simulation of acoustic logging measurements in horizontal and deviated wells

The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.     

于桥水库风季流场的计算与分析

本文依据水流运动的基本理论,并考虑风对水体的剪切作用,应用改进的移步ADI法,建立了风作用下流场的数值模似系统,其结果与连续四年的实测资料验证相吻合。计算预报了于桥水库三个典型水文年定常风不同风速、风向的流场,并由计算机绘制出彩色流场图。文中还针对计算应用成果进行了分析。     

一个分层水库温跃层的模拟与验证

以西安金盆水库为例,建立了分层水库水温结构的数值模拟方法,并以实测数据进行模型验证.运用Fluent软件数值研究了不同短波辐射强度及短波辐射衰减系数条件下温跃层的形成过程与特性.水库水面总传热量在春、夏季为正值,在秋、冬季为负值,长波辐射是水面总传热量的主要影响因素,短波辐射则是温跃层形成的主要影响因素.随短波辐射衰减系数的降低,温跃层厚度增加,温跃层内温度梯度减小,短波辐射衰减系数值与实测的藻类浓度存在良好的正相关性.水库具有极限短波辐射强度,温跃层内温差随水面短波辐射强度的增加呈现先增加后减小的变化趋势;但水面短波辐射强度过高时,难以达到热平衡而形成稳定的温跃层.     

青藏块体东北缘地壳水平应变场数值模拟研究

基于青藏块体东北缘1999～2001年GPS结果,分别采用块体整体旋转与线性应变模型和弹性力学有限元法这两类地壳形变数值模拟方法,分析了该区地壳水平应变场特征.结果表明:(1)两类方法在研究地壳形变时各具一定的优势,前者对块体整体运动变形及块体与块体之间的相关性研究具有一定的优势,而后者则较强体现出了应变高值区与深大断层在空间分布的紧密结合性;(2)两类方法所获得的应变高值区具有良好的空间分布一致性,主要集中在阿尔金断裂、祁连山断裂中东段、东昆仑断裂和海原断裂这些深大断裂处及其附近,面压缩值达到了-3×10-8以上,最大剪应变值达到了10×10-8以上;(3)应变高值区的空间分布与中强地震的发生具有一定的对应关系.     

Ultrasonic method for measuring water holdup of low velocity and high-water-cut oil-water two-phase flow

Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China’s onshore oil fields are mostly in the high-water-cut production stage. This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration (of oil droplets) in oil-water two-phase flow, which makes it difficult to measure water holdup in oil wells. In this study, we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in low-velocity and high water-cut conditions. First, we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling. Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor. Based on the results, we then investigate the effects of oil-droplet diameter and distribution on the ultrasonic field. To further understand the measurement characteristics of the ultrasonic sensor, we perform a flow loop test on vertical upward oil-water two-phase flow and measure the responses of the optimized ultrasonic sensor. The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow (D OS/W flow), but the resolution is favorable for dispersed oil in water flow (D O/W flow) and very fine dispersed oil in water flow (VFD O/W flow). This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.     

海洋拖曳式水平电偶源数值模拟与电场接收机研制

海洋天然气水合物资源调查中,电磁法已成为地震勘探最有效的补充手段之一,而我国海洋水合物电磁探测中的拖曳式电场接收机此前是空白.本文从方法原理、仪器设备和海洋试验等方面,对海洋拖曳式水平电偶极-偶极方法进行研讨.建立三维电阻率模型,进行数值模拟,讨论收发距、发射频率、拖曳装置离海底高度等参数与电场异常响应的关系,针对特定模型给出了最佳观测参数以指导仪器硬件设计与海洋测试.在此基础上,采用高级嵌入式计算机的智能化控制、前端低噪声斩波放大、高精度同步采集等技术,研制了由拖体、承压舱、采集电路、低噪声电场传感器等组成的海洋拖曳式水平轴向电场接收机,其系统噪声达到2 nV/√HZ@(1~100Hz).将自主研制的电场接收机用于国内首次拖曳式电偶极-偶极方法的海洋试验与评估,并对海试资料进行了处理与分析.研究结果表明,所研发的仪器达到了设计指标,能够用于海洋天然气水合物调查.