Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

A machine learning approach for predicting computational intensity and domain decomposition in parallel geoprocessing

ABSTRACT

High performance computing is required for fast geoprocessing of geospatial big data. Using spatial domains to represent computational intensity (CIT) and domain decomposition for parallelism are prominent strategies when designing parallel geoprocessing applications. Traditional domain decomposition is limited in evaluating the computational intensity, which often results in load imbalance and poor parallel performance. From the data science perspective, machine learning from Artificial Intelligence (AI) shows promise for better CIT evaluation. This paper proposes a machine learning approach for predicting computational intensity, followed by an optimized domain decomposition, which divides the spatial domain into balanced subdivisions based on the predicted CIT to achieve better parallel performance. The approach provides a reference framework on how various machine learning methods including feature selection and model training can be used in predicting computational intensity and optimizing parallel geoprocessing against different cases. Some comparative experiments between the approach and traditional methods were performed using the two cases, DEM generation from point clouds and spatial intersection on vector data. The results not only demonstrate the advantage of the approach, but also provide hints on how traditional GIS computation can be improved by the AI machine learning.     

江汉盆地钻孔测井资料确定地应力最大水平主压应力方向

利用钻孔测井资料并运用地层倾角测量信息分析法，给出了江汉盆地地应力最大水平主压应力方向为ＮＥ６０～６５°     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

四川石棉独立碲矿床蚀变岩的地质地球化学研究

四川石棉大水沟独立碲矿床堪称世界首例，其矿体两侧的围岩蚀变较窄，且蚀变带与矿脉及未蚀变岩呈清楚而截然的接触关系。主要围岩蚀变为白云石化、白（绢）云母化、云英岩化及电气石化等。这些蚀变作用分别发生于１７７．７～１６５．１Ｍａ和９１．７１～８０．１９Ｍａ。蚀变作用过程中ＳｉＯ２大量析出，Ａｌ２Ｏ３、ＴｉＯ２、ＣａＯ、ＭｇＯ、Ｋ２Ｏ、Ｈ２Ｏ、ＣＯ２、Ｆｅ２Ｏ３、ＭｎＯ则有不同程度地带入。与此同时，ＲＥＥ在此过程中包集体迁出，Ｔｅ、Ｂｉ、ＡＳ等元素的含量则决定于蚀变作用及其强度，而与原岩无关。     

A reinterpretation of the crustal N‐isotope record: evidence for a 15N‐enriched Archean atmosphere?

A new compilation of N‐isotope and abundance data for metasedimentary rocks, and hyrdothermal micas that proxy for bulk crust, show systematic patterns. (1) δ¹⁵N values of kerogen in Precambrian cherts are more negative relative to siliciclastic counterparts, probably due to a mantle hydrothermal component. (2) There is a secular trend from average δ¹⁵N 15.3 ± 1.8‰ in Archean shales, through intermediate values in the Proterozoic, to Phanerozoic counterparts where δ¹⁵N averages +3.5‰. (3) Hydrothermal micas in metamorphic hydrothermal systems of Palaeozoic and Mesozoic age that proxy for crust have δ¹⁵N within the range of contemporaneous sedimentary rocks. (4) Hydrothermal micas track the secular trend of δ¹⁵N for kerogen from 2.7 Ga to the Phanerozoic. (5) Within Precambrian datasets δ¹⁵N does not increase with decreasing N content; accordingly, high δ¹⁵N values cannot stem either from metamorphism or form Rayleigh fractionation. (6) Previous studies show isotopic shifts during metamorphism are only +1 to +3‰ up to amphibolite facies. Values of 10–24‰ are attributed to a high δ¹⁵N Archean atmosphere, a residual signature of CI carbonaceous chondrites where δ¹⁵N is +30‰ to + 42‰.     

PROPAGATION OF 30—60 DAY LOW FREQUENCY OSCILLATIONS AND THEIR INFLUENCE UPON THE SUBTROPICAL WESTERLIES JET STREAM DURING NORTHERN HEMISPHERE WINTER

