多层油资源运移聚集的数值模拟和实际应用

研究了三维问题的石油地质和渗流力学特征。提出了数学模型,构造了新的算子分裂隐式迭代格式,并成功地对胜利油田东营凹陷(厚层)和惠民凹陷(多层)带断层和通道的实际问题进行了大规的数值模拟计算。计算结果在油藏位置等方面和实际情况相吻合,完整地解决了这一石油地质和渗流力学的国际着名问题。     

三维非线性对流扩散问题的数值方法在渗流力学的应用

对三维非线性对流扩散问题提出一类适合并行计算的二阶迎风分数步差分格式,采用分数步技术,将三维问题化为连续解3个一维问题计算．利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、微分方程先验估计的理论和技巧,得到收敛性的最佳阶的误差估计．该方法已成功的应用油资源运移聚集渗流力学数值模拟计算、海水入侵预测和防治的工程实践中．     

非线性多层渗流系统的数值方法及其应用

对多层非线性渗流系统提出适合并行计算的迎风分数步差分格式，利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧，得到收敛性的最佳阶的误差估计．该方法已成功的应用到油资源渗流力学运移聚集数值模拟的生产实际中．     

渗流耦合系统动边值问题特征分数步差分方法

对多层渗流方程耦合系统动边值问题,提出适合并行计算的一类特征分数步差分格式,利用区域变换、变分形式、粗细网格配套、乘积型叁二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的最佳阶l2误差估计．该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果．     

三维油资源运移聚集的模拟和应用

盆地发育的运移聚集史数值模拟,其功能是重建油气盆地的运移聚集演化史,它对于油资源的勘探,确定油藏位置和计算油藏贮量,寻找新的油田,具有极其重要的价值.本文从地质科学实际出发,研究了三维问题的地质和渗流力学特征,提出二阶修正算子分裂隐式迭代格式,对着名油水运移聚集实验进行了数值模拟,结果基本吻合,并对胜利油田东营凹陷的实际问题进行数值模拟试验,结果和实际地质情况(油田位置等)基本吻合,成功解决了这一着名问题.     

滩海地区运移聚集的精细数值模拟和分析

对滩海地区三层油资源运移聚集进行高精度精细平行数值模拟,提出数学模型和精细平行算子分裂隐式迭代格式,设计了并行计算程序,提出了并行计算的信息传递和交替方向网格剖分方法. 并对不同的CPU组合进行并行计算和分析,对滩海地区数值模拟结果和实际情况吻合.对模型问题进行数值分析,得到最佳阶误差估计,成功地解决了这一困难问题．     

三维渗流耦合系统动边值问题迎风差分方法的理论和应用

对多层渗流方程耦合系统动边值问题,提出适合并行计算的两类迎风差分格式,利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的l~2误差估计.该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果.     

三维非线性多层渗流方程耦合系统的差分方法

对三维多层非线性渗流方程耦合系统提出适合并行计算的二阶迎风分数步差分格式, 利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l²误差估计. 该方法已成功地应用到油资源运移聚集数值模拟的生产实际中.     

非线性渗流耦合系统动边值问题二阶迎风分数步差分方法

对多层非线性渗流方程耦合系统三维动边值问题, 提出适合并行计算的一类二阶迎风分数步差分格式, 利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶l² 误差估计. 该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中, 得到了很好的数值模拟结果.     

多层非线性渗流耦合系统的特征分数步差分方法

对多层非线性渗流耦合系统提出适合并行计算的特征分数步差分格式, 利用变分形式、能量方法、粗细网格配套、分片双二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l₂误差估计. 该方法已成功的应用到多层油资源评估的生产实际中.     

Mixed algorithm in searches of mechanical system steady-state conditions for low precision of the state estimation

This paper focuses on algorithms of numerical solution of nonlinear system of equations. Analytical formulas of their nonlinear functions my not be calculated with the requested precision. Additionally, analytical formulas of the partial derivatives are unknown. They are evaluated numerically by finite differences method. It effects in erroneous estimations. Described situation is critical when steady state conditions are searched for mechanical systems. According to the precision of the numerical procedures, their dynamic equations are known with limited precision. This same stands for the system's final conditions (obtained by a numerical integration). It the actual case, the classical Newton-Raphson algorithm can be ineffective. As an alternative, a mixed search algorithm is proposed in the paper. By contrast to the classical algorithm, within the search direction defined by the Newton-Raphson indication, (potentially erroneous) a local minimum is searched. Both the algorithms are tested on analytical functions; on randomized functions and on a model of a mechanical system. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

