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

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

考虑二次梯度项影响的双渗模型的动态特征

在传统试井模型的非线性偏微分方程中根据弱可压缩液体的假设忽略了二次梯度项,已经知道在试井较长时间忽略二次梯度项将产生误差,故对于双重介质流动系统保留了非线性偏微分方程中的二次梯度项,建立了双渗流动模型.采用Douglas-Jones预估校正法获得了无限大地层定产量生产时和定压生产时双渗模型的数值解,分别讨论了液体压缩系数和双重介质参数变化时压力变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

部分射开压裂直井渗流压力动态分析

建立符合非均质油藏部分射开地层实际的物理模型和三维各向异性矩形油藏的不稳定渗流数学模型,考虑不渗透顶、底边界和定压顶、底边界等边界条件的组合,通过无因次量纲变换、拉普拉斯变换、傅里叶余弦变换和分离变量等方法,得到拉普拉斯域的解析解,利用斯蒂芬森数值反演方法,得出实数域的压力数值解.绘制压力动态曲线,并进行敏感性分析.计算结果与数值模拟基本吻合,证实方法的可靠性.敏感性分析表明:压力动态曲线可分为早期线性流、中期径向流、晚期球形流、边界控制流四个流动期.裂缝长度主要影响早期线性流,渗透率各向异性主要影响中期径向流,储层射开程度和裂缝方位主要影响晚期球形流,边界条件和油藏宽度主要影响边界控制流.该方法可以确定最优射开程度、垂向渗透率等参数,为油藏工程分析和压裂工艺设计提供指导.     

混浊介质多光子激发荧光显微成像的蒙特卡罗模拟

将蒙特卡罗方法与几何光线追踪技术相结合用于研究样品混浊特性对多光子激发蒌光显微成像的影响。推导了模拟多层样品显微成像的数学模型,给出了用程序框图表示的基本仿真系统骨架,对传统蒙特卡罗方法的改进有效地增强了对聚焦光束的模拟能力,初步的模拟结果说明了混浊介质散射特性对显微成像的显著影响。     

带有垂直裂缝定压油气藏模型的数值解

压裂作为低渗油气藏开发和油田增产的措施之一,被油田广泛应用.为使压裂方案经济合理,达到少投入多产出的目的,并提高压裂设计水平,我们对油气藏定压生产条件下压裂井产量预测的数学模型进行了数值求解,指出了各种模型的适用条件和要考虑的因素.本文主要介绍区域有裂缝存在、系数还依赖未知函数的非线性抛物型油气藏方程的数值求解。     

Nonlinear dispersion equation and guided modes in a slab waveguide composed of a negative-index medium

We consider a slab waveguide made of a material with negative permeability and negative permittivity, the surrounding medium is conventional. By using the dispersion relation, we show the dispersion curves from the variation of the core width and other physical parameters. The associated energy flow for several thicknesses of the nonlinear layer waveguide is presented. We discuss possible new types of modes that can occur. It is found that the dispersion curves and the energy flow can be controlled by varying the thickness of the guide, the filling factor and the effective plasma frequency. A numerical simulation is done to well illustrate.     

考虑渗透率张量的非均质油藏有限元数值模拟方法

针对具有混合边界的非均质油藏,考虑全张量形式的渗透率,建立弹性微可压缩单相流体不稳定渗流问题的数学模型.根据变分原理,将压力微分方程的边值问题转化为泛函的极值问题,建立渗流模型的有限元方程.针对典型的均质和非均质渗流问题进行模拟计算,得到油藏内压力动态分布曲线,并分析曲线特征.研究表明,有限元法计算精度很高,适用于求解利用渗透率张量表征的非均质油藏渗流问题.为非均质油藏的开发和精细油藏数值模拟提供了理论依据.     

增强型地热系统地下渗流场的模拟分析

采用自行开发的数值模型,模拟分析循环流体在双井增强型地热系统(EGS)地下热储中的渗流过程,系统分析不同水力渗透条件热储中流场形态的变化规律.结果表明,重力作用和流体动静压转换造成的"基础压差"与循环流体在地下裂隙岩体中的沿程阻力之间的相对大小是决定热储中流场分布形态的重要因素,而循环流量大小对热储内渗流场形态影响非常有限.基于此,我们提出水平井和多井环绕布局两种抑制流体短路的方案,为EGS的建设提供理论指导.     

Chirped and chirpfree soliton solutions of generalized nonlinear Schrödinger equation with distributed coefficients

Using the F-expansion method, we systematically present exact solutions of the generalized nonlinear nonlinear Schrödinger equation with varying intermodal dispersion and nonlinear gain or loss. This approach allows us to obtain large variety of solutions in terms of Jacobi-elliptical and Weierstrass-elliptical functions. The chirped and unchirped spatiotemporal soliton solutions and trigonometric-function solutions have been also obtained as limiting cases. The dynamics of these spatiotemporal soliton is discussed in context of optical fiber communication. To visualize the propagation characteristics of chirp and unchirped dark-bright soliton solutions, few numerical simulations are given. It is found that wave profile of solitons depend on the group velocity dispersion and the gain or loss functions.     

Binary Level Set Methods for Dynamic Reservoir Characterization by Operator Splitting Scheme

In this paper, operator splitting scheme for dynamic reservoir characterization by binary level set method is employed. For this problem, the absolute permeability of the two-phase porous medium flow can be simulated by the constrained augmented Lagrangian optimization method with well data and seismic time-lapse data. By transforming the constrained optimization problem in an unconstrained one, the saddle point problem can be solved by Uzawas algorithms with operator splitting scheme, which is based on the essence of binary level set method. Both the simple and complicated numerical examples demonstrate that the given algorithms are stable and efficient and the absolute permeability can be satisfactorily recovered.     

