三维MANDELBROT过程的渗流下临界状态

本文研究了三维Ｍａｎｄｅｌｂｒｏｔ过程渗流下临界形态，证明了渗流函数θ（Ｐ）在Ｐ＿ｃ点是不连续的，同时得到了当Ｐ＜Ｐ＿ｅ时，三维Ｍａｎｄｅｌｂｒｏｔ渗流Ａ的最大连通分支是点，即Ａ是“尘状”的。     

带吸收项的渗流方程的柯西牛曼问题

本文在半无界域中考察带吸收项的渗流方程的柯西牛曼问题．主要结果为：（１）证明该问题弱解的存在唯一性，（２）证明边值条件不影响初始扰动传播的有限性．     

三维单域弹性问题虚边界无网格伽辽金法分析

针对更具一般性的三维问题,虚边界无网格伽辽金法被进一步推广研究,提出了三维单域弹性问题虚边界无网格伽辽金法,包括ANSYS有限元软件提取面单元、节点数据信息,给出的命令流具有动态数组的优点,输出的节点坐标达到28bit。详细推导了三维单域弹性问题虚边界无网格伽辽金积分方程,并为了便于编程实现进行数值离散,得到积分方程对应的离散格式。最后计算三维混凝土立方体受压试块应力分析,取中间四个截面上的应力进行验证。结果表明,三维单域弹性问题虚边界无网格伽辽金法计算可行、精确性好;给出的ANSYS命令流,能够提前准备编程数据,通用性强,大大简化任意面的单元划分与节点坐标信息提取工作,利于建议方法的推广应用。     

双墩柱波浪渗流力

计算了在线性波或椭圆余弦波作用下双圆墩柱波浪渗流力，讨论了双墩 柱沿波向串列和并列布置时墩柱上渗流力随墩距的变化．计算中在内域采用 有限元方法，在外域使用了本文给出的外场解．进行了双墩柱波浪渗流压力 测量的模型实验．渗流压力的计算结果与实验结果符合较好．     

无限导流水平井动态压力分布研究

针对水平井渗流中一系列问题,以单井作为研究对象,从渗流力学基本理论出发,围绕物理数学模型,综合渗流力学、油藏工程、现代数学等多学科交叉的方法,研究了无限导流水平井的非稳态渗流理论。在此基础上推导了顶底封闭油藏中无限导流水平井的线源渗流模型及体源渗流模型。得到其在时间域、空间域上的动态压力分布。     

真空预压法中渗流体力和压差面力联合作用机理浅析

从渗流问题的基本方程出发，导出了真空预压的固结解析解，认为在地基中存在负压边界引起的瞬时稳定渗流与固结渗流2种不同特性的渗流、瞬时稳定渗流体力使膜下地基竖向拉伸、侧向受压，接近纯剪切应力状态，体积压缩变形较小，难以有效地加固地基．柔性膜上压差面力产生的偏应力与渗流体力的偏应力方向相反，大致可以互相抵消，使地基处于接近球应力状态．因此，膜上压差是加固地基的主要荷载，对地基的固结特性和最终加固效果起着控制作用．     

基于QSGS法3D重构土体渗流场的LBM数值模拟

基于四参数随机生长法（QSGS）以及改进QSGS方法等效重构三维（3D）各向同性/各向异性土体内部孔隙结构，结合格子玻尔兹曼方法（LBM），采用MATLAB自编程序进行细观渗流场数值模拟，分析不同孔隙参数及模型各向异性等因素对渗流特性的影响规律. 结果表明：3D模型尺寸100×100×100格点时，连通孔隙率nc达到最大增幅19.23%. 流体易在连通性好、孔喉直径大的部位形成主渗流通道，且存在“指进效应”，流体越靠近孔道中轴线流速越快. 3D重构土体的计算渗透率k随模型孔隙率n增大而增大，随土颗粒尺寸减小而减小，当分布概率pc = 0~0.05（大、中颗粒）时渗透率显著降低，降低幅度达42.86%. 考虑各向异性的重构模型出口边界渗流速度波动幅度较大、变化趋势规律性不明显，且流线分布稀疏、相互交错、主渗流通道不明显. 该研究能揭示流体在3D土体孔隙结构中的流动规律，并为3D土体重构和细观渗流模拟研究方法提供借鉴.     

渗流分析的非正交网格有限体积法

将有限体积法的体积积分用Ｇｒｅｅｎ公式转化为单元边界积分，使之能应用于非正交网格的渗流分析计算；保留有限体积法计算简单实用、计算精度高、省机时省内存的特点，又有网格划分灵活方便，对复杂计算域适应性强的优点。     

组合分形油藏非达西低速渗流问题

将分形几何结合渗流力学建立了一类组合分形油世故非达西低速渗流的数学模型，即在两个同心圆柱莴心球地层中，中心一口生产井生产时内域为分形油散发非达低速渗流上域为分形油藏的达西渗流。讨论了内边界定流量和内边界定压两种情况，并利用拉普拉斯变换求出了两种情况在两个区域内压力分布的精确解。     

二维渗流的两个结果

本文给出了二维渗流的两个结果：（１）在临界概率时，对包含原点的开串在边长为２ｎ的正方形中点的平均个数给出了一个估计；（２）对平面上无穷子区域，讨论了渗流存在及临界概率为１／２的条件．     

致密油藏有限元数值模拟

 目前数值模拟方法均基于达西渗流,无法对致密油藏非达西渗流进行有效模拟。致密油藏内在复杂,为实现致密油藏非达西渗流的科学模拟,建立了油水2相微可压缩非达西流的油藏数值模拟数学模型。同时考虑到致密油藏渗透率低、非均质性强的特点并兼顾模拟精度采用有限单元法进行离散,从而建立了有限元数值模型,形成了新的致密油藏数值模拟方法,最后将该方法进行实例检验,验证了方法的正确性。     

