基于BISQ模型双相各向同性介质中地震波数值模拟

基于BISQ机制,推导了双相各向同性介质中弹性波数值模拟方程,并借助交错网格有限差分方法进行了波场模拟。证实了三种波(快纵波、慢纵波、横波SV)的存在和波场特征;在合成地震记录中,观测列三种波的反射,由于各种波相互之间的转换,转换波在地震记录中显示也很清楚,这使得波场变得更为复杂化。     

基于ANSYS与OpenGL的高速公路结构物仿真建模及三维地层地表修剪方法

现代三维视景建模与仿真技术在高速公路工程选线、施工以及运营管理中具有良好的应用前景。结构物的仿真建模以及能反映施工前后地表形态变化的三维地层模型研究是高速公路信息化中的重要问题。介绍了一种借用ANSYS有限元分析软件提取结构物三维模型表面几何信息的简便方法;结合OpenGL纹理贴图技术,实现了结构物的仿真建模;针对结构物施工后对地表形态改变较大的问题,提出利用结构物表面数据对三维地层模型地表网格进行修剪的方法--地表网格竖直投影修剪法。把上述方法应用到广东省某高速公路高边坡的三维视景仿真建模中,结果表明,边坡的仿真模型较逼真。三维地层模型修剪后的地表形态与真实情况接近,验证了方法的合理性、可靠性。     

山东大尹格庄金矿床深部三维定量成矿预测

大尹格庄金矿是胶西北招平断裂带中段的大型金矿床,但随着逐渐的开采其保有资源量不断下降,急需在矿床的深边部区域开展接续资源找矿工作。本文以大尹格庄金矿为研究对象,以区域成矿理论和找矿勘查模型为指导,利用三维地质建模技术构建了相关地质体的三维模型,在此基础上,采用定量化的方法提取地质体与成矿相关的控矿因素指标和矿化指标,并建立了三维定量成矿预测模型,并应用该模型圈定了3个找矿靶区。研究表明,该方法能够适用于大中型矿山的深边部找矿工作,其三维及定量化的预测结果能够为后续的勘探及开采工作提供重要参考。     

火山岩储层建模初探

火山岩储层建模研究是一个世界性难题。近年来 ,国内外在含油气盆地的火山岩中发现了大量的油气藏。为了探索火山岩储层物性的空间分布特征 ,作者选取了松辽盆地徐家围子地区火山岩发育的营城组为研究目的层 ,以火山岩相分析为基础 ,以基于相分析的神经网络技术为手段 ,通过多种参数的选取和计算 ,尝试性地建立了该区火山岩储层地质模型。该模型表明 :层状火山机构比块状火山机构的物性好 ;多相火山机构比单相火山机构物性好 ;裂缝发育的火山机构比裂缝不发育的火山机构物性好。基底涌流相、火山空落相、火山沉积相是有利的储集相带。     

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.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

Fractal and Multifractal Modeling of Hydrothermal Mineral Deposit Spectrum： Application to Gold Deposits in Abitibi Area, Ontario, Canada

A number of fractal/multifractal methods are introduced for quantifying the mineral de-lmsit spectrum which include a number-size, grade-tonnage model, power spectrmn model,multi-fractal model and an eigeavalue spectrmn model The first two models characterize mineral deposits spec-tra based on relationships among the measures of mineral deposits.These include the number of deposits,size of deposits,concentration and volume of mineral deposits.The last three methods that deal with the spatial-temporal spectra of mineral deposit studies are all expected to be popularized in near future.A case study of hydrothermal gold deposits from the Abitibi area,a world-class mineral district is used to demonstrate the principle as well as the applications of methods proposed in this paper,It has been shown that fractal and multifractal models are generally applicable to modeling of mineral deposits and occurrences.Clusters of mineral deposits were identified by several methods including the power spectral eral deposits in the Timmins and Kirkland Lake camps.     

基于中等复杂地球系统模式的全新世 气候突变事件模拟试验*

