High Precision Time Domain Forward Modeling for Crosshole Electromagnetic Tomography

To improve the resolution of crosshole electromagnetic tomography, high precision of forward modeling is necessary. A pseudo-spectral time domain (PSTD) forward modeling was used to simulate electromagnetic wave propagation between two boreholes. The PSTD algorithm is based on the finite difference time domain (FDTD) method and uses the fast Fourier transform (FFT) algorithm for spatial derivatives in Maxwell's equations. Besides having the strongpoint of the FDTD method, the calculation precision of the PSTD algorithm is higher than that of the FDTD method under the same calculation condition. The forward modeling using the PSTD method will play an important role in enhancing the resolution of crosshole electromagnetic tomography.     

Analysis of rock drainage and cooling experiments for underground cryogenic LNG storage

Hydrogeological monitoring was conducted around a pilot cavern for underground cryogenic LNG (Liquefied Natural Gas) storage. The monitoring was mainly focused on the operation of a drainage and recharge system. After the operation of the drainage system commenced, the drainage rate decreased rapidly in the initial stages and then decreased gradually. Hydrogeological monitoring revealed that the rock drainage system operated effectively. During drainage, the water table was maintained below the cavern roof. The recharge for ice-ring formation was performed in two phases. The first phase involved the cessation of pumping in downward-drainage holes and the second involved the closure of upward boreholes. Since the water table was maintained below the cavern roof, artificial recharge was planned at first. However, it was not implemented due to heavy rainfall in the recharge stage. On the basis of hydrogeological monitoring and hydraulic tests, it was found that the fractures above the roof and on the right wall of the pilot cavern mainly affected seepage into the cavern and thermal variation due to the storage of liquid nitrogen. Thermal variation was examined by the thermometers installed around the pilot cavern. The cooling and thawing processes reveal the characteristics of thermal distribution in the rock and the 0 °C isotherm. The cooling phase lasted for six months, and the 0 °C isotherm progressed in time after the injection of liquid nitrogen into the cavern. The isotherm propagated up to about 4 m from the floor and the sidewall of the cavern and about 3 m from the cavern roof. The cooling rate of the rock mass above the cavern roof was lower than that of the other cavern sides due to the gaseous space in the upper part of the containment. The fractures were analyzed and considered for thermal modeling. A two-dimensional finite element analysis was performed to compare the field monitoring at the pilot cavern. The numerical modeling shows the distance between the ice ring and heat transfer pattern of the fractures around the pilot cavern. The propagation of the measured and calculated 0 °C isotherm reveals that the water-conveying joint on the right wall might affect thermal propagation through a thermal pipe.     

Constraining the stepwise migration of the eastern Tibetan Plateau margin by apatite fission track thermochronology

Granites sampled from Garzê-Litang thrust, Longmen Shan thrust, Garzê and Litang strike-slip faults in the eastern Tibetan Plateau have been analyzed with apatite fission track thermochronological method in this study. The measured fission track apparent ages, combined with the simulated annealing mod- eling of the thermal history, have been used to reconstruct the thermal evolutionary histories of the samples and interpret the active history of the thrusts and faults in these areas. Thermal history mod- eling shows that earlier tectonic cooling occurred in the Garzê-Litang thrust in Miocene (~20―16 Ma) whereas the later cooling occurred mainly in the Longmen Shan thrust since ~5 Ma. Our study sug- gests that the margin of eastern Tibetan Plateau was extended by stages: through strike-slip faults deformations and related thrusts, the upper crust formed the Garzê-Litang margin in the Miocene epoch and then moved to the Longmen Shan margin since ~5 Ma. During this process, the deformations of different phases in the eastern Tibetan Plateau were absorbed by the thrusts within them and conse- quently the tectonic events of long-distance slip and extrusion up to hundreds of kilometers have not been found.     

Calcium isotope fractionation in coccoliths of cultured Calcidiscus leptoporus, Helicosphaera carteri, Syracosphaera pulchra and Umbilicosphaera foliosa

Four species of marine calcifying algae, the coccolithophores Calcidiscus leptoporus, Helicosphaera carteri, Syracosphaera pulchra and Umbilicosphaera foliosa were grown in laboratory cultures under temperatures varying between 14 and 23 °C, and one species, C. leptoporus, under varying [CO₃²⁻], ranging from 105 to 219 μmol/kg. Calcium isotope compositions of the coccoliths resemble in both absolute fractionation and temperature sensitivity previous calibrations of marine calcifying species e.g. Emiliania huxleyi (coccolithophores) and Orbulina universa (planktonic foraminifera) as well as inorganically precipitated CaCO₃, but also reveal small species specific differences. In contrast to inorganically precipitated calcite, but similar to E. huxleyi and O. universa, the carbonate ion concentration of the medium has no statistically significant influence on the Ca isotope fractionation of C. leptoporus coccoliths; however, combined data of E. huxleyi and C. leptoporus indicate that the observed trends might be related to changes of the calcite saturation state of the medium. Since coccoliths constitute a significant portion of the global oceanic CaCO₃ export production, the Ca isotope fractionation in these biogenic structures is important for defining the isotopic composition of the Ca sink of the ocean, one of the key parameters for modelling changes to the marine Ca budget over time. For the present ocean our results are in general agreement with the previously postulated and applied mean value of the oceanic Ca sink (Δ_sed) of about − 1.3‰, but the observed inter- and intra-species differences point to possible changes in Δ_sed through earth history, due to changing physico-chemical conditions of the ocean and shifts in floral and faunal assemblages.     

三维地震勘探激发参数选择与分析

