大型双槽渡槽地震反应分析

为了得到比较精确的大型双槽渡槽结构的地震反应，根据渡槽结构复杂、影响动力因素较多等特点，采用8结点空间块体单元对渡槽槽身和槽墩进行有限元离散，分别用耦合自由度和多个弹簧单元分别模拟渡槽槽身横向拉杆、加劲肋和盆式橡胶支座，建立了大型双槽渡槽结构地层反应分析的有限元模型，计算了南水北调水利工程中的双洎河大型双槽渡槽的振动特性，并分别用反应谱法和时程分析法对该渡槽进行了地震反应分析，计算结果可为该渡槽的抗震设计提供参考。     

Practical VTI approximations: a systematic anatomy

Transverse isotropy (TI) with a vertical symmetry axis (VTI) often provides an appropriate earth model for prestack imaging of steep-dip reflection seismic data. Exact P-wave and SV-wave phase velocities in VTI media are described by complicated equations requiring four independent parameters. Estimating appropriate multiparameter earth models can be difficult and time-consuming, so it is often useful to replace the exact VTI equations with simpler approximations requiring fewer parameters. The accuracy limits of different previously published VTI approximations are not always clear, nor is it always obvious how these different approximations relate to each other. Here I present a systematic framework for deriving a variety of useful VTI approximations. I develop first a sequence of well-defined approximations to the exact P-wave and SV-wave phase velocities. In doing so, I show how the useful but physically questionable heuristic of setting shear velocities identically to zero can be replaced with a more precise and quantifiable approximation. The key here to deriving accurate approximations is to replace the stiffness a₁₃ with an appropriate factorization in terms of velocity parameters. Two different specific parameter choices lead to the P-wave approximations of Alkhalifah (Geophysics 63 (1998) 623) and Schoenberg and de Hoop (Geophysics 65 (2000) 919), but there are actually an infinite number of reasonable parametrizations possible. Further approximations then lead to a variety of other useful phase velocity expressions, including those of Thomsen (Geophysics 51 (1986) 1954), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Harlan (Stanford Exploration Project Report 89 (1995) 145), and Stopin (Stopin, A., 2001. Comparison of v(θ) equations in TI medium. 9th International Workshop on Seismic Anisotropy). Each P-wave phase velocity approximation derived this way can be paired naturally with a corresponding SV-wave approximation. Each P-wave or SV-wave phase velocity approximation can then be converted into an equivalent dispersion relation in terms of horizontal and vertical slownesses. A simple heuristic substitution also allows each phase velocity approximation to be converted into an explicit group velocity approximation. From these, in turn, travel time or moveout approximations can also be derived. The group velocity and travel time approximations derived this way include ones previously used by Byun et al. (Geophysics 54 (1989) 1564), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Tsvankin and Thomsen (Geophysics 59 (1994) 1290), Harlan (89 (1995) 145), and Zhang and Uren (Zhang, F. and Uren, N., 2001. Approximate explicit ray velocity functions and travel times for P-waves in TI media. 71st Annual International Meeting, Society of Exploration Geophysicists, Expanded Abstracts, 106–109).     

压裂井经济分析的优化方法

水力压裂是油气井增产的一项重要技术措施。压裂效果的经济评价，实际上是一个压裂设计的优化过程。研究经济问题即为研究压裂优化问题。提出了压裂井的施工费用计算方法和贴现率、净收入、净现值的计算公式，并对压裂井进行了经济分析，给出了压裂设计优化的步骤。建立了八屋—后五家户气田压裂经济评价模型，并应用到实际施工现场。使油气单井产量比过去增加2～3倍。     

电法在濒海及岛屿地区寻找地下淡水中的应用

采用电法的基本原理，研究各种不同岩性和成淡水的电性差畀，划分地下成淡水分界面。在严重缺水的濒海及岛屿地区找到了较好的地下淡水，解决了守岛部队吃水难的问题。     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

小湾电站高边坡系统锚固与排水的优化设计方法

据小湾电站地下厂房进水口开挖高边坡的岩体结构与工程地质特征,采用基于结构面网络模拟的随机楔体稳定分析方法对边坡的稳定性进行了三维分析。进一步运用风险分析理论对边坡系统锚杆与预应力锚索的布置进行了优化设计,最后运用结构面网络模拟理论对排水方案进行了优化设计。     

废弃防空洞及软弱地基的加固处理

针对工程填土地层及废弃防空洞,拟采用锤击沉管挤密碎石桩法进行加固处理。但在处理部分防空洞后,因震动大、噪音高,对周边居民和民居影响严重而停工。通过多次试验,摸索出钻孔夯填法并结合灌浆法回填防空洞,喷浆桩法加固软弱地基等经验,并对废弃防空洞的回填加固机理进行了简要分析。经检测,加固处理效果显著。至今,住宅楼建成4年使用正常。     

Researches on Cartographic Database-Based Interactive Three-Dimensional Topographic Map

With the development of computer graphics, the three-dimensional (3D) visualization brings new technological revolution to the traditional cartography. Therefore, the topographic 3D-map emerges to adapt to this technological revolution, and the applications of topographic 3D-map are spread rapidly to other relevant fields due to its incomparable advantage. The researches on digital map and the construction of map database offer strong technical support and abundant data source for this new technology, so the research and development of topographic 3D-map will receive greater concern. The basic data of the topographic 3D-map are rooted mainly in digital map and its basic model is derived from digital elevation model (DEM) and 3D-models of other DEM-based geographic features. In view of the potential enormous data and the complexity of geographic features, the dynamic representation of geographic information becomes the focus of the research of topographic 3D-map and also the prerequisite condition of 3D query and analysis. In addition to the equipment of hardware that are restraining, to a certain extent, the 3D representation, the data organization structure of geographic information will be the core problem of research on 3D-map. Level of detail (LOD). space partitioning, dynamic object loading (DOL) and object culling are core technologies of the dynamic 3D representation. The objectselection, attribute-query and model-editing are important functions and interaction tools for users with 3D-maps provided by topographic 3D-map system, all of which are based on the data structure of the 3D-model. This paper discusses the basic theories, concepts and cardinal principles of topographic 3I)-map,expounds the basic way to organize the scene hierarchy of topographic 3D-map based on the node mechanism and studies the dynamic representation technologies of topographic 3D-map based on LOD, space partitioning, DOL and object culling. Moreover. such interactive operation functions are explored, in this paper, as spatial query, scene editing and management of topographic 3D-map. Finally, this paper describes briefly the applications of topographic 3D-map in its related fields.     

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

闪电首次回击过程的光谱特性

用无狭缝摄谱仪，拍摄了4次云对地闪电首次回击过程390一660nm波长范围的光谱，并将原子结构的理论应用于闪电光谱的研究，用多组态Dirac—Fock方法，系统考虑相对论效应、电子关联、弛豫效应等重要贡献，对光谱的精细结构进行了计算，给出了谱线的波长、振子强度以及相应的上能级激发能量。通过与试验观测谱线的比较分析发现，大多数谱线对应的上激发能在23ev左右。首次回击过程中，闪电通道内NII离子的密度最大，并且大多数处于n＝3的激发态。