全新世暖湿气候背景下的代表性气候突变是发生在约8.2kaB.P. 和4.2kaB.P. 的气候突然变冷变干的事件。对于气候突变事件的出现,海洋环流被普遍认为是其中的关键环节。科学家们利用不同类型的气候模式进行过个例和机理模拟试验,但还未能对气候突变机理做出合理的解释。本研究使用一个中等复杂程度的地球系统模式EMIC ( Earth-system Models Intermediate Complexity) MPM-2,模拟研究了全新世气候背景下北大西洋区域,淡水强迫机制作用对温盐环流以及全球气候的可能影响、气候系统响应和恢复时间特征以及南北半球的差异。模拟结果表明海洋对北大西洋淡水通量异常的强度响应敏感,海洋环流对热量和盐份的输送与淡水通量异常强度之间的平衡决定了温盐环流能否对淡水强迫产生不可逆转变,当淡水通量异常达到某一阈值后可能引起海洋和气候平衡态的改变。同时还发现,在气候变化过程中,两次淡水通量异常对海洋和大气系统的强迫作用是不同的,首次淡水通量异常对气候平衡态的影响要远大于其后的淡水通量异常。模拟结果还表明,北大西洋海洋的异常状况可以通过海洋、大气等各圈层相互作用传递到不同圈层和空间区域,在南北半球出现不同的响应特征。上述模拟结果对于认识全新世气候变化特征、不同突变事件之间的关系及其对现代气候的可能影响具有启发意义。     

青海共和盆地下白垩统烃源岩地球化学特征及其生油意义

我国西北大部分地区在早白垩世处于挤压为主的构造环境,但共和地区处于伸展的构造环境,发育了巨厚的下白垩统湖相沉积。我们首次在共和盆地下白垩统发现了一套２５８．１５ｍ 厚的烃源岩。这一发现改变了以往西北白垩系烃源岩仅分布在祁连－秦岭以北的观点, 这对于研究我国西北地区白垩系烃源岩的分布和油气田的形成具有重大意义。共和盆地下白垩统从下到上分为万一段、万二段、万三段、万四段。万一段和万四段是炎热干旱气候条件下形成的红色沉积, 基本没有生烃能力。万二段烃源岩是一套炎热潮湿气候条件下形成的浅湖－半深湖相沉积, 烃源岩有机质丰度达到了中到好烃源岩的丰度标准, 干酪根类型以Ⅰ２ 型为主。该烃源岩在Ｋ２ 早期（９０．０１Ｍａ） 进入生油门限, 开始生油, 从Ｋ２ 晚期（７４．１１Ｍａ） 以来,一直处于大量生油的成熟阶段。下白垩统万三段沉积时期, 水体时深时浅, 暗色泥岩和红色泥岩交替出现。暗色泥岩属于中等烃源岩, 其有机质类型主要为Ⅱ型。它在晚白垩世 （７８．０２Ｍａ） 进入生油门限, 但一直到第四纪 （１．５Ｍａ） 才进入大量生油的成熟阶段。共和盆地下白垩统烃源岩发育, 具有形成中小型油气田的可能性。     

Simultaneous transport of water and solutes under transient unsaturated flow conditions — A case study

The imbalance between incoming and outgoing salt causes salinization of soils and sub-soils that result in increasing the salinity of stream-flows and agriculture land. This salinization is a serious environmental hazard particularly in semi-arid and arid lands. In order to estimate the magnitude of the hazard posed by salinity, it is important to understand and identify the processes that control salt movement from the soil surface through the root zone to the ground water and stream flows. In the present study, Malaprabha sub-basin (up to dam site) has been selected which has two distinct climatic zones, sub-humid (upstream of Khanapur) and semi-arid region (downstream of Khanapur). In the upstream, both surface and ground waters are used for irrigation, whereas in the downstream mostly groundwater is used. Both soils and ground waters are more saline in downstream parts of the study area. In this study we characterized the soil salinity and groundwater quality in both areas. An attempt is also made to model the distribution of potassium concentration in the soil profile in response to varying irrigation conditions using the SWIM (Soil-Water Infiltration and Movement) model. Fair agreement was obtained between predicted and measured results indicating the applicability of the model.