南秦岭勉略构造带高压基性麻粒岩变质作用及其锆石U-Pb年龄

勉县-略阳地区是勉略蛇绿构造混杂岩带的代表区段,本文在勉县北部徐家坪地区确定了主要矿物为Grt+Cpx+Pl和具有典型"白眼圈"反应结构的两类高压基性麻粒岩,分别对其进行细致的岩相学研究,并利用THERMOCALC3.33程序进行P-T视剖面图计算。一类高压基性麻粒岩的峰期矿物组合为Grt1+Cpx+Pl1+Qz,对应温压条件为T=800～860℃,P=12.4～14.6kbar,晚期退变质矿物组合为Grt2+Hbl+Pl2+Qz。另一类是具有典型"白眼圈"反应结构的高压基性麻粒岩,"白眼圈"结构中斜长石为富Na的钠-更长石,以此推断该高压基性麻粒岩早期矿物组合中含绿辉石,因此其变质峰期矿物组合可能为Grt1+Omp(?)+Qz或Grt1+Cpx(?)+Pl+Qz,对应温压条件分别为T=775～900℃,P>19.2kbar和T=750～850℃,P=16.5～19.8kbar;该岩石后期经历了以矿物组合Grt2+Opx+Hbl1+Pl1+Qz为代表的麻粒岩相及以Grt3+Hbl2+Pl2+Qz为代表的角闪岩相两期退变质作用。造成这两种高压基性麻粒岩峰期变质矿物组合及其温压条件存在差异的原因可能是岩石原始成分的不同。对高压基性麻粒岩及其中的浅色脉体分别进行了LA-ICP-MS锆石U-Pb年代学分析,得到高压基性麻粒岩214±11Ma的变质年龄及脉体215±5Ma的结晶时代,并结合锆石微量元素特征分析,认为214±11Ma的年龄值代表该高压基性麻粒岩角闪岩相的退变质时代,同时获得该高压基性麻粒岩原岩形成时代可能为477Ma。综合两件高压基性麻粒岩的P-T演化轨迹及变质时代,建立高压基性麻粒岩的P-T-t演化轨迹,据此反映秦岭造山带在印支期沿勉略构造带发生俯冲-碰撞造山过程。     

Fast imaging and inversion of on and off shore electromagnetic data

We consider thin layer algorithms for inversion of electromagnetic data acquired on large experimental grids. These algorithms do not involve targeted search of the model parameters. Instead, they determine an integrated characteristic of the model, i.e. heterogeneous conductance or transverse resistance, using a filtration technique developed for thin layer models, which adequately represent the geo-electric situation by a sequence of laterally homogeneous and heterogeneous layers. The corresponding filter reflects impedance/admittance relations between different components of the electromagnetic field in the stratified part of the model. Convolution of such filters with the data obtained from experimental measurements represents the most time consuming part of the computation. The convolution is carried out using a high performance algorithm, which makes the process of inversion extremely fast. The choice of a particular inversion algorithm is dictated by the specific geo-electric situation in the area of interest and the electromagnetic data subject to interpretation. Typically, most computations can be carried out in up to several minutes using a regular desktop or laptop PC.We restrict our consideration to two simplest algorithms for determination of the conductance of heterogeneous conductive layers (S-algorithm) and the transverse resistance of heterogeneous resistive layer (T-algorithm). The algorithms can be useful for interpretation of data acquired in on- and off-shore conditions with natural and controlled sources. Inversion can be jointly carried out for data sets collected for different source locations, at different frequencies, etc. Likewise, the algorithms can be used for interpreting the data collected using a moving source.Performance of the algorithms is demonstrated on a set of test inversions of numerically simulated data sets. The first group of the examples is typical for environmental applications. It is studied using the S-algorithm applied to natural and controlled source measurements for shallow and deep targets. The T-algorithm is applied for delineation of a gas hydrate type of reservoir using the data numerically simulated for such a problem.In geo-electrical situations that can be addressed using thin layer models, results of inversion hardly require any further elaboration. For instance, this happens when parameters of the heterogeneous layer satisfy thin layer conditions and its thickness can be determined from seismic measurements. When this is not the case, the model found by a thin layer inversion usually reveals correct location of the anomalous body and provides reasonable estimate of its integrated electric properties. In such situations the algorithm works as an imaging technique, which together with seismic information on the layer boundaries provides a justified initial model for a full scale 3D inversion and hopefully allows for avoiding pitfalls associated with a straightforward application of more traditional 3D inversion algorithms.     

Isotherms for the sorption of zinc and copper onto kaolinite: comparison of various error functions

In this research, the equilibrium sorption of Zn(II) and Cu(II) by kaolinite was explained using the Freundlich, Langmuir and Redlich–Peterson isotherms, via both linear and non-linear regression analyses. In the case of non-linear regression method, the best-fitting model was evaluated using six different error functions, namely coefficient of determination (r ²), hybrid fractional error function (HYBRID), Marquardt’s percent standard deviation (MPSD), average relative error (ARE), sum of the errors squared (SSE) and sum of the absolute errors (EABS). The examination of error estimation methods showed that the Langmuir model provides the best fit for the experimental equilibrium data for both linear and non-linear regression analyses. The SSE function was found to be a better option to minimize the error distribution between the experimental equilibrium data and predicted two-parameter isotherms. In the case of three-parameter isotherm, HYBRID was found to be the best error function to minimize the error distribution structure between experimental equilibrium data and theoretical isotherms. Non-linear method was found to be more appropriate method for estimating the isotherm parameters.     

