Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Hydrothermal Mineralization on the Mesoproterozoic Passive Continental Margins of China： A Case Study of the Langshan-Zha＇ertaishan Belt, Inner Mongolia, China

Most ore-forming characteristics of the Langshan-Zha‘ertaishan hydrothermal exhalation belt, which consists of the Dongshengmiao, Huogeqi, Tanyaokou and Jiashengpan large-superlarge Zn-Pb-Cu-Fe sulfide deposits, are most similar to those of Mesoproterozoic SEDEX-type provinces of the world. The characteristics include: (1) All deposits of this type in the belt occur in third-order fault-basins in the Langshan-Zha‘ertaishan aulacogen along the northern margin of the North China Platform; (2) these deposits with all their orebodies hosted in the Mesoproterozoic impure dolomite-marble and carbonaceous phyllite (or schists) have an apparent stratabound nature; ores display laminated and banded structures, showing clear depositional features; (3) there is some evidence of syn-sedimentary faulting, which to a certain extent accounts for the temporal and spatial distribution and the size of the orebodies in all deposits and the formation of intrabed conglomerates and breccias; (4) they show lateral and vertical zonation of sulfides; (5) The Cu/(Pb Zn Cu) ratio of the large and thick Pb Zn Cu orebodies gradually decreases from bottom to top; and (6) barite is interbedded with pyrites and sometimes with sphalerite. However, some characteristics such as the Co/Ni radio of the pyrites, the volcanism, for example, of the Langshan-Zha‘ertalshan metallogenic belt, are different from those of the typical SEDEX deposits of the world. The meta-basic volcanic rock in Huogeqi, the sodic bimodal volcanic rocks in the Dongshengmiao and potassic bimodal-volcanic rocks with blastoporphyfitic and blasto-glomeroporphyritic texture as well as blasto-amygdaloidal structure in the Tanyaokou deposits have been discovered in the only ore-bearing second formation of the Langshan Group in the past 10 years. The metallogeny of some deposits hosted in the Langshan Group is closely related to syn-sedimentary volcanism based on the following facts: most of the lead isotopes in sphalerite, galena, pyrite, pyrrhotite and chalcopyrite plot on both sides of the line for the mantle or between the lines for the mantle and lower crust in the lead isotope composition diagram; cobalt content of some pyrites samples is much higher than the nickel content (Co/Ni= 11.91-12.19). Some volcanic blocks and debris have been picked out from some pyritic and pyrrhotitic ores. All Zn-Pb-Cu-Fe sulfide orebodies in these deposits occur in the strata overlying metamorphic volcanic rocks in the only ore-bearing second formation. In the Jiashengpan deposit that lacks syn-sedimentary volcanic rocks in the host succession only Pb and Zn ores occur without Cu ore, but in the Dongshengmiao, Tanyaokou and Huogeqi deposits with syn-sedimentary volcanic rocks in the host succession Cu ores occur. This indicates a relatively higher ore-forming temperature. The process of synsedimentary volcanic eruption directly supplied some ore-forming elements, and resulted in secular geothermal anomaly favorable for the circulation of a submarine convective hydrothermal system, which accounts for the precipitation of deep mineralizing fluids exhaling into anoxidic basins along the syn-sedimentary fault system in the Langshan-Zha‘ertai rift. The Dongshengmiao, Tanyaokou, and Huogeqi deposits hosted in the Langshan Group appear to be a transitional type of mineral deposit between SEDEX and VMS-types but with a bias towards SEDEX, while the Jiashengpan deposit hosted in the Zha‘ertai Group is of a characteristic SEDEX type. This evidence, together with other new discoveries of Mesoproterozoic volcanic rocks and the features of lithogeny and metallogeny of the Bayun Obo deposit in the neighborhood emphasize the diversity, complexity and uniqueness of the Mesoproterozoic Langshan-Zha‘ertal-Bayun Obo ore belt.     

Stream channel erosion in a rapidly urbanizing region of the US–Mexico border: documenting the importance of channel hardpoints with Structure‐from‐Motion photogrammetry