基于Schwarz Christoffel变换的曲流河地质建模方法

为了解决曲流河物性参数沿着河道方向多变问题。通过封闭多边形区域到矩形区域的Schwarz Christoffel变换基本原理，建立不规则带状封闭曲流河边界映射到规则矩形区域边界的数学模型。采用高斯雅克比积分方法和牛顿迭代法，完成了从曲流河边界到矩形区域边界映射数学模型的数值计算。根据映射边界参数，采用二维粒子群方法完成了平面井位映射数值计算。在地质模型的建立过程中，考虑了映射规则和缩放比例问题。最后以X砂岩油藏为例，选择了该区域中38口井，建立的映射前后的地质模型，并对物性参数进行了分析。通过分析对比，采用Schwarz Christoffel映射后建立的地质模型及分析的物性参数更加符合地质规律。     

岩石面糙率的分布及其对裂隙水流的影响

就岩石央糙率的分布规律及其对裂隙水流的影响进行了计算分析,通过天然大理石板人工构成开裂隙模型所进行的大应力高水头耦合试验研究,对大理石板在糙度的实测数据进行了数值分析,分析表明,该岩面糙率符合斯分布,由此得出了数值岩面,同时,用数学模型确定了裂隙面接触状态与应力的关系,用应力场与渗流场耦合地描述了裂隙流渗随应力变化的渗流特性,从计算途径验证了试验结果。     

管流与渗流耦合流动理论研究初探

针对缝洞型储层系统的特点,提出了研究这类流动规律的基本思路:将缝洞型储层系统分解为溶洞单元和裂缝单元,它们交替连接组成了整个储层流动系统;溶洞单元中的流动看成是管道流动,裂缝单元中的流动遵循达西定律.根据这个思路,研究了平面内缝洞耦合线性流动机理,建立了平面内线性流动的物理模型、数学模型及其数值求解方法.     

一个垂直平均水流运动的边界通用程式：Ⅰ.程式结构

采用一种破开算子的特征线方法，将垂直平均二维水流运动方程转化为一系列一维方程式，从而由特征有限差分法求解。该模式适用于河流、河口和海湾的平面二维水流或潮流模拟。为使模式在这些水域模拟中具有边界通用性，设计了一个包含实际水域可能出现的各种不同类型边界（共２８种）的边界信息库，分别给出了计算域内点和各种边界的计算公式。     

Artificial Boundary Method for Calculating Ship Wave Resistance

The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of the physical domain in the potential flow problem, we introduce one vertical side as an artificial upstream boundary and two vertical sides as the artificial downstream boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the potential flow problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. The solution of the variational problem by the finite element method gives the numerical approximation of the potential flow around the ship, which was used to calculate the wave resistance. The numerical examples show the accuracy and efficiency of the proposed numerical scheme.     

天然气水合物渗流特征及其描述

天然气水合物渗流是一个多相多组分非等温的物理化学渗流过程,这个过程包含了相变、能量变化以及储层介质的物理性质的变化。提出了天然气水合物藏孔隙度、渗透率、饱和度等参数的定义方式。应用渗流力学理论,结合水合物开采方式,对其渗流过程进行了描述,建立了天然气水合物藏开采的基本数学模型,并探讨了目前天然气水合物渗流研究存在的问题,为下一步天然气水合物藏的物理模拟和数值计算提供了理论基础。     

节点虚流量法的某抽水蓄能电站地下厂房渗流场分析

利用求解有自由面渗流场的节点虚流量法,对某抽水蓄能电站施工期地下厂房渗流场分布进行三维有限元仿真模拟,着重研究施工期厂房区洞室的渗流特性,并对施工期厂房区渗控方案进行评价.研究表明,施工期间排水设施尚未形成时,各已经开挖的洞室成为排水的主要边界;在保证各洞室的渗透水能及时排出的提前下,渗流场可以满足施工环境要求;顶排型排水幕排水效果并不明显.     

潮流数值模拟中紊动粘性系数的研究

从流体运动基本方程出发，对Boussinesq假设的合理性进行了初步探讨，阐述了紊动粘性系数的物理含义及其在数值计算中所起的作用，就一内设正交丁坝的半封闭海域中的平面紊动粘性系数对二维潮流数值计算结果的影响进行了讨论。结果表明，平面紊动粘性系数对潮位和平均流速的影响较小；随着紊动粘性系数的增加，丁坝背流侧回波的强度和影响范围均明显增大，就风生流和环岛水流两个理论算例中的垂向紊动粘性系数对三维潮流数值计算结果的影响进行了讨论。结果表明，垂向紊动粘性系数越大，流速分层越不明显，流速的绝对值也越小；垂向紊动粘性系数的梯度变化则对水平流速的垂向分布起着决定性的作用。综合理论分析和数值试验的成果，垂向紊动粘性系数宜采用抛物线分布的结构型式。     

基于颗粒流方法的碎石渗透系数分析

渗透系数是级配碎石渗流特性的重要参数,通过颗粒流数值方法建立碎石土渗流计算模型,利用流动方程和水压方程以及碎石之间的应力变形响应进行耦合计算,得到相应级配碎石的渗流系数,数值计算结果表明:颗粒流数值分析方法(PFC)适用于模拟任意形状、大小的单元集合体,允许单元在受力过程中发生开裂和运动.颗粒流不需要给出材料的宏观参数,只需给出材料的细观参数,通过细观结果的累加,PFC可以模拟材料在宏观状态下的表象.给出的9种级配碎石的渗透性远大于规范要求的级配集料的渗透系数,渗透系数均满足工程要求,但综合比较,孔隙率较大的碎石3类较优.