流通量间断守恒高阶交通流模型及其数值模拟

在非均匀道路条件下，推广了各向异性守恒高阶交通流模型(CHO模型)，获得流通量间断CHO模型，并基于其Riemann不变量性质，运用局部简化方法及δ映射算法，设计了求解流通量间断CHO模型的一阶Godunov、EO(Engquist-Osher)和LF(Lax-Friedrichs)等数值格式.通过数值模拟表明流通量间断CHO模型是合理有效的，它可以描述平衡态和非平衡态交通流，相对于流通量间断LWR(Lighthill-Whitham-Richards)模型，其能更好地刻画实际交通现象.     

Convergence of Dirichlet quotients and selective decay of 2D magnetohydrodynamic flows

The selective decay phenomena have been observed by physicists for many dynamic flows such as Navier-Stokes flows, barotropic geophysical flows, and magnetohydrodynamic (MHD) flows in either actual physical experiments or numerical simulations. In the previous paper (M.-Q. Zhan, 2010 [20]), the author showed the validity of the selective decay principle for the 2D magnetohydrodynamic (MHD) flows in the case of small magnetic Prandtl number. In this paper, we shall show the validity of the selective decay principle for the 2D magnetohydrodynamic (MHD) flows for any magnetic Prandtl number with periodic boundary conditions.     

Non-axisymmetric consolidation of poroelastic multilayered materials with anisotropic permeability and compressible constituents

This paper presents an analytical layer-element solution to non-axisymmetric consolidation of multilayered poroelastic materials with anisotropic permeability and compressible constituents. By applying Fourier expansions, Hankel transforms and Laplace transforms to the state variables involved in the governing equations of poroelasticity with respect to the circumferential, radial and time coordinates, respectively, the analytical layer-element (i.e. a symmetric stiffness matrix) is derived, which describes the relationship between the transformed generalized stresses and displacements at the surface (z = 0) and those at an arbitrary depth z, considering the corresponding boundary and continuity conditions at the layer interfaces, the global stiffness matrix of a multilayered system is assembled in the transformed domain. The actual solutions in the physical domain are acquired by applying numerical quadrature schemes for the inversion of the Laplace–Hankel transform. Finally, numerical calculation is presented to investigate the influence of layering and poroelastic material parameters on consolidation process.     

Semi infinite Orthotropic Cantilevered Strips and the Foundations of Plate Theories

Elastostatic problems of semiinfinite orthotropic cantilevered strips with traction-free edges and loading at infinity are reduced to the solution of a single scalar Fredholm integral equation of the first kind with a generalized Cauchy kernel. The known complex variable method for equations with a Cauchy type kernel is extended to handle the singularities in the solution for the generalized Cauchy kernel. The reduced problem lends itself to a more efficient numerical solution scheme than all existing methods. Moments of stresses at the root of the cantilever are accurately evaluated and used for the correct formulation of displacement boundary conditions for a plate theory solution (or the actual interior solution) of the elastostatics of thin flat bodies.     

采用高分辨率数值方法重构与实测数据一致的实际爆炸波的部分流场

本文采用球面爆炸的活塞推进模型和对间断解具有高分辨率的二阶Godunov型有限差分方法,以两类化爆试验值作依据,通过恰当调整运动活塞边界条件,数值重构了实际半球形爆炸波的部分流场.方法简单、可靠,所得结果可用于评估任意同类型炸药一维爆炸波的远场效应.     

Equilibrium attractivity of Krylov-W-methods for nonlinear stiff ODEs

Since implicit integration schemes for differential equations which use Krylov methods for the approximate solution of linear systems depend nonlinearly on the actual solution a classical stability analysis is difficult to perform. A different, weaker property of autonomous dissipative systemsy′=f(y) is that the norm ‖f(y(t))‖ decreases for any solutiony(t). This property can also be analysed for W-methods using a Krylov-Arnoldi approximation. We discuss different additional assumptions onf and conditions on the Arnoldi process that imply this kind of attractivity to equilibrium points for the numerical solution. One assumption is general enough to cover quasilinear parabolic problems. This work was supported by Deutsche Forschungsgemeinschaft.