Urbanization can lead to accelerated stream channel erosion, especially in areas experiencing rapid population growth, unregulated urban development on erodible soils, and variable enforcement of environmental regulations. A combination of field surveys and Structure‐from‐Motion (SfM) photogrammetry techniques was used to document spatial patterns in stream channel geometry in a rapidly urbanizing watershed, Los Laureles Canyon (LLCW), in Tijuana, Mexico. Ground‐based SfM photogrammetry was used to map channel dimensions with 1 to 2 cm vertical mean error for four stream reaches (100–300 m long) that were highly variable and difficult to survey with a differential GPS. Regional channel geometry curves for LLCW had statistically larger slopes and intercepts compared with regional curves developed for comparable, undisturbed reference channels. Cross‐sectional areas of channels downstream of hardpoints, such as concrete reaches or culverts, were up to 64 times greater than reference channels, with enlargement persisting, in some cases, up to 230 m downstream. Percentage impervious cover was not a good predictor of channel enlargement. Proximity to upstream hardpoint, and lack of riparian and bank vegetation paired with highly erodible bed and bank materials may account for the instability of the highly enlarged and unstable cross‐sections. Channel erosion due to urbanization accounts for approximately 25–40% of the total sediment budget for the watershed, and channel erosion downstream of hardpoints accounts for one‐third of all channel erosion. Channels downstream of hardpoints should be stabilized to prevent increased inputs of sediment to the Tijuana Estuary and local hazards near the structures, especially in areas with urban settlements near the stream channel. Copyright © 2017 John Wiley & Sons, Ltd.     

Flood frequency analysis of Turkey using L‐moments method

The index flood procedure coupled with the L‐moments method is applied to the annual flood peaks data taken at all stream‐gauging stations in Turkey having at least 15‐year‐long records. First, screening of the data is done based on the discordancy measure (D_i) in terms of the L‐moments. Homogeneity of the total geographical area of Turkey is tested using the L‐moments based heterogeneity measure, H, computed on 500 simulations generated using the four parameter Kappa distribution. The L‐moments analysis of the recorded annual flood peaks data at 543 gauged sites indicates that Turkey as a whole is hydrologically heterogeneous, and 45 of 543 gauged sites are discordant which are discarded from further analyses. The catchment areas of these 543 sites vary from 9·9 to 75121 km² and their mean annual peak floods vary from 1·72 to 3739·5 m³ s^?1. The probability distributions used in the analyses, whose parameters are computed by the L‐moments method are the general extreme values (GEV), generalized logistic (GLO), generalized normal (GNO), Pearson type III (PE3), generalized Pareto (GPA), and five‐parameter Wakeby (WAK). Based on the L‐moment ratio diagrams and the |Z_dist|‐statistic criteria, the GEV distribution is identified as the robust distribution for the study area (498 gauged sites). Hence, for estimation of flood magnitudes of various return periods in Turkey, a regional flood frequency relationship is developed using the GEV distribution. Next, the quantiles computed at all of 543 gauged sites by the GEV and the Wakeby distributions are compared with the observed values of the same probability based on two criteria, mean absolute relative error and determination coefficient. Results of these comparisons indicate that both distributions of GEV and Wakeby, whose parameters are computed by the L‐moments method, are adequate in predicting quantile estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

Response of detrital zircon and monazite,and their U–Pb isotopic systems,to regional metamorphism and host‐rock partial melting,Cooma Complex,southeastern Australia

Progressive Early Silurian low‐pressure greenschist to granulite facies regional metamorphism of Ordovician flysch at Cooma, southeastern Australia, had different effects on detrital zircon and monazite and their U–Pb isotopic systems. Monazite began to dissolve at lower amphibolite facies, virtually disappearing by upper amphibolite facies, above which it began to regrow, becoming most coarsely grained in migmatite leucosome and the anatectic Cooma Granodiorite. Detrital monazite U–Pb ages survived through mid‐amphibolite facies, but not to higher grade. Monazite in the migmatite and granodiorite records only metamorphism and granite genesis at 432.8 ± 3.5 Ma. Detrital zircon was unaffected by metamorphism until the inception of partial melting, when platelets of new zircon precipitated in preferred orientations on the surface of the grains. These amalgamated to wholly enclose the grains in new growth, characterised by the development of {211} crystal faces, in the migmatite and granodiorite. New growth, although maximum in the leucosome, was best dated in the granodiorite at 435.2 ± 6.3 Ma. The combined best estimate for the age of metamorphism and granite genesis is 433.4 ± 3.1 Ma. Detrital zircon U–Pb ages were preserved unmodified throughout metamorphism and magma genesis and indicate derivation of the Cooma Granodiorite from Lower Palaeozoic source rocks with the same protolith as the Ordovician sediments, not Precambrian basement. Cooling of the metamorphic complex was relatively slow (average ～12°C/10⁶y from ～730 to ～170°C), more consistent with the unroofing of a regional thermal high than cooling of an igneous intrusion. The ages of detrital zircon and monazite from the Ordovician flysch (dominantly composite populations 600–500 Ma and 1.2–0.9 Ga old) indicate its derivation from a source remote from the Australian craton.     

影响低渗透油藏低速非线性渗流的实验研究

