复杂模型重力位场快速正反演方法及应用

提要 本文详细讨论了一种三维重力位场快速正反演方法。作者在前人工作的基础上,对算法作了行之有效的改进,通过对反演中的不稳定因素进行各种理论模型试算,得出保证迭代反演稳定收敛的准则,编制出可在微型机IBM—PC上运行的人机对话式自动正反演程序。本文还对各种不均质模型进行了模似计算并将该方法应用于某含油气沉积盆地的双层界面构造研究,揭示出了储油有利地段。     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

Design of stable drainage networks in upland forestry plantations

Erosion of man-made, forestry drainage channels occurs when the plough cuts through the surface peat layer into the underlying erodible mineral soil. A procedure is developed, based on hydraulic considerations, which will allow the drainage engineer to design stable drainage networks in upland forestry plantations. An example design chart is given for an erodible sandy loam type soil.     

Dissipative-acoustic instability in accretion disks at a nonlinear stage

Nonstationary hydrodynamic models of a viscous accretion disk around a central compact object were constructed. Two different numerical methods (TVD and SPH) are used to study the dynamics of dissipatively unstable acoustic perturbations at the nonlinear stage in terms of the standard α-disk model. The standard disk accretion in the Shakura-Sunyaev model is unstable against acoustic waves for various parameters of the system. If the α parameter, which specifies the level of turbulent viscosity, exceeds α?0.03, then a complex nonstationary system of small-scale weak shock waves is formed. The growth rate of the perturbations is higher in the central disk region. For α?0.2, the relative shock amplitude can exceed 50% of the equilibrium disk parameters. The reflection of waves from the disk boundaries and their nonlinear interaction are important factors that can produce unsteady accretion. The luminosity of such a disk undergoes quasi-periodic oscillations at a level of several percent (?5%) of the equilibrium level.     

REVIEWS OF GEOGRAPHIC SOFTWARE

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)     

Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.     

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.