A FINITE ELEMENT SHALLOW WATER CIRCULATION MODEL FOR TOLO HARBOUR, HONG KONG

Numerical solution of two-dimensional shallow water equations has been achieved by the finite element method. In this work, the governing equations were first converted into different equations suitable for explicit, forward marching, integration in the time domain. A second order Taylor's series expansion was used for the continuity equation and an "upwind" differencing was incorporated in the nonlinear convective terms in the momentum equations.The model was applied to Tolo Harbour in Hong Kong. Velocity vector plots and surface elevations obtained from the calculations are reasonably consistent with the expected tidal response of the systems.     

瞬变电磁法三维模拟计算研究进展

瞬变电磁法应用于矿产资源、环境工程等领域.目前主流的瞬变电磁三维正演模拟方法包括积分方程法、有限差分法、有限体积法和有限元法.随着观测环境的复杂化以及探测精度要求的提高,有必要研究瞬变电磁法高精度三维模拟计算,以便推动数据处理解释方法的进步.本文系统介绍了瞬变电磁三维正演计算研究进展,分析了准静态差分方程的构建和发射源...     

Vector form Intrinsic Finite Element Method for the Two-Dimensional Analysis of Marine Risers with Large Deformations

Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static and dynamic analyses of marine risers. The proposed approach uses the vector form intrinsic finite element (VFIFE) method, which is based on vector mechanics theory and numerical calculation. In this method, the risers are described by a set of particles directly governed by Newton’s second law and are connected by weightless elements that can only resist internal forces. The method does not require the integration of the stiffness matrix, nor does it need iterations to solve the governing equations. Due to these advantages, the method can easily increase or decrease the element and change the boundary conditions, thus representing an innovative concept of solving nonlinear behaviors, such as large deformation and large displacement. To prove the feasibility of the VFIFE method in the analysis of the risers, rigid and flexible risers belonging to two different categories of marine risers, which usually have differences in modeling and solving methods, are employed in the present study. In the analysis, the plane beam element is adopted in the simulation of interaction forces between the particles and the axial force, shear force, and bending moment are also considered. The results are compared with the conventional finite element method (FEM) and those reported in the related literature. The findings revealed that both the rigid and flexible risers could be modeled in a similar unified analysis model and that the VFIFE method is feasible for solving problems related to the complex behaviors of marine risers.     

Adaptive finite element modeling of marine controlled-source electromagnetic fields in two-dimensional general anisotropic media

In this paper, we extend the scope of numerical simulations of marine controlled-source electromagnetic (CSEM) fields in a particular case of anisotropy (dipping anisotropy) to the general case of anisotropy by using an adaptive finite element approach. In comparison to a dipping anisotropy case, the first order spatial derivatives of the strike-parallel components arise in the partial differential equations for generally anisotropic media, which cause a non-symmetric linear system of equations for finite element modeling. The adaptive finite element method is employed to obtain numerical solutions on a sequence of refined unstructured triangular meshes, which allows for arbitrary model geometries including bathymetry and dipping layers. Numerical results of a 2D anisotropic model show both anisotropy strike and dipping angles have great influence on the marine CSEM responses.     

Joint inversion of magnetotelluric and fullwaveform seismic data based on alternating cross-gradient structural constraints

Magnetotelluric(MT)inversion and seismic inversion are important methods for the interpretation of subsurface exploration data,but separate inversion of MT and seismic produces ambiguous and non-unique results due to various factors.In order to achieve accurate results,the authors propose a joint inversion method of two-dimensional MT and seismic data in the frequency domain.The finite element method is used for numerical simulation of electromagnetic data in the forward modelling,and the Gauss-Newton method is used for the inversion.The 9-point-finite-difference method is used to solve the seismic wave field in the acoustic wave equation,and the inverse problem of seismic data is solved by full waveform inversion with a conjugate gradient,a simple and fast method.Cross gradient functions are used to provide constraint structure between resistivity and velocity parameters to carry out the joint inversion.The joint inversion algorithm is tested by double-rectangular model synthesis data,and the accuracy of the algorithm is verified.The results show that the joint inversion results are better than those from separate inversion.The algorithm is applied to a geophysical model of a metalliferous deposit in Jinchuan and is compared with the separate inversion results.It shows that the results obtained with joint inversion are much closer to the real model.     

