STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

According to the research theory of improved black oil simulator, a practical mathematical model for C02 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance fast, streamline method is employed to solve this model as an alternative to traditional finite difference methods. Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method (BEM), an explicit total variation diminishing (TVD) method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less information need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an effective tool for early stage screening andmiscible oil field.reservoir dynamic management of the CO2 miscible oil field.     

Dynamic solution of hyperbolic cooling tower with nonsymmetric imperfections by finite element method

In this paper, an approximate method is established to solve the couple problem between different harmonic waves in terms of the needs of nonsymmetric imperfections of cooling tower. By means of this method, the frequencies and response of rotational shell with local geometric imperfections are analysed and calculated.     

复杂边界条件下渗流场流线分布研究

流线分布研究已受到油藏工程师们的普遍关注。本文从无限大油藏稳态渗流场基本解出发,结合边界元方法求解出复杂边界条件下稳态渗流场的势分布,以此为基础提出了流线场的生成办法,并给出了应用实例。本文方法的优点在于求解过程中将问题的维数降低了一维,减少了计算量;计算精度较高且具有一定的普遍应用性;适应于求解任意形状的包括定压、定流量或混合边界在内的组合边界问题。实例表明,利用本文方法产生的流线分布因能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

煤层瓦斯流动的固－气耦合数学模型及数值解法的研究

本文根据煤层瓦斯流动时，固体变形与瓦斯流动的基本理论，提出了耦合数学模型，并证明该耦合数学模型解的存在性和光滑性，给出了这一非线性数学模型的数值解法。     

复杂边界非均质渗流场流线分布研究

建立了考虑源(汇)影响的含有不渗透区域复杂边界条件下非均质油藏稳定渗流的数学模型。利用扰动边界元方法求解数学模型,获得了地层中任意一点的压力公式.在此基础上,提出了流线场的生成方法。绘制了考虑非均质性、复杂边界和不渗透区域影响的流线分布图,并分析了流线分布的特征。通过分析表明,渗流场的非均质性和不渗透区域的存在都对流线分布存在较大的影响。利用本文方法产生的流线分布图能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

ENHANCED OIL RECOVERY BY FLOODING WITH HYDROPHILIC NANOPARTICLES

In this paper, the mechanism of enhanced oil recovery using lipophobic and hydrophilic polysilicon （LHP） nanoparticles ranging in size from 10 to 500 nm for changing the wettability of porous media was analysed theoretically. A one-dimensional two-phase mathematical model considering the migration and adsorption of LHP and wettability change in reservoir rock was proposed, and a simulator was developed to quantitatively predict the changes in relative and effective permeability of the oil and water phases and the oil recovery in sandstone after water driving. Numerical simulations were conducted to study the distribution of the particle concentration, the reduction in porosity and absolute permeability, the LHP volume retention on pore walls and in pore throats along a dimensionless distance, and oil production performance. In conclusion, oil recovery can obviously be improved by flooding with hydrophilic nanometer powders though permeability declines for the retention of nanoparticles in porous media. It is suggested that an LHP concentration ranging from 0.02 to 0.03 is preferable to enhance oil recovery.     

二维复合材料结构模态分析的边界元法

将边界元方法用于分析二维复合材料结构的自由振动模态，利用特解处理体积力（惯性力）仅需静态基本解就可求解问题，对一各向同性悬臂梁，用该法得到的结果与用有限元或各向性边界元法得到的结果符合得很好，但该法可解各向异性问题，对层状复合材料简支梁，用该法得到了数值结果与用一维层状复合材料梁的理论解的比较表明，当结构的长厚比大于２０时，二者符合得很好，当结构的长厚比小于２０时，一维层状材料梁的理论将产生很大的     

冻土水热力耦合作用的数学模型及数值模拟

将冻土体视为空间弹性体,提出了土体在冻结过程中水分场、温度场、应力场三场耦合的一般数学模型,并给出了相应的离散方程及其解法,最后给出了数值算例,并与实测值比较,证明了该模型和算法的正确性。     

不可压饱和多孔弹性梁、杆动力响应的数学模型