Remediation of colloid-facilitated contaminant transport in saturated porous media treated by nanoparticles

Facilitation of contaminant transport in porous media due to the effect of indigenous colloidal fine materials has been widely observed in laboratory and field studies. It has been explained by the increase in the apparent solubility of low soluble contaminants as a result of their adsorption on the surface of fine particles. Attachment of colloidal fine particles onto the rock surface could be a promising remedy for this challenge. In this experimental study, the effect of five types of metal oxide nanoparticles, γ-Al₂O₃, ZnO, CuO, MgO, and SiO_2, on suspension transport was investigated. In several core flooding tests, different nanofluids were used to saturate the synthetic porous media. Subsequently, after sufficient soaking time, the suspension was injected into the treated porous media. Analysis of the effluent samples’ concentration by Turbidimeter apparatus demonstrated that the presence of nanoparticles on the rock surface resulted in a significant reduction in fine concentrations in the effluent samples compared with non-treated media; ZnO and γ-Al₂O₃ demonstrated the best scenarios among the tests performed in this study. In order to characterize the surface properties of the treated porous media, the zeta potential of the surface was measured. Results showed that the treated porous media acts as a strong adsorbent of fine particles, which are the main carrier of contaminants in porous media. These findings were quantitatively confirmed by calculation of the total energy of interaction between the fine particles and rock surface using the Derjaguin–Landau–Verwey–Overbeek theory.     

Perspectives of the Reduction of Nutrient Export from River Watersheds through the Introduction of Best Available Technologies for Agricultural Production: Based on Modeling Results

Water Resources - Estimate is made for the possible reduction of nutrient export from the drainage basins of tributaries of the Kuibyshev Reservoir through the introduction of the best available...     

Relative location and source mechanism of inland earthquakes in Northern Egypt

We have relocated 259 inland earthquakes in northern Egypt using the double-difference hypocenter technique. Among this dataset we are able to determine source mechanisms of 200 events using P-wave polarities and amplitude ratios as well. The studied earthquakes have been recorded by the Egyptian National Seismological Network from October 1997 to December 2006 with local magnitude (M_L) varies between 1.5 and 5.0. Three earthquake dislocations have been defined namely: Dahshour, southeast Beni-Suef, and Cairo-Suez district. Earthquake activities tend to occur in clusters along the first dislocation (Dahshour) however, relatively scattered along the second (southeast Beni-Suef) and the third (Cairo-Suez district) dislocations. At Dahshour dislocation three distinct clusters have been distinguished. Source mechanism solutions of Dahshour earthquakes displayed normal faulting with a strike-slip component to strike-slip faulting with a minor normal dip-slip component. Most of earthquake focal mechanism orientations are varying from NE-SW to NW-SE. The fault plane solutions of Beni-Suef earthquakes represented normal faulting with a strike-slip component. If the NNW-SSE striking plane has been chosen to be the actual fault plane, some solutions would indicated normal faulting with a sinistral strike-slip motion and other reflect normal faulting with a dextral strike-slip component. The fault plane solutions of Cairo-Suez district earthquakes are compatible with E-W to ENE-WSW striking normal fault with a dextral strike-slip motion.     

2020年龙泉山断裂带MS5.1地震成因综合分析

2020年2月3日四川省成都市青白江区发生M_S5.1地震,震中烈度为Ⅵ度。该地震事件震中位于龙泉山断裂带上,距离成都市中心38 km,是龙泉山断裂带历史上非常罕见的5.0级以上地震事件。针对该事件成因进行了综合分析与研究,具体内容包括:(1)通过收集历史地震资料讨论龙泉山断裂带的地震活动性;(2)利用高质量的波形数据对主震位置进行重定位;(3)根据地震层析成像获得的三维vP、vS以及泊松比(σ)模型分析了孕震构造和流体影响,以及(4)利用固体潮理论模型分析了固体潮与地震触发的相关性。结果表明,本次M_S5.1地震发生在龙泉山断裂带北段,震中坐标为(30.732°N,104.529°E),震源深度为15.12 km;震源位于高-低泊松比过渡带附近,并伴随着大范围的低速异常,初步推断与深部流体有关;同时,固体潮在断层面上产生的剪切应力变化,也可能与本次地震的触发密切相关,暗示着在地震发生前龙泉山北段的地震危险性已经达到了较高水平。因此深部流体侵入作用、强震同震效应以及特定孕震构造环境的综合影响可能是导致本次地震触发的主要因素。     

Refraction of the principal stress trajectories and the stress jumps on faults and contact surfaces: Part 1. Non-constrained regular trajectories

