内蒙古赤峰地区早二叠世花岗岩LA-ICP-MS锆石U-Pb年龄及地球化学特征   总被引：1，自引：1，他引：0  

李斌  陈井胜  刘淼  杨帆  李伟  吴振  陈敏华  武昌胜 《地质通报》2018,37(9):1671-1681

以赤峰地区铭山复式岩体、尖山子岩体为研究对象,通过锆石U-Pb年龄、地球化学研究,确定了其形成时代,探讨了岩石成因和岩浆源区性质及其形成的构造背景。LA-ICP-MS锆石~(206)Pb/~(238)U测年结果表明,铭山复式岩体中灰白色斑状黑云二长花岗岩、尖山子岩体中浅肉红色细粒正长花岗岩分别形成于284.4±7.9Ma、294.7±8.5Ma,为早二叠世。根据地质体间接触关系可知,铭山复式岩体中的二长花岗岩岩体和尖山子岩体均为复式岩体,有待进一步解体。斑状黑云二长花岗岩SiO_2含量较高,K_2O+Na_2O含量较高,A/CNK1,属于高钾钙碱性I型花岗岩。细粒正长花岗岩Al_2O_3含量为14.32%~15.14%,Na_2O/K_2O=0.71~0.99,A/CNK=1.17~1.20,在标准矿物中出现刚玉分子。二者微量元素特征相似,富集大离子亲石元素,亏损高场强元素Nb、Ta,具Sr、P、Ti的负异常,表现出岛弧岩浆岩类特征。结合岩石地球化学、区域地质特征,认为赤峰地区早二叠世岩体形成于岛弧/活动大陆边缘构造背景,其形成与古亚洲洋向南俯冲有关。  相似文献   

Discriminating characters of ore-forming intrusions in the super-large Chalukou porphyry Mo deposit,NE China     

Peixin Duan  Cui Liu  Xuanxue Mo  Jinfu Deng  Jinhua Qin  Yu Zhang  Shipan Tian 《地学前缘(英文版)》2018,9(5):1417-1431

The Chalukou porphyry Mo deposit, located in the Great Hinggan Range, is the largest Mo deposit in northeast China, although the age and genesis of the associated magmatic intrusions remain debated.Here we report zircon U-Pb ages and trace elements, whole rock geochemistry and Sre Nd isotope data with a view to understand the relationship between the magmatism and molybdenum mineralization.Zircon U-Pb analysis yield an age of 475 Ma for rhyolite in the older strata, 168 Ma for the premineralization monzogranite, and 154 Ma for the syn-mineralization granite porphyry. The granite porphyry and quartz porphyry are considered as the ore-forming intrusions. These rocks are peraluminous, alkali-calcic, and belong to high-K to shoshonitic series with a strong depletion of Eu. They also display characteristics of I-type granites. The rocks exhibit wide variations of(87 Sr/86 Sr)iin the range of 0.705426 -0.707363, and ε_(Nd)(t) of -3.7 to 0.93. Zircon REE distribution patterns show characteristics between crust and the mantle, implying magma genesis through crust-mantle interaction. The Fe_2O_3/FeO values(average 1) for the whole rock and EuN/Eu*Nvalues(average 0.45), Ce~(4+)/Ce~(3+) values(average 301)for zircon grains from the granite porphyry are higher than those from other lithologies. These features suggest that the ore-forming intrusions(syn-mineralization porphyry) had higher oxygen fugacity conditions than those of the pre-mineralization and post-mineralization rocks. The Chalukou Mo deposit formed in relation to the southward subduction of the Mongol-Okhotsk Ocean. Our study suggests that the subduction-related setting, crust-mantle interaction, and the large-scale magmatic intrusion were favorable factors to generate the super-large Mo deposits in this area.  相似文献   

The Fragmentation in a Protoplanetary Disk and the Properties of the Resultant Self-Gravitating Clumps     

Chunjian Liu  Zhen Yao  Lianhua Tian 《Astronomy Reports》2018,62(11):713-723