基于多孔介质理论,首先建立了饱和多孔弹性杆件弯曲与轴向变形时动力响应的数学模型.其次,基于多孔弹性梁弯曲变形的数学模型,利用Laplace变换,分析了两端可渗透的饱和多孔弹性悬臂梁在自由端受阶梯载荷作用下的动静力响应,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等物理量随时间的响应曲线.发现不可压多孔弹性梁的拟静态响应亦存在Mandel-Cryer现象,多孔弹性梁的挠度具有与粘弹性梁挠度类似的蠕变特征,然而,其应力响应不同于粘弹性梁,随着时间的增加,梁拟静态响应的弯矩逐渐增加,并达到一个稳态值.这些结果有助于揭示植物根茎等力学行为的机理.     

应变局部化分析的嵌入强间断多尺度有限元法

固体材料的应变局部化行为是导致结构破坏失效的重要因素之一,开展相关数值模拟分析对于结构安全性评估具有重要意义.然而由于材料的非均质和多尺度特性,采用传统数值方法进行求解时通常需要从最小特征尺度离散求解的结构,这将大幅度增加计算规模和成本.针对这一问题,本文提出了一种基于嵌入强间断模型的多尺度有限元方法.该方法从粗细两个尺度离散求解模型,首先在细尺度单元上引入嵌入强间断模型来描述单元间断特性,所附加的跳跃位移自由度则通过凝聚技术进行消除,从而保持细尺度单元刚度阵维度不变.其次,提出了一种增强多节点粗单元技术,其可根据局部化带与粗单元边界相交情况自适应动态地增加粗节点,新构造的增强数值基函数可以捕捉细尺度间断特性,完成物理信息从细单元到粗单元的准确传递以及宏观响应的快速分析;再次,在细尺度解的计算中,将细尺度解分解为降尺度解与单胞局部摄动解,从而消除弹塑性分析时单胞内部的不平衡力.最后,通过两个典型算例分析,并与完全采用细单元的嵌入有限元结果进行对比,验证了所提出算法的正确性与有效性.     

城市河网洪水过程的一维数值模拟

基于邻接矩阵和邻接表概念,提出了新的拓扑结构,即利用net、edge和node对象实现河网拓扑信息及河道断面和结点信息的储存;基于虚设河段法,推导出城市河网内边界即集中入流(雨水管出口)、闸堰的追赶系数,从而实现了基于GIS平台的城市河网非恒定流一维数值模拟程序的通用性和可移植性.最后,将温州市区“9.4”暴雨期间的河网洪水过程的计算值与实测值进行了对比,结果表明本模型可较好地模拟城市河网的洪水过程.     

THEORETICAL RELATIONS OF BOUNDARY DISPLACEMENT DERIVATIVES AND TRACTIONS AT SINGULAR BOUNDARY POINT FOR 2D ISOTROPIC ELASTIC PROBLEMS

Although boundary displacement and traction are independent field variables in boundary conditions of an elasticity problem at a non-singular boundary point, there exist definite relations of singularity intensities between boundary displacement derivatives and tractions at a singular boundary point. The analytical forms of the relations at a singular smooth point for 2D isotropic elastic problems have been established in this work. By using the relations, positions of the singular boundary points and the corresponding singularity intensities of the unknown boundary field variables can be determined a priori. Therefore, more appropriate shape functions of the unknown boundary field variables in singular elements can be constructed. A numerical example shows that the accuracy of the BEM analysis using the developed theory is greatly increased.     

含缺陷结构的塑性极限分析

结合极限分析中的数学规划理论和有限元技术,提出了三维含缺陷结构极限分析的数学规划方法,并采用罚函数法引入塑性不可压条件.对于考虑多组独立变化载荷联合作用的情况,提出了加载路径射线辐射求解方案,并基于这种射线辐射状的加载路径,推导了多组载荷联合作用下结构塑性极限上限分析的数崐学规划格式,编制了相应的有限元程序.文中的数值结果表明了该方法的正确性与有效性.     

平面问题等价边界积分方程的三次边界轮廓法

基于弹性力学平面问题等的边界积分方程,给出了三次单元的边界轮廓法。根据平面问题解的复变函数表示,构造了三次形函数。给出了对于混合边值问题求解系统方程确定的边界轮廓方程配置和三次单元界轮廓法的实施。     

考虑介质微观非均质特征的乳状液渗流模型