Rock masses contain ubiquitous multiscale heterogeneities, which (or whose boundaries) serve as the surfaces of discontinuity for some characteristics of the stress state, e.g., for the orientation of principal stress axes. Revealing the regularities that control these discontinuities is a key to understanding the processes taking place at the boundaries of the heterogeneities and for designing the correct procedures for reconstructing and theoretical modeling of tectonic stresses. In the present study, the local laws describing the refraction of the axes of extreme principal stresses T ₁ (maximal tension in the deviatoric sense) and T ₃ (maximal compression) of the Cauchy stress tensor at the transition over the elementary area n of discontinuity whose orientation is specified by the unit normal n are derived. It is assumed that on the area n of discontinuity, frictional contact takes place. No hypotheses are made on the constitutive equations, and a priori constraints are not posed on the orientation on the stress axes. Two domains, which adjoin area n on the opposite sides and are conventionally marked + and ?, are distinguished. In the case of the two-dimensional (2D) stress state, any principal stress axis on passing from domain ? to domain + remains in the same quadrant of the plane as the continuation of this axis in domain +. The sign and size of the refraction angle depend on the sign and amplitude of the jump of the normal stress, which is tangential to the surface of discontinuity. In the three-dimensional (3D) case, the refraction of axes T ₁ and T ₃ should be analyzed simultaneously. For each side, + and ?, the projections of the T ₁ and T ₃ axes on the generally oriented plane n form the shear sectors S ⁺ and S ^?, which are determined unambiguously and to whose angular domains the possible directions p ⁺ and p ^? of the shear stress vectors belong. In order for the extreme stress axes T ₁ ⁺ ,T ₃ ⁺ and T ₁ ^? , T ₃ ^? to be statically compatible on the generally oriented plane n, it is required that sectors S ⁺ and S ^? had a nonempty intersection. The direction vectors p ⁺ and p ^? are determined uniquely if, besides axes T ₁ ^? , T ₃ ^? and T ₁ ⁺ , T ₃ ⁺ , also the ratios of differential stresses R ⁺ and R ^? (0 ≤ R ^± ≤ 1) are known. This is equivalent to specifying the reduced stress tensors T _R ⁺ and T _R ^? The necessary condition for tensors T _R ⁺ and T _R ^? being statically compatible on plane n is the equality p ⁺ = p ^?. In this paper, simple methods are suggested for solving the inverse problem of constructing the set of the orientations of the extreme stress axes from the known direction p of the shear stress vector on plane n and from the data on the shear sector. Based on these methods and using the necessary conditions of local equilibrium on plane n formulated above, all the possible orientations of axes T ₁ ⁺ , T ₃ ⁺ are determined if the projections of axes T ₁ ^? , T ₃ ^? axes on side — are given. The angle between the projections of axes T ₁ ⁺ , T ₁ ^? and/or T ₃ ⁺ , T ₃ ^? on the plane can attain 90°. Besides the general case, also the particular cases of the contact between the degenerate stress states and the special position of plane n relative to the principal stress axes are thoroughly examined. Generalization of the obtained results makes it possible to plot the local diagram of the orientations of axes T ₁ ⁺ , T ₃ ⁺ for a given sector S ^?. This diagram is a so-called stress orientation sphere, which is subdivided into three pairs of areas (compression, tension, and compression-extension). The tension and compression zones cannot contain the poles of T ₃ ⁺ and T ₁ ⁺ axes, respectively. The compression-extension zones can contain the poles of either T ₁ ⁺ or T ₃ ⁺ axis but not both poles simultaneously. In the particular case when the shear stress vector has a unique direction p ^? on side ?, the areas of compression-extension disappear and the diagram is reduced to a beach-ball plot, which visualizes the focal mechanism solution of an earthquake. If area n is a generally oriented plane and if the orientation of the pairs of the statically compatible axes T ₁ ^? , T ₃ ^? and T ₁ ⁺ , T ₃ ⁺ is specified, then, the stress values on side + are uniquely determined from the known stress values on side ?. From the value of differential stress ratio R ^?, one can calculate the value of R ⁺, and using the values of the principal stresses on side ?, determine the total stress tensor T ⁺ on side +. The obtained results are supported by the laboratory experiments and drilling data. In particular, these results disclose the drawbacks of some established notions and methods in which the possible refraction of the stress axes is unreasonably ignored or taken into account improperly. For example, it is generally misleading to associate the slip on the preexisting fault with the orientation of any particular trihedron of the principal stress axes. The reconstruction should address the potentially statically compatible principal stress axes, which are differently oriented on opposite sides of the fault plane. The fact that, based on the orientation of the intraplate principal stresses at the base of the lithosphere, one cannot make a conclusion on the active or passive influence of the mantle flows on the lithospheric plate motion is another example. The present relationships linking the stress values on the opposite sides of the fault plane on which the orientations of the principal stress axes are known demonstrate the incorrectness of the existing methods, in which the reduced stress tensors within the material domains are reconstructed without allowance for the dynamic interaction of these domains with their neighbors. In addition, using the obtained results, one can generalize the notion of the zone of dynamical control of a fault onto the case of the existence of discontinuities in this region and analyze the stress transfer across the system of the faults.