We study the fragmentation properties in the protoplanetary disk and properties of the resultant self-gravitating clumps using our newly constructed disk model. Our disk model includes the mass inflall term from a molecular cloud core and the photoevaporation winds effect. We adopt the conventional fragmentation criterion to judge whether a protoplanetary disk can fragment. In this work, we follow our previous work to investigate the properties of the resultant self-gravitating clumps. In our calculation, the initial masses of the resultant self-gravitating clumps lie in the range of tens of M_J to more than one hundred of M_J, where M_J is the Jupiter mass. These initial masses can seemingly account for the masses of extrasolar planets in magnitude. We also calculate the subsequent gas accretion of clumps in 1.27 × 10⁴ yr after the formation of self-gravitating clumps. We find that the subsequent gas accretion of self-gravitating clumps is very efficient, and the clump masses grow to hundreds of M_J and the physical radii R_c of clumps increase to about 10 AU. Additionally, we also calculate the orbital migration of clumps. We find that most clumps have short migration timescale to be accreted onto the protostar, and only a small fraction of clumps have long migration timescale (>10⁶ yr) to successfully become gas giant planets. These results are consistent with previous studies.  相似文献   

Adjoint analysis of Buckley-Leverett and two-phase flow equations     

Savithru Jayasinghe  David L. Darmofal  Marshall C. Galbraith  Nicholas K. Burgess  Steven R. Allmaras 《Computational Geosciences》2018,22(2):527-542

This paper analyzes the adjoint equations and boundary conditions for porous media flow models, specifically the Buckley-Leverett equation, and the compressible two-phase flow equations in mass conservation form. An adjoint analysis of a general scalar hyperbolic conservation law whose primal solutions include a shock jump is initially presented, and the results are later specialized to the Buckley-Leverett equation. The non-convexity of the Buckley-Leverett flux function results in adjoint characteristics that are parallel to the shock front upstream of the shock and emerge from the shock front downstream of the shock. Thus, in contrast to the behavior of Burgers’ equation where the adjoint is continuous at a shock, the Buckley-Leverett adjoint, in general, contains a discontinuous jump across the shock. Discrete adjoint solutions from space-time discontinuous Galerkin finite element approximations of the Buckley-Leverett equation are shown to be consistent with the derived closed-form analytical solutions. Furthermore, a general result relating the adjoint equations for different (though equivalent) primal equations is used to relate the two-phase flow adjoints to the Buckley-Leverett adjoint. Adjoint solutions from space-time discontinuous Galerkin finite element approximations of the two-phase flow equations are observed to obey this relationship.  相似文献   

High-dimensional geostatistical history matching     

João Carneiro  Leonardo Azevedo  Maria Pereira 《Computational Geosciences》2018,22(2):607-622

Hydrocarbon reservoir modelling and characterisation is a challenging subject within the oil and gas industry due to the lack of well data and the natural heterogeneities of the Earth’s subsurface. Integrating historical production data into the geo-modelling workflow, commonly designated by history matching, allows better reservoir characterisation and the possibility of predicting the reservoir behaviour. We present herein a geostatistical-based multi-objective history matching methodology. It starts with the generation of an initial ensemble of the subsurface petrophysical property of interest through stochastic sequential simulation. Each model is then ranked according the match between its dynamic response, after fluid flow simulation, and the observed available historical production data. This enables building regionalised Pareto fronts and the definition of a large ensemble of optimal subsurface Earth models that fit all the observed production data without compromising the exploration of the uncertainty space. The proposed geostatistical multi-objective history matching technique is successfully implemented in a benchmark synthetic reservoir dataset, the PUNQ-S3, where 12 objectives are targeted.  相似文献   

Constraint energy minimizing generalized multiscale finite element method in the mixed formulation     

Eric Chung  Yalchin Efendiev  Wing Tat Leung 《Computational Geosciences》2018,22(3):677-693

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast.  相似文献   

Interior boundary-aligned unstructured grid generation and cell-centered versus vertex-centered CVD-MPFA performance     

Shahid Manzoor  Michael G. Edwards  Ali H. Dogru  Tareq M. Al-Shaalan 《Computational Geosciences》2018,22(1):195-230