油水乳化渗流是三元复合驱和热力采油过程中常见的现象,地层介质的微观孔隙结构特征对乳状液流动有着重要影响.现有描述乳状液渗流的理论模型都属于确定性方法,只能反映出介质孔隙结构的体积平均效果.当介质内部微观非均质性相比其尺寸不能被忽略时,采用确定性方法描述会与实验结果存在偏差.基于连续时间随机游走理论建立了描述乳状液渗流的随机理论模型.该模型引入反映液滴微观运动特征的跃迁时间和跃迁位移两个概率分布函数来反映多孔介质微观非均质特征.研究结果表明该模型能很好地刻画实验曲线中出现的与介质尺度相关的拖尾现象,可作为更一般的过滤模型.     

一种有限元模型修正中的参数选择方法

有限元模型修正面临的首要问题是待修正参数的选择.目前主要是基于灵敏度分析的结果进行参数选择.本文从结构固有频率分析的能量法出发,论证了某阶固有频率主要受到同阶振型中发生弹性变形部位刚度的影响,据此提出一种待修正参数的选择方法,以振型中弹性变形为依据,选择弹性变形较大部位的参数作为待修正参数.用一个梁模型仿真验证了该方法的正确性,并将其应用于系杆拱桥的有限元模型修正中,修正后有限元计算结果与试验结果的最大误差缩小至3.3%.     

NEURO FUZZY MODEL FOR PREDICTING THE DYNAMIC CHARACTERISTICS OF BEAMS

An adaptive neuro-fuzzy inference system （ANFIS） is introduced to predict the dynamic behavior of beams. The effects of axial forces and large displacements are considered in the analysis. A database of tests for the dynamic characteristics of beams is developed from the experimental tests. The responses of nonlinear vibration force for the single and multiple-stepped beams are calculated from the finite element method （FEM）, experimental tests and neuro-fuzzy model for comparison. The neuro-fuzzy model provides a general framework for the combination of neural networks and fuzzy logic. It is more flexible with more options of incorporating the fuzzy nature of the real-world system and is an useful estimation tool for the dynamic characteristics of beams. Therefore, ANFIS can be a useful tool for dynamic behaviour analysis of multiple-stepped beams subjected to axial loads and large displacement.     

MODEL VERIFICATION ON 3D TIDAL CURRENT ANALYSIS IN TOKYO BAY

The results of a research project to verify the newly improved multiple- level model for 3D tidal current analysis in Tokyo Bay are presented. The improved multiple-level model includes additional effects due to Coriolis force, river inflows and wind shear stresses. Furthermore, a new numerical treatment of the open boundary condition was applied which effectively eliminated the spurious reflective waves often generated by various numerical methods simulating free surface flows. The mean (time-averaged or residual) and tidal currents in Tokyo Bay were simulated as examples to demonstrate the validity and capability of the newly improved multiple-level model. A series of numerical experiments was conducted to carefully examine the tidal circulations affected by the forcing factors of Coriolis force, river inflows and wind shears, both individually and combined. The numerical results demonstrated that the effects of each forcing term are physically reasonable, with the wind shear effect being the most significant and the case including all forcing terms being in best overall agreement with the field data collected in Tokyo Bay by the Ministry of Transportation. This study has contributed not only to the verification of the newly improved multiple-level model but also to the enhancement of the accuracy of numerical simulations of three-dimensional flow in coastal waters by this model.     

变分多尺度方法在一维力学模型中的应用

将变分多尺度方法应用于一维缆索模型,导出受力缆索的宏观有限元模型并求得细观位移解析解,总结出变分多尺度方法应用于具体模型的关键点和缺陷. 假定刚度为常值,数值模拟一定边界和受力下的缆索,得到宏观和细观位移. 将细观与宏观位移叠加,相比于精确位移得出：细观位移可视为常规有限元模型的后验误差. 变分多尺度方法在一维力学模型中的成功应用,推进了其实用性,为其在更多力学及工程问题中的运用和发展提供了参考.     

孔洞裂纹问题的一种数值方法

研究孔洞与裂纹的相互作用问题,通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系,将原问题分解为承受远处载荷不含裂纹不含孔洞的均匀问题,和在远处不承受载荷但在裂纹面上和孔洞表面上承受面力的多孔洞多裂纹问题.于是,以应力强度因子作为参量的问题可以通过考虑后者来解决,而利用笔者提出的杂交位移不连续法,这种多孔洞多裂纹问题是容易数值求解的.算例说明该数值方法对分析平面弹性介质中孔洞与裂纹的相互作用既简单又有效.