裂变中子角关联分布的数值模拟

为了实现对自发裂变中子源能量、多重性和中子角关联分布的数值模拟,开发了根据裂变过程,逐步给出裂变中子多重性、能量和角分布等参数的NESF数值模拟程序,并将该程序集成到现有中子输运软件中.重点介绍我们开发的用于模拟自发裂变中子源能量、多重性和中子角关联分布的程序的背景、物理依据、程序流程,以及基于该程序开展的中子角关联分布特征与核部件特征关联的数值模拟.     

自由电子激光的计算机模拟

本文使用1(2/2)维电磁相对论粒子程序对自由电子激光进行计算机模拟。在线性段,计算机模拟与用解析理论计算的增长率符合很好。与没有轴向引导磁场的自由电子激光相比,加一弱轴向引导磁场时,Wiggler场的强度可降得较低。加一强轴向引导磁场时,波能量达到饱和水平所需时间较短。     

基于JMCT软件的中子活化数值模拟程序的开发和检验北大核心CSCD

中子活化产物和辐射特征的数值模拟程序是研究材料活化效应的重要工具。在JMCT软件的基础上开发了具备材料中子活化效应模拟能力的数值模拟程序,并将其命名为“中子活化数值模拟程序”,旨在将其应用于军控核查、核安全等领域的研究中。对该程序在核弹头内部中子输运和活化计算的准确性进行了验证,发现该程序对核弹头内部中子输运和活化的计算精度优良。利用该程序研究了混凝土地面核素在裂变核材料的裂变中子辐照下的活化效应,计算结果进一步验证了中子活化数值模拟程序的功能。     

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

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

Computational Model of Thermal‐Fluid Flow in a Radio‐Frequency Plasma Torch

The paper aims to clarify the modelling results concerning the heat transfer and fluid flow in a radio‐frequency plasma torch with argon at atmospheric pressure. Fluid numerical simulation requires the coupling of magnetohydrodynamics (MHD) and thermal phenomena. This model combines Navier–Stokes equations with the Maxwell's equations for compressible fluid and electromagnetic phenomena successively. A numerical formulation based on the finite element method is used. In this study, fluid flow and temperature equations are simultaneously solved (direct method, instead of using the indirect method) using a finite elements method (FEM) for optically thin argon plasmas under the assumptions of local thermodynamic equilibrium (LTE) and laminar flow. Appropriate boundary conditions are given, and nonlinear parameters such as the thermal and electrical conductivity of the gas and input power used in the simulation are detailed. We have found that the source of power is located on the torch wall in this type of inductive discharge. The center can be heated by conduction and convection via electromagnetic phenomena (power loss and Lorentz force). (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

聚变堆液态包层准二维磁流体湍流模型研究

为研究聚变堆氚增殖包层中液态金属湍流磁流体动力学(MHD)效应,开发了一种新的准二维单方程 MHD 湍流模型,并进行了相关数值模拟程序的编制及验证。对于矩形管道中的准二维 MHD 湍流流动,三维流 动主要发生在哈德曼层中,中心的主流区呈现出二维流动。为了反映这种特殊的流动特征,新湍流模型在标准 k-ε 模型的基础之上去掉了传统的耗散项,代之以电磁耗散项来模拟湍流 MHD 效应。同时,采用 Bradshaw 假设来对 湍流涡粘系数进行模化。为验证该湍流模型是否合理,编制了相关数值模拟程序,并利用直接数值模拟(DNS)结 果对该程序进行了校正,数值模拟结果与 DNS 结果吻合较好。计算结果表明,该湍流模型可应用于聚变堆液态 包层 MHD 湍流流动的数值模拟。     

Nonlinear Analysis of Elliptical Waveguide FEM Amplifier

A set of nonlinear three-dimension equations for the FEM amplifier with elliptical waveguide and planar wiggler is derived. Based on these equations, a universal numerical simulation program is developed independently. The characterisitics of this FEM including the evolution of power, saturation efficiency and bandwidth are investigated. Electron momentum spread and wiggler taper are also considered.     

椭圆波导FEM放大器的非线性特性研究

 以椭圆波导、平板摇摆器为FEM放大器的模型，导出了自洽的注波互作用三维非线性方程组。在此基础上编制了相应的程序。利用该程序分析了FEM放大器的饱和功率、效率以及带宽等高频特性。考虑了摇摆器渐变、空间电荷波、电子初始速度零散对高频特性的影响。     

Passive harmonic mode-locking in a fiber laser with nonlinear polarization rotation

On basis of numerical simulation and analytical treatment, we investigate the multiple pulse passive mode-locking of a fiber laser with nonlinear loss due to the nonlinear polarization rotation technique. Various additional mechanisms resulting in the ordering of ultrashort pulses inside the laser resonator and in the realization of harmonic passive mode-locking are analyzed. Among which are the active harmonic modulation of the intracavity loss or of the refractive index, the passive modulation of the index due to the inertial properties of intracavity elements, and the amplification modulation due to the depletion and the relaxation of the gain. The velocity of pulses relative to each other is determined, it results in attractive or repulsive pulse-pulse interaction, and allows to evaluate the conditions of ordering of the pulses inside the cavity.     

An FDTD Analysis of Nonlinear Photonic Crystal Waveguides

We present a modelling technique for analyzing dispersion characteristics of nonlinear photonic crystal waveguides. This technique combines the nonlinear finite-difference time-domain method, based on the numerical simulation of oscillating dipole radiation, with the super cell approach and the periodic boundary conditions. The technique presented makes it possible to calculate the dispersion characteristics of Kerr-like nonlinear photonic crystals waveguides made by removing some scatterers. A numerical tool based on this technique can be used to design various components for integrated optical circuits.