平面地质图的三维地质体建模方法研究

长期以来,三维地质体建模研究主要依赖于钻孔、剖面等勘察数据进行精细化建模。然而,在缺乏大区域宏观地质背景的前提下,仅考虑局部地区的精细建模,往往会导致以偏盖全。鉴此,本文提出了一种以平面地质图为基础数据,以图切地质剖面为中介的三维地质体建模方法。该方法通过间接法加密了研究区域的产状数据,用三次样条拟合与人工修编结合的方法,对褶皱等复杂构造作处理,实现了平面地质图的图切剖面自动绘制方法;以规则体元为三维空间数据模型,通过自动绘制研究区域一系列相互平行的图切剖面,构建了栅格化的三维地质体模型;以星岗地区为研究区域,构建了原型实验系统,验证了该方法的有效性与实用性。实例证明,在其他地质数据匮乏的条件下,平面地质图的三维地质体建模方法是构建区域三维地质体模型的一种有效解决方案。     

INVENTOR用于斜齿轮造型和有限元分析

利用任意转角位置的渐开线齿轮齿廓数学模型在Inventor中进行编程二次开发,实现了斜齿轮的三维建模,并进行了有限元分析,结果是造型精度高,有限元分析较准确,易于掌握。     

特征约束六面体网格生成的改进栅格法

目前,针对地学对象的非结构化六面体网格生成问题仍未得到很好的解决。为满足相关地学研究对六面体网格的需求,促进地学模拟及传统GIS空间分析的进一步发展,本文对传统栅格法进行扩展,提出了一种面向地学研究的六面体网格生成的改进栅格法。通过提取实体模型表面的几何特征,考虑内部孔、洞约束生成骨干网格,在此基础上构建一套特征约束模板来处理复杂的特征约束,使之既能保持地学研究对象的几何形态特征,又能实现加密区域与非加密区域的平滑过渡,从而满足地学分析与模拟的需求。     

综合沉积正演与多点地质统计模拟辫状河三角洲: 以塔河T区为例

为有效提高三角洲沉积相模拟的精度,以塔河T区为例,应用地震、岩心及测井资料,通过定量评估可容纳空间、物源供给和沉积搬运之间复杂的关系,建立该地区辫状河三角洲砂体沉积正演模型,并将此转化为多点地质统计模拟的三维训练图像,进行研究区辫状河三角洲多点地质统计模拟。研究表明:辫状河三角洲砂体分布受沉积正演模拟控制参数影响,其中砂质供应体积分数、沉积物流入比例、水流载荷量及洪水期间隔均属于强敏感性参数,对研究区三角洲砂体分布影响较大;砂岩含量沉积正演模型体现了辫状河三角洲退积式沉积过程,符合研究区沉积特征,将其转化为三维训练图像,多点地质统计模拟结果在三维空间具有受训练图像约束的特征,体现了训练图像所反映的相带接触关系,并且与钻井认识一致。所提出的这种新的获取三角洲沉积体三维训练图像的方法,综合了沉积正演模拟与多点地质统计模拟的优势,取得了一定的应用效果,对类似沉积体三维地质建模具有一定的借鉴作用。      

高精度曲面建模的一种快速算法

 高精度曲面建模(HASM)是一种全新的曲面建模方法,其整个过程可分为偏微分方程的离散、采样方程建立和代数方程组求解3个阶段。目前所采用的求解对称正定方程组的方法主要是共轭梯度法。为了解决HASM的计算速度问题,本文给出了2种新的预处理共轭梯度算法,分别为不完全Cholesky分解共轭梯度法和对称逐步超松弛预处理共轭梯度法。实验表明,不完全Cholesky分解共轭梯度法收敛速度最快,且这2种预处理方法均比其他方法具有更快的收敛速度。     

���ڶ�߶ȱ�Ե�������쳣����

????λ??С????????????????????????????????????????????????Щ???????????????????????????????????????????????????????÷???????????????κμ??輴??????????β?????????????????н????????????????????????????????????????÷??????????????     

对流边界层理论中不稳定非自治系统的数值解（英文）

为了求解对流边界层理论中一个非自治微分方程系统,作者采用伽略金有限元方法,此方法是通过将无限区间上的三阶非线性微分方程转化成有限区间上的二阶微分形式,并构造出相应的伽略金有限元方程来求得数值解,该数值解与先前一些作者的结果一致,并且计算效率显高于其它数值方法.     

FAST APPROXIMATE INVERSION OF MAGNETOTELLURIC(FAIMT) DATA FOR A HORIZONTALLY LAYERED EARTH WITH NEAR-SURFACE INHOMOGENEITIES

IwrRODUcrIONNearsurface2-Dand3-DinhomogeneousdistortionsinMTdatamakeitdimculttointerpretMTdataevenwhentheregionalorbasinstruCturescanbetreatalas1-dfornsionalnyohmann,l975,Weide1t,l975,Wannamakeretal.,l984,LaTormcaetal.,l986,Zhang,l987,TorresWhrdinetal.,l992).Inthelasttwentyyears,2-Dand3-DmodelinghavebeenaanmPlishedbyusingnumricalmethods.Afteranalyringthesefor-wardmodelingresultS,wenowdeariyunderstandthemaorphySica1meCanisrnsasultingfromthesedistortionsofnear6urfaceinhomogeneities(Pa…     

地区级实景三维模型色彩均衡方法实验

自2020年起各省相继启动"实景三维中国"建设,小面积单一场景实景三维建模向地区级全域实景三维建模转变成为一种常态,因量变引起的小区域建模中易被忽略的技术问题也日益凸显,如三维模型成果色彩不均衡问题。由于采集设备、天气情况、地形地貌及航摄状态等诸多要素的影响,地区级实景三维模型成果整体的目视效果往往存在差异。本文以济宁市主城区实景三维建模项目为例,探讨地区级实景三维模型成果色彩均衡的技术方法。结果表明:采用基于模型纹理贴图色彩均衡的方法可以有效提高模型成果整体的目视效果,同时相较于基于原始影像预处理的模型色彩均衡方法,该方法效率提高了9倍。     

一种基于GTP的复杂地质体真3维建模方法:岩柱体分区法

在总结分析现有复杂断层及倒转褶皱3维建模方法的基础上,提出一种基于广义三棱柱(GTP)的3维交互建模方法(岩柱体分区RPBP建模法)。文中定义了岩柱体、界限三角形两个重要概念,RPBP法可以通过构建断层界限三角形,将复杂断层系转化为多个简单断层系,然后采用构建简单断层系模型的方法来实现复杂断层系的3维建模;而对于倒转褶皱,RPBP法则通过构建褶皱轴界限三角形避免数据局部重复和冲突,使GTP建模过程得以顺利进行。实验表明,该方法能便捷可靠地实现含多断层相互切割、倒转褶皱等复杂地质体的真3维建模。     

基于钻井数据的地质体3维可视化研究

地质体3维可视化是地质工程中急需解决的一个热点问题,利用可视化的结果,可以将地质体的结构特征直观形象地展现在用户面前,使之可以非常直观地、清楚地对地质体进行各种观察和分析。笔者在综合考虑各种建模方法的基础之上,结合本文的建模数据——钻井数据的特点,采用广义三棱柱(Generalized tri-prism,简称GTP)作为建模的体元对地质体进行3维建模,并对传统的基于钻井数据的3维地质建模的方法做了一些改进,使GTP可以自动地、快速地剖分成四面体,以此来解决由于GTP的侧面不是平面所带来的问题。最后,利用VTK(Visualization Toolkit,简称VTK)可视化技术,开发了一个3维地质可视化系统,并利用该系统对本文建模方法的实际建模效果进行了展示和验证。     

地壳块体应力场动力学演变过程仿真的初步研究

运用三维非线性动态有限元方法,研究华北地壳块体对渤海、黄海潮汐水位变化的动态响应,得出该块体中各力学量随时间演变的数值结果以及动态、直观的三维仿真映象,并以此探索计算力学与地学综合研究的途径.从仿真模拟结果可以看出:华北地区近期发生的几次较大地震均在应力较大或应力变化剧烈的地点.由于海洋动态载荷的量级小,计算结果的应力绝对值也很小,因此,这种吻合表明较小的动态载荷下的动态响应揭示了孕震构造的特征,显示了应力集中的地点.     

Two-dimensional numerical simulation of flow around three-stranded rope

Three-stranded rope is widely used in fishing gear and mooring system. Results of numerical simulation are presented for flow around a three-stranded rope in uniform flow. The simulation was carried out to study the hydrodynamic characteristics of pressure and velocity fields of steady incompressible laminar and turbulent wakes behind a three-stranded rope. A three-cylinder configuration and single circular cylinder configuration are used to model the three-stranded rope in the two-dimensional simulation. The governing equations, Navier-Stokes equations, are solved by using two-dimensional finite volume method. The turbulence flow is simulated using Standard κ-ε model and Shear-Stress Transport κ-ω(SST) model. The drag of the three-cylinder model and single cylinder model is calculated for different Reynolds numbers by using control volume analysis method. The pressure coefficient is also calculated for the turbulent model and laminar model based on the control surface method. From the comparison of the drag coefficient and the pressure of the single cylinder and three-cylinder models, it is found that the drag coefficients of the three-cylinder model are generally 1.3–1.5 times those of the single circular cylinder for different Reynolds numbers. Comparing the numerical results with water tank test data, the results of the three-cylinder model are closer to the experiment results than the single cylinder model results.     

面向GTP的三维地质模型空间剖切方法与应用

广义三棱柱（GTP）是近年来提出并被广泛应用于三维地质建模领域的一种较为成熟的空间数据模型,可满足大多数地质钻孔数据的三维地质建模的要求。然而,现有针对GTP模型的三维空间剖切分析方法,仍无法适用于偏斜钻孔数据所建三维地质模型高效多次任意剖切、不能支持较为复杂的地质模型空间分析的问题。本文通过对GTP剖切算法改进研究,提出了采用动态四面体剖分法减小数据冗余;鉴于GTP形态较为复杂,可对剖切后保留多面体进行拓扑关系的重组,解决任意多次剖切问题;将多次“面-体”形式的单剖切运算组合为一次“体-体”形式的复合剖切运算的方法,可高效实现对三维地质模型空间分析。本文着重面向GTP体元的多重任意切剖平面的剖切问题,在分析了GTP体元特征的基础上,弥补了三维地质建模领域中对GTP体元任意、多重剖切方法研究的不足。应用实例表明,上述改进算法可以提高GTP剖切计算速度,并能够快速实现巷道掘进模拟,以及空间开挖模拟等复杂的空间分析功能。     

3D GIS下地下水流有限元数值模拟方法与应用

目前缺乏针对较强空间特性水文地质计算参数的高效采集方法与管理机制,而且现有孔隙地下水流有限元数值模拟在三维可视化与空间分析等方面仍存在很多不足。因此,本文以有限元数值计算方法和三维地理信息平台为基础,结合GIS空间分析算法和计算机图形学理论,针对孔隙承压地下水有限元分析过程的关键步骤（概念模型构建、空间离散、水文地质参数提取与赋值等）,提出了孔隙地下水有限元数值模拟过程在3D GIS下的实现方法和技术框架。基于水文地质钻孔数据和GTP（Generalized Tri-Prism）空间数据模型,以盐城市滨海平原水文地质区为研究区域,实现了地下水系统空间结构特征的三维表达,同时与地下水流有限元数值模型相耦合,提出了3D GIS模型与地下水非稳定流有限元数值模型耦合的核心处理算法,并在此基础上开发了3D GIS地下水流数值模拟系统。应用实例表明,该技术框架可有效简化有限元分析流程、优化模型计算效率、实现地下水有限元数值模拟过程及计算结果的三维可视化。