文中介绍了山区进行三维地震勘探的重要性及其方法原理。以山西晋城煤业集团某勘探区为实例进行实地勘探和试验,并对试验结果进行了详细的分析,进一步验证了该方法的可行性。     

线性粘弹介质中地震波场数值模拟

本文将以往两种粘弹介质中地震波模拟方法的优点结合起来，以模型理论和积分本构方程为基础，从理论上分析了模型对地震波场的影响；采用交错网格有限差分法对粘弹介质中的地震波进行数值模拟.数值计算结果表明该方法不仅便于计算，同时也便于从力学的角度来分析地震波的传播.数值计算结果与理论分析一致，说明这种方法可以更为有效地模拟粘弹介质中地震波的传播.     

对门源盆地"红层"形成时代的新认识

门源盆地位于青藏高原东北缘北祁连山腹地,为狭长的新生代走滑拉分山间宽谷盆地。其中充填一厚约407~960 m、干旱气候条件下的河湖相沉积"红层",长期以来各项研究工作均未能于其中采获任何生物化石。因此,尽管各家均认同其属"第三纪"无疑,但仍对其确切时代存在着较大的争论,而对其时代的确定有助于正确认识门源盆地形成时的构造背景,具有较为重要的意义。通过区域地层对比、"红层"形成气候条件与南侧西宁盆地及北邻酒泉—张掖盆地古气候环境所作的对比,深入探讨认为其时代属中中新世,并初步分析了门源盆地形成的新构造运动大陆动力学机制背景。     

三维地震数据可视化原理及方法

可视化是实现三维地震数据真三维解释的有效方法,本文介绍了三维地震数据可视化的基本原理,论述了用光线投射算法实现可视化的主要步骤,详细推导了合成可视化图像的计算公式,分析了阻光度曲线的物理意义及对可视化图像的调节作用,并用两个实际三维地震数据进行了试验.结果表明该方法能够由原始三维地震数据直接对地质体进行成像,能将分散的、孤立的信息互相联系起来,揭示隐藏在数据中的地质现象和规律.     

Analysis and modeling of the seasonal South China Sea temperature cycle using remote sensing

The present paper describes the analysis and modeling of the South China Sea (SCS) temperature cycle on a seasonal scale. It investigates the possibility to model this cycle in a consistent way while not taking into account tidal forcing and associated tidal mixing and exchange. This is motivated by the possibility to significantly increase the model’s computational efficiency when neglecting tides. The goal is to develop a flexible and efficient tool for seasonal scenario analysis and to generate transport boundary forcing for local models. Given the significant spatial extent of the SCS basin and the focus on seasonal time scales, synoptic remote sensing is an ideal tool in this analysis. Remote sensing is used to assess the seasonal temperature cycle to identify the relevant driving forces and is a valuable source of input data for modeling. Model simulations are performed using a three-dimensional baroclinic-reduced depth model, driven by monthly mean sea surface anomaly boundary forcing, monthly mean lateral temperature, and salinity forcing obtained from the World Ocean Atlas 2001 climatology, six hourly meteorological forcing from the European Center for Medium range Weather Forecasting ERA-40 dataset, and remotely sensed sea surface temperature (SST) data. A sensitivity analysis of model forcing and coefficients is performed. The model results are quantitatively assessed against climatological temperature profiles using a goodness-of-fit norm. In the deep regions, the model results are in good agreement with this validation data. In the shallow regions, discrepancies are found. To improve the agreement there, we apply a SST nudging method at the free water surface. This considerably improves the model’s vertical temperature representation in the shallow regions. Based on the model validation against climatological in situ and SST data, we conclude that the seasonal temperature cycle for the deep SCS basin can be represented to a good degree. For shallow regions, the absence of tidal mixing and exchange has a clear impact on the model’s temperature representation. This effect on the large-scale temperature cycle can be compensated to a good degree by SST nudging for diagnostic applications.     

INTEGRATION OF THE MODELS OF ANNAGNPS AND REMM TO ASSESS RIPARIAN BUFFER SYSTEM FOR SEDIMENT REDUCTION

The United States Department of Agriculture (USDA) Annualized Agricultural Non-Point Source Pollution model (AnnAGNPS) is used to help evaluate a watershed response to agricultural management practices to control water quality. However, AnnAGNPS version 3.5 does not contain features to estimate the effect of a riparian buffer (RB) system on water quality. The Riparian Ecosystem Management Model (REMM) is used to simulate the impact of riparian buffer systems on water quality. However, frequently the lack of measured upland loadings that are required by REMM simulation limits the application of REMM. To address this data gap, a study was conducted to integrate AnnAGNPS with REMM for RB system simulation. AnnAGNPS was used to simulate water and sediment loadings from an upland field into a three-zone RB system at the Gibbs Farm located in the Georgia coastal plain. These AnnAGNPS outputs were used as the inputs to REMM. REMM was used to simulate water and sediment movement along the riparian buffers. The AnnAGNPS simulated amount of annual runoff at the edge of the field was close to observed amounts (Nash-Sutcliffe efficiency of 0.92). It is believed that a substantial portion of sand was removed from the runoff one meter into the grass buffer where the samplers were located; therefore, sand was excluded from the AnnAGNPS simulation for comparison with observed sediment. Excluding sand, the AnnAGNPS predicted amount of annual sediment matches the observed amount fairly well (Nash-Sutcliffe efficiency of 0.46). In addition, based on evaluating the percent reduction of sediment at each zonal interface, the AnnAGNPS/REMM model well simulated the function of the RB system to reduce sediment.