Based on daily ECMWF gridpoint data of two winters during 1981—1983 including an ENSOyear,propagation of low frequency oscillations(LFO)during Northern Hemisphere winters andtheir influences upon 30—60 day oscillations of the subtropical jet stream are studied with the sta-tistical methods as complex empirical orthogonal function(CEOF)and so on.Results show that inthe winter of a normal year(1981—1982),30—60 day oscillations in the subtropical zone aremainly in the northern and southern flanks of exit region of jet stream.In the ENSO year(1982—1983),they are mainly in the vicinity of entrance and exit regions of jet stream.Intraseasonalchanges of subtropical jet stream manifested themselves as latitudinal fluctuation or longitudinalprogression or regression of about 40 day period.There are marked differences between propagat-ing passages of low frequency modes responsible for changes of subtropical jet stream in the normalyear(1981—1982)and in the ENSO year(1982—1983).Changes of oscillation amplitude showobvious phases.In general,the one in late winter is stronger than that in early winter,strongestone occurs in February.     

Silicon limitation on primary production and its destiny in Jiaozhou Bay, China——Ⅳ：Study on cross-bay transect from estuary to ocean

The authors analyzed the data collected in the Ecological Station Jiaozhou Bay from May 1991 to November 1994, including 12 seasonal investigations, to determine the characteristics, dynamic cycles and variation trends of the silicate in the bay. The results indicated that the rivers around Jiaozhou Bay provided abundant supply of silicate to the bay. The silicate concentration there depended on river flow variation. The horizontal variation of silicate concentration on the transect showed that the silicate concentration decreased with distance from shorelines. The vertical variation of it showed that silicate sank and deposited on the sea bottom by phytoplankton uptake and death, and zooplankton excretion. In this way, silicon would endlessly be transferred from terrestrial sources to the sea bottom. The silicon took up by phytoplankton and by other biogeochemical processes led to insufficient silicon supply for phytoplankton growth. In this paper, a 2D dynamic model of river flow versus silicate concentration was established by which silicate concentrations of 0.028–0.062 μmol/L in seawater was yielded by inputting certain seasonal unit river flows (m³/s), or in other words, the silicate supply rate; and when the unit river flow was set to zero, meaning no river input, the silicate concentrations were between 0.05–0.69 μmol/L in the bay. In terms of the silicate supply rate, Jiaozhou Bay was divided into three parts. The division shows a given river flow could generate several different silicon levels in corresponding regions, so as to the silicon-limitation levels to the phytoplankton in these regions. Another dynamic model of river flow versus primary production was set up by which the phytoplankton primary production of 5.21–15.55 (mgC/m²·d)/(m³/s) were obtained in our case at unit river flow values via silicate concentration or primary production conversion rate. Similarly, the values of primary production of 121.98–195.33 (mgC/m²·d) were achieved at zero unit river flow condition. A primary production conversion rate reflects the sensitivity to silicon depletion so as to different phytoplankton primary production and silicon requirements by different phytoplankton assemblages in different marine areas. In addition, the authors differentiated two equations (Eqs. 1 and 2) in the models to obtain the river flow variation that determines the silicate concentration variation, and in turn, the variation of primary production. These results proved further that nutrient silicon is a limiting factor for phytoplankton growth. This study was funded by NSFC (No. 40036010), and the Director's Fund of the Beihai Sea Monitoring Center, the State Oceanic Administration.     

基于Chandra的活动星系核观测

活动星系核（AGN）是宇宙中最奇特的天体之一。它是真正意义上的全波天体，其中X射线波段的发射功率占到全波段功率的50％左右。AGN的X射线辐射研究涉及天体物理中的最基本问题，例如能量产生、辐射机制和宇宙论等，而Chandra X射线卫星的高分辨率图像和光谱对这一研究有着重要作用。以Chandra卫星的部分观测结果为例，简要介绍了几类不同类型AGN的X射线辐射研究进展：（1）宽吸收线类星体APM08279＋5255（z=3．91）的X射线谱分析，以及高红移类星体的观测概况；（2）Seyfert星系NGC 4151延展X射线发射问题的解决，及NGC 1068X射线辐射与光学波段的高激发态发射线（[OⅢ]λ5007）有很强相关性的发现；（3）6个BLLac天体样本的X射线环境分析；（4）射电星系X射线喷流的观测等。