Thermal and rheological structures of the Xisha Trough, South China Sea

The Xisha Trough, located in the northwest of the South China Sea (SCS) mainly rifted 30 Ma ago, has been a failed rift since the cessation of the seafloor spreading of the NW subbasin. Based on the velocity–depth model along Profile OBH-4 across the Xisha Trough, a seven-layer density–depth model is used to estimate density structure for the profile. The relationship between seismic velocity and radiogenic heat production is used to estimate the vertical distribution of heat sources in the lower crust. The 2-D temperature field is calculated by applying a 2-D numerical solution of the heat conduction equation and the thermal lithosphere thickness is obtained from the basalt dry solidus (BDS). The rheology of the profile is estimated on the basis of frictional failure in the brittle regime and power-law steady-state creep in the ductile regime. Rheological model is constructed for a three-layer model involving a granitic upper crust, a quartz diorite lower crust and an olivine upper mantle. Gravity modeling supports basically the velocity–depth model. The Moho along Profile OBH-4 is of relatively high heat flow ranging from 46 to 60 mW/m² and the Moho heat flow is higher in the trough than on the flanks. The depth of the “thermal” lithospheric lower boundary is about 54 km in the center, deepens toward two sides, and is about 75 km at the northern slope area and about 70 km at the southern Xisha–Zhongsha Block. Rheological calculation indicates that the two thinnest ductile layers in the crust and the thickest brittle layer in the uppermost mantle lie in the central region, showing that the Xisha Trough has been rheologically strengthened, which are mainly due to later thermal relaxation. In addition, the strengthening in rheology during rifting was not the main factor in hampering the breakup of the Xisha Trough.     

Computational modeling of cyclic mobility and post-liquefaction site response

Under seismic excitation, liquefied clean medium to dense cohesionless soils may regain a high level of shear resistance at large shear strain excursions. This pattern of response, known as a form of cyclic mobility, has been documented by a large body of laboratory sample tests and centrifuge experiments. A plasticity-based constitutive model is developed with emphasis on simulating the cyclic mobility response mechanism and associated pattern of shear strain accumulation. This constitutive model is incorporated into a two-phase (solid–fluid), fully coupled finite element code. Calibration of the constitutive model is described, based on a unique set of laboratory triaxial tests (monotonic and cyclic) and dynamic centrifuge experiments. In this experimental series, Nevada sand at a relative density of about 40% is employed. The calibration effort focused on reproducing the salient characteristics of dynamic site response as dictated by the cyclic mobility mechanism. Finally, using the calibrated model, a numerical simulation is conducted to highlight the effect of excitation frequency content on post-liquefaction ground deformations.     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.     

用ANSYS模拟结构-地基动力相互作用振动台试验的建模方法

本文以结构－地基动力相互作用振动台模型试验为基础，结合通用有限元软件ANSYS，对结构－地基动力相互作用体系进行有限元计算建模的一些问题作了研究，主要包括柔性土容器的模拟、土体动力本构模型的选用、土体与结构交界面上的状态非线性模型、网格划分、重力的考虑、结构中钢筋的处理以及对称性的应用等。文中给出了利用上述建模方法对结构－地基动力相互作用体系进行计算的一些加速度时程结果，并与试验结果相对照，吻合较好。通过计算分析，验证了简化处理方法的合理性和计算模型的可行性。     

An oases dynamics model and its application in the Manas oasis

This study applies theories of systems and system dynamics to ecological economic systems in the oasis environment, with the intention of finding the basic characteristics of the oasis environment and its feedback structures. This study explores the inner mechanisms of economic development in the oasis environment and its relationship to policy variables. This model is applied to economic development in the Manas oasis of the Xinjiang Uygur Autonomous Region of China. System dynamics method starts from studying feedback structure within the system study and then searches for the best solution through identifying the right "policy point" and through model simulation in the computer. It embodies the cream of "systems cybernetics theory". This study is justified because the nature of the oasis system and the simulation results are scientifically sound and practically feasible.     

基于OpenGL的多层三维空间实体叠加建模

多层三维空间实体普遍存在于现实生活中。借助于目前国际最广泛应用的三维图形接口OpenGL将其真实地显示出来，具有重要的现实意义。其中多层三维空间实体建模技术是显示过程的关键步鼻。在实现单层空间实体三维建模的基础上，利用OpenGL实现了多层三维空间实体叠加建模，叠加实现的过程及其属性数据的表示方法。     

不同电阻率测深方法对旁侧不均匀体的反映

利用三维有限元数值模拟,重点研究了旁侧三维不均匀体所产生的视电阻率的影响问题,展示了4种常规电阻率测深法在旁侧不均匀体上的异常分布特点;探讨了测深点正下方目标体异常与旁侧不均匀体异常的区别问题,指出旁侧不均匀体的影响是产生电阻率法假异常的主要原因之一.     

生态地理建模中的多尺度问题

本文在分析生态地理建模内涵的基础上,讨论了生态地理建模中的尺度转换问题、跨尺度相互作用问题、空间尺度与时间尺度的关联问题和多尺度数据处理问题.由于生态地理问题的非线性、生态环境的异质性和随机事件,简单的线性尺度转换方法远不能满足生态地理建模的要求.为了从根本上解决生态地理建模中的时空尺度问题,除需要运用微分几何学和等级理论等经典方法外,还需要引入格点生成法和网格计算等现代理论和技术手段.     

安塞油田坪桥水平井区沉积微相三维建模研究

以安塞油田坪桥水平井区为例,探讨了应用水平井和露头资料及标点过程(示性点过程)方法进行沉积微相三维建模的思路与方法,同时对研究区三角洲前缘分流河道的定量结构模式进行了深入的研究。研究表明,坪桥水平井区三角洲分流河道单砂体宽度一般为40～100m左右。砂体结构呈两种模式,其一为“迷宫状”,砂体离散地分布于湖相泥岩中,其二为“拼合板状”,多个砂体在侧向上叠置而形成连片分布的砂体,但部分单砂体间可被泥质或钙质薄层侧积层所隔挡。应用水平井资料和露头资料,可有效地提取三维建模所必须的地质统计学参数。研究表明,采用标点过程方法,可建立符合地质实际的沉积微相三维模型。     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.