以膨胀性的粘土为主要研究对象,通过粘土浆体的流变性实验和填砂管渗流实验,并使用一种新研制的微流量测量工具,研究了低渗岩心在较低流速下的渗流特性;利用微尺度流动理论,Einstein粘度定律等进行了理论分析,对引起低渗透油藏低速非线性渗流的影响因素进行了初步的探索.结论认为,粘土的膨胀性以及其水化后形成浆体的流变性是引起低速非线性渗流的主要因素之一.     

基于VisIt的全球科学数据并行可视化——以大气温度场为例

针对当前大规模全球科学数据可视化中存在的单机可视化数据量有限、从底层开发并行可视化系统难度大等问题,该文基于分布式环境和VisIt,提出了一种简便、开放而又有效的大规模全球科学数据可视化方法.介绍了VisIt的体系结构及运行机制,给出了自定义数据的并行可视化方法;并基于NCEP数据集及全球空间格网,在小规模集群环境下实现了小粒度适应性球体退化八叉树格网(SDOG)下的全球大气温度场的并行可视化.VisIt的并行可视化性能测试结果表明:通过增加计算节点,VisIt能有效摆脱传统单机可视化对数据量的限制,可实现大规模全球科学数据的并行可视化.     

Geology and Geochemistry of the Bianbianshan Au‐Ag‐Cu‐Pb‐Zn Deposit,Southern Da Hinggan Mountains,Northeastern China

The Bianbianshan deposit, the unique gold-polymetal (Au-Ag-Cu-Pb-Zn) veined deposit of the polymetal metallogenic belt of the southern segment of Da Hinggan Mountains mineral province, is located at the southern part of the Hercynian fold belt of the south segment of Da Hinggan Mountains mineral province, NE China. Ores at the Bianbianshan deposit occur within Cretaceous andesite and rhyolite in the form of gold-bearing quartz veins and veinlet groups containing native gold, electrum, pyrite, chalcopyrite, galena and sphalerite. The deposit is hosted by structurally controlled faults associated with intense hydrothermal alteration. The typical alteration assemblage is sericite + chlorite + calcite + quartz, with an inner pyrite - sericite - quartz zone and an outer seicite - chlorite - calcite - epidote zone between orebodies and wall rocks. δ34 S values of 17 sulfides from ores changing from –1.67 to +0.49‰ with average of –0.49‰, are similar to δ34 S values of magmatic or igneous sulfide sulfur. 206Pb/204Pb, 207Pb/204Pb and 208Pb/ 204Pb data of sulfide from ores range within 17.66–17.75, 15.50–15.60, and 37.64–38.00, respectively. These sulfur and lead isotope compositions imply that ore-forming materials might mainly originate from deep sources. H and O isotope study of quartz from ore-bearing veins indicate a mixed source of deep-seated magmatic water and shallower meteoric water. The ore formations resulted from a combination of hydrothermal fluid mixing and a structural setting favoring gold-polymetal deposition. Fluid mixing was possibly the key factor resulting in Au-Ag-Cu-Pb-Zn deposition in the deposit. The metallogenesis of the Bianbianshan deposit may have a relationship with the Cretaceous volcanic-subvolcanic magmatic activity, and formed during the late stage of the crust thinning of North China.     

地基微波辐射计资料在降水分析中的应用

利用地基微波辐射计的观测资料研究了2016年8月昆明长水国际机场出现的两次不同性质雷暴过程在不同发展阶段的温湿参量分布和演变特征,结果表明：（1）天气尺度强迫下雷暴前40 min,温湿廓线表现为低层温度降低、水汽密度升高,雷暴结束后趋于稳定,该变化可反映对流活动的起止。而局地热力雷暴前8 h~40 min,低层温度持续升高,降水前半小时水汽密度迅速升高,反映出了大气层不稳定能量积蓄过程和对流爆发前的空中水汽变化。（2）积分水汽含量（IWV）和云内液态水路径（LWP）在降水前有明显的陡增。天气尺度强迫下雷暴降水期间,IWV和LWP随着温度的增减相互转化,低层的剧烈降温,加速了水汽凝结饱和到降落的转化,造成降水强度进一步增大,对应IWV减小,转化为LWP和地面降水。热力雷暴降水期间,温度均匀变化,IWV、LWP均匀增长,水汽转化量极小,IWV增量远高于实际降水量。（3）大气可降水量（PW）在两类降水过程前半小时显著上升,天气尺度强迫下雷暴降水前总增幅达14.6 mm,局地热力雷暴降水过程前PW均匀增长7.6 mm。

     

一次秋季局地强降水过程分析

对常规气象要素、卫星云图和T106物理量诊断产品进行分析,得出2001年10月2～3日发生在桂东南局地强降水成因及中尺度系统特征。