Grid generation for reservoir simulation must honor classical key constraints and be boundary aligned such that control-volume boundaries are aligned with geological features such as layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. An unstructured grid generation procedure is proposed that automates control-volume and/or control point boundary alignment and yields a PEBI-mesh both with respect to primal and dual (essentially PEBI) cells. In order to honor geological features in the primal configuration, we introduce the idea of protection circles, and to generate a dual-cell feature based grid, we construct halos around key geological features. The grids generated are employed to study comparative performance of cell-centred versus cell-vertex control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulations using equivalent degrees of freedom. The formulation of CVD-MPFA schemes in cell-centred and cell-vertex modes is analogous and requires switching control volume from primal to dual or vice versa together with appropriate data structures and boundary conditions. The relative benefits of both types of approximation, i.e., cell-centred versus vertex-centred, are made clear in terms of flow resolution and degrees of freedom required.  相似文献   

Ensemble clustering for efficient robust optimization of naturally fractured reservoirs     

Zhe Liu  Fahim Forouzanfar 《Computational Geosciences》2018,22(1):283-296

In the development of naturally fractured reservoirs (NFRs), the existence of natural fractures induces severe fingering and breakthrough. To manage the flooding process and improve the ultimate recovery, we propose a numerical workflow to generate optimal production schedules for smart wells, in which the inflow control valve (ICV) settings can be controlled individually. To properly consider the uncertainty introduced by randomly distributed natural fractures, the robust optimization would require a large ensemble size and it would be computationally demanding. In this work, a hierarchical clustering method is proposed to select representative models for the robust optimization in order to avoid redundant simulation runs and improve the efficiency of the robust optimization. By reducing the full ensemble of models into a small subset ensemble, the efficiency of the robust optimization algorithm is significantly improved. The robust optimization is performed using the StoSAG scheme to find the optimal well controls that maximize the net-present-value (NPV) of the NFR’s development. Due to the discrete property of a natural fracture field, traditional feature extraction methods such as model-parameter-based clustering may not be directly applicable. Therefore, two different kinds of clustering-based optimization methods, a state-based (e.g., s _w profiles) clustering and a response-based (e.g., production rates) clustering, are proposed and compared. The computational results show that the robust clustering optimization could increase the computational efficiency significantly without sacrificing much expected NPV of the robust optimization. Moreover, the performance of different clustering algorithms varies widely in correspondence to different selections of clustering features. By properly extracting model features, the clustered subset could adequately represent the uncertainty of the full ensemble.  相似文献   

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

Halvor M. Nilsen  Jan Nordbotten  Xavier Raynaud 《Computational Geosciences》2018,22(1):233-260

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.  相似文献   

A numerical method for the solution of the nonlinear turbulent one-dimensional free surface flow equations     

Sujit K Bose 《Computational Geosciences》2018,22(1):81-86

Free surface flow of an incompressible fluid over a shallow plane/undulating horizontal bed is characteristically turbulent due to disturbances generated by the bed resistance and other causes. The governing equations of such flows in one dimension, for finite amplitude of surface elevation over the bed, are the Continuity Equation and a highly nonlinear Momentum Equation of order three. The method developed in this paper introduces the “discharge” variable q = η U, where η = elevation of the free surface above the bed level, and U = average stream-wise forward velocity. By this substitution, the continuity equation becomes a linear first-order PDE and the momentum equation is transformed after introduction of a small approximation in the fifth term. Next, it is shown by an invertibility argument that q can be a function of η: q = F(η), rendering the momentum equation as a first order, second degree ODE for F(η), that can be be integrated by the Runge-Kutta method. The continuity equation then takes the form of a first order evolutionary PDE that can be integrated by a Lax-Wendroff type of scheme for the temporal evolution of the surface elevation η. The method is implemented for two particular cases: when the initial elevation is triangular with vertical angle of 120 ^° and when it has a sinusoidal form. The computations exhibit the physically interesting feature that the frontal portion of the propagating wave undergoes a sharp jump followed by tumbling over as a breaker. Compared to other discretization methods, the application of the Runge-Kutta and an extended version of the Lax-Wendroff scheme is much easier.  相似文献