东昆仑东段布青山地区上二叠统格曲组物源分析及其构造意义

对东昆仑东段布青山得力斯坦地区出露的上二叠统格曲组砾岩层进行砾石成分、砾度统计及系统的LA-ICP-MS锆石UPb年龄谱分析。结果表明,格曲组砾岩层砾石成分以石英岩和花岗岩为主,硅质岩和基性岩次之,砂岩和灰岩较少,砾石的分散系数为1.54~2.02,该套砾岩为近源快速堆积的产物。砾岩碎屑锆石U-Pb年龄可分为3组:1早奥陶世—晚志留世年龄组为499~409Ma,峰值年龄为426Ma,对应早古生代末期原特提斯洋向北俯冲碰撞产生的一系列构造岩浆事件;2新元古代年龄组为744~619Ma,峰值年龄为744Ma,对应于全球Rodinia超大陆裂解事件;3古元古代年龄组为2443Ma,对应东昆仑地区古元古代构造岩浆热事件。结合碎屑锆石年龄及沉积学特征综合分析,花岗质砾石来源于北侧东昆仑造山带加里东期岩浆弧,沉积岩砾石则可能来自造山带早期的沉积地层,石英岩及其他变质岩砾石则多来自东昆仑基底变质岩系。综合判别,格曲组为一套沉积于活动大陆边缘环境的滨浅海相磨拉石建造,代表南侧古特提斯洋向北俯冲开始的构造阶段,是初始俯冲的沉积构造响应。     

Contribution of phytoliths to total biogenic silica volumes in the tropical rivers of Malaysia and associated implications for the marine biogeochemical cycle

The contribution of phytoliths to total biogenic silica (BSi) volumes in rivers worldwide, and the associated implications for the biogeochemical cycle, require in-depth study. Based on samples from rivers in Peninsular Malaysia, this project investigated the source and characteristics of BSi found in Asian tropical rivers, as well as the process of reverse weathering taking place in these fluvial systems. Results indicated that BSi samples collected in sediments consisted of phytolith, diatom and sponge spicules. Phytoliths, predominantly of the elongate form, comprised 92.8%–98.3% of BSi in the Pahang River. Diatom BSi in this river consisted mainly of pennatae diatoms, but represented a relatively small proportion of the total BSi volume. However, diatom BSi (predominantly of the Centricae form) was more prevalent in the Pontian and Endau Rivers with shares of 68.8% and 79.3% of the total BSi volumes, respectively, than Pahang River. Carbon contents of the BSi particulates ranged from 1.85% to 10.8% with an average of 4.79%. These values are higher than those recorded in other studies to date, and indicate that BSi plays a major role in controlling permanent carbon burial. This study suggests that phytoliths from terrestrial plants are the primary constituents of BSi in the rivers of Peninsular Malaysia, and therefore represent a significant proportion of the coastal silica budget.     

钻探-物探测井技术在覆盖区填图中的应用研究

钻探作为获取深部地质信息的最直接手段,可查明第四纪地层序列,确定基岩顶面埋深与起伏、验证推测模型、研究地质构造、圈定地质体间的重要接触关系,服务于覆盖区填图。物探测井是一种科学有效的井内探测手段,利用特定的岩性或岩性组合对应着特定的测井参数组合的特点,利用测井参数组合可识别地层的岩性、确定地层的深度与厚度,可进行地层划分对比等。钻探-物探测井技术的综合利用,将提高覆盖区勘探效率,降低勘探成本,为下一步开展的覆盖区填图工作提供技术方法支撑。     

油铀兼探的找矿思路在松辽盆地的应用——以中央拗陷区南部为例

以"油铀兼探"的找矿思路,带动松辽盆地中南部开展砂岩型铀矿调查选区及钻探验证工作.将异常强度250 API和450 API确定为划分潜在铀矿孔、潜在铀矿化孔和正常孔的关键技术指标,对中央拗陷区南部1572个钻孔进行异常筛查工作.排查发现潜在铀矿孔362个,潜在铀矿化孔223个.对191号钻孔和36号钻孔进行钻探验证,发现1个工业孔、1个异常孔.认为区内寻找砂岩型铀矿的有利地段为油田区反转构造形成的隆起区边部四方台组河流相砂体,并总结了中央拗陷区南部6种主要控矿因素.     

Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Hybrid Gaussian-cubic radial basis functions for scattered data interpolation

Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent; however, for the datasets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large datasets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

Influence of physico-mechanical properties of Indian dimension stones on cutting rate of diamond wire saw

Cutting performance of diamond wire saw is a key factor influencing mine planning, production scheduling, and equipment selection for dimension stone quarries. It is normally measured in terms of cutting rate. Rock samples collected from various granite and marble quarries in India were tested in laboratory to determine their physico-mechanical properties. Cutting rate of diamond wire saw was measured in the field studies during the actual cutting process in quarries. Using these laboratory determined properties and the cutting rate, a multiple linear regression model has been developed to predict the cutting rate of diamond wire saw. Physico-mechanical properties of rocks determined in laboratory are used as independent variables and cutting rate as predictor variable in the regression model. The study indicates that the cutting rate increases with a decrease in most of the hardness and strength parameters of rock. The final model is tested for its goodness of fit indicating a significant linear relation between cutting rate and physico-mechanical properties, namely tensile strength, slake durability index, and Cerchar hardness index with regression coefficient of 94%. The resulting model can be used suitably for different types of hard to medium hard and soft dimension stones. The generalized model for estimating the cutting rate becomes a handy tool for mining engineers to work out operating efficiency, expenses, planning etc. of the dimension stone block cutting.