深层黏性土变形破坏的影响因素分析

河南刘河煤矿矿二1煤层上为厚度20~40m的第四系覆盖层,其中黏土层平均厚度为6.6m。由于第四系孔隙水含水组的四个含水段均为富水性中等及以上含水层,所以该层黏性土在受到上覆高压影响后能否控制顶板突水显得尤为关键,是解决新生界含水层下的安全开采的关键问题。通过室内实验测得该矿深度约在150m的4个钻孔土样的物理参数和土颗粒级配分析,按实验结果模拟深层黏性土,考虑受力面积、含水率和厚度3个影响黏性土破坏的主要因素,利用正交实验原理安排实验。结果表明,受力面积和含水率越大,黏性土的破坏程度增大;黏性土的厚度越大,破坏程度相应减小;3个因素中,受力面积的影响最大;当黏性土含水率在10%~25%,受力面积越小,厚度越大,黏性土不易发生破坏。该研究对预防和保证矿山安全生产具有重要意义。     

Robust and reliable metamodels for mechanized tunnel simulations

The main objective of this paper is to construct a robust and reliable metamodel for the mechanized tunnel simulation in computationally expensive applications. To accomplish this, four metamodeling approaches have been implemented and their performance has been systematically evaluated through a comparative study utilizing pure mathematical test functions. These metamodels are quadratic polynomial regression, moving least squares, proper orthogonal decomposition with radial basis functions, and an extended version of the latest approach. This extended version has been proposed by the authors and named proper orthogonal decomposition with extended radial basis functions. After that, a system identification study for mechanized tunneling has been conducted through the back analysis of synthetic measurements. In this study, the best performing metamodel, that is the one suggested by the authors, has been employed to surrogate a complex and computationally expensive 3D finite element simulation of the mechanized tunnel. The obtained results demonstrate that the proposed metamodel can reliably replace the finite element simulation model and drastically reduce the expensive computation time of the back analysis subroutine.     

Time-Expanded Sampling Approach for Ensemble Kalman Filter:Experiment Assimilation of Simulated Soundings

In the Ensemble Kalman Filter (EnKF) data assimilation-prediction system, most of the computation time is spent on the prediction runs of ensemble members. A limited or small ensemble size does reduce the computational cost, but an excessively small ensemble size usually leads to filter divergence, especially when there are model errors. In order to improve the efficiency of the EnKF data assimilation-prediction system and prevent it against filter divergence, a time-expanded sampling approach for EnKF based on the WRF (Weather Research and Forecasting) model is used to assimilate simulated sounding data. The approach samples a series of perturbed state vectors from Nb member prediction runs not only at the analysis time (as the conventional approach does) but also at equally separated time levels (time interval is △t) before and after the analysis time with M times. All the above sampled state vectors are used to construct the ensemble and compute the background covariance for the analysis, so the ensemble size is increased from Nb to Nb+2M￡Nb=(1+2M)×Nb) without increasing the number of prediction runs (it is still Nb). This reduces the computational cost. A series of experiments are conducted to investigate the impact of △t (the time interval of time-expanded sampling) and M (the maximum sampling times) on the analysis. The results show that if △t and M are properly selected, the time-expanded sampling approach achieves the similar effect to that from the conventional approach with an ensemble size of (1+2M)×Nb, but the number of prediction runs is greatly reduced.     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

异心球形电容器电容的精确解

提出了导体系统电容问题求解的一种新方法，借助双球面坐标系和平行板电容公式，求出了异心球形电容器电容的解析式。此外，对另外几种类型异体系统的电容也作了讨论。     

STUDY ON THE MIGRANT LABOUR FORCE OF CHINA IN RECENT YEARS──A Case Study of Nanhai City of Guangdong Province

I.INT~IOWFOralOngthe,thesyStemofPermanentresidencempistrationhasdividedChineseampleintOthO8ewith"tirbanresidence"andthOSehaving"rUralresidence".TherationingsyStem,underwhichedibleoilandgrainwasdistributedonlytOthOSewith"urbanresidence"Calds,hasbeendondnant.Bothsystems,beingcentraltorestrictingthegrowthofurbanpopulation,havestagnatedrUral-tirbaningration.Duringthe30yearsbefore1978,thestatedidsPOnsorpopulationmovementsfromrUraltourbanareasbymeansofenrollingnewworkers,butingeneral,thiss…     

Phase Equilibria Constraints on the Chemical and Physical Evolution of the Campanian Ignimbrite

The Campanian Ignimbrite is a > 200 km³ trachyte–phonolitepyroclastic deposit that erupted at 39·3 ± 0·1ka within the Campi Flegrei west of Naples, Italy. Here we testthe hypothesis that Campanian Ignimbrite magma was derived byisobaric crystal fractionation of a parental basaltic trachyandesiticmelt that reacted and came into local equilibrium with smallamounts (5–10 wt%) of crustal rock (skarns and foid-syenites)during crystallization. Comparison of observed crystal and magmacompositions with results of phase equilibria assimilation–fractionationsimulations (MELTS) is generally very good. Oxygen fugacitywas approximately buffered along QFM + 1 (where QFM is the quartz–fayalite–magnetitebuffer) during isobaric fractionation at 0·15 GPa ( 6km depth). The parental melt, reconstructed from melt inclusionand host clinopyroxene compositions, is found to be basaltictrachyandesite liquid (51·1 wt% SiO₂, 9·3 wt%MgO, 3 wt% H₂O). A significant feature of phase equilibria simulationsis the existence of a pseudo-invariant temperature, 883 °C,at which the fraction of melt remaining in the system decreasesabruptly from 0·5 to < 0·1. Crystallizationat the pseudo-invariant point leads to abrupt changes in thecomposition, properties (density, dissolved water content),and physical state (viscosity, volume fraction fluid) of meltand magma. A dramatic decrease in melt viscosity (from 1700Pa s to 200 Pa s), coupled with a change in the volume fractionof water in magma (from 0·1 to 0·8) and a dramaticdecrease in melt and magma density acted as a destabilizingeruption trigger. Thermal models suggest a timescale of 200kyr from the beginning of fractionation until eruption, leadingto an apparent rate of evolved magma generation of about 10^–3km³/year. In situ crystallization and crystal settling in density-stratifiedregions, as well as in convectively mixed, less evolved subjacentmagma, operate rapidly enough to match this apparent volumetricrate of evolved magma production. KEY WORDS: assimilation; Campanian Ignimbrite; fractional crystallization; magma dynamics; phase equilibria     

Trace element and Nd isotope analyses of apatite in granitoids and metamorphosed granitoids from the eastern Central Asian Orogenic Belt: Implications for petrogenesis and post-magmatic alteration

We present in situ trace element and Nd isotopic data of apatites from metamorphosed and metasomatized (i.e., altered) and unaltered granitoids in the Songnen and Jiamusi massifs in the eastern Central Asian Orogenic Belt, with the aim of fingerprinting granitoid petrogenesis, including both the magmatic and post-magmatic evolution processes. Apatites from altered granitoids (AG) and unaltered granitoids (UG) are characterized by distinct textures and geochemical compositions. Apatites from AG have irregular rim overgrowths and complex internal textures, along with low contents of rare earth elements (REEs), suggesting the re-precipitation of apatite during epidote crystallization and/or leaching of REEs from apatite by metasomatic fluids. ε_Nd(t) values of the these apatites are decoupled from zircon ε_Hf(t) values for most samples, which can be attributed to the higher mobility of Nd as compared to Sm in certain fluids. Apatites from UG are of igneous origin based on their homogeneous or concentric zoned textures and coupled Nd-Hf isotopic compositions. Trace element variations in igneous apatite are controlled primarily by the geochemical composition of the parental melt, fractional crystallization of other REE-bearing minerals, and changes in partition coefficients. Sr contents and Eu/Eu* values of apatites from UG correlate with whole-rock Sr and SiO₂ contents, highlighting the effects of plagioclase fractionation during magma evolution. Apatites from UG can be subdivided into four groups based on REE contents. Group 1 apatites have REE patterns similar to the host granitoids, but are slightly enriched in middle REEs, reflecting the influence of the parental melt composition and REE partitioning. Group 2 apatites exhibit strong light REE depletions, whereas Group 3 apatites are depleted in middle and heavy REEs, indicative of the crystallization of epidote-group minerals and hornblende before and/or during apatite crystallization, respectively. Group 4 apatites are depleted in heavy REEs, but enriched in Sr, which are features of adakites. Some unusual geochemical features of the apatites, including the REE patterns, Sr contents, Eu anomalies, and Nd isotopic compositions, indicate that inherited apatites are likely to retain the geochemical features of their parental magmas, and thus provide a record of small-scale crustal assimilation during magma evolution that is not evident from the whole-rock geochemistry.     

基于CMIP6模式的黄河上游地区未来气温模拟预估

利用第六次国际耦合模式比较计划（CMIP6）提供的5个气候模式,并结合基于地面气象站的CN05.1气象资料,评估了CMIP6模式对黄河上游地区1961—2014年气温变化的模拟能力。基于7个共享社会经济路径及代表性浓度路径（SSP-RCP）组合情景,结合多模式集合平均预估了2015—2100年黄河上游地区年均气温和季平均气温的时空变化规律。结果表明：多模式集合平均能较好地模拟黄河上游地区历史平均气温的空间分布格局与年变化。7个未来情景一致表明,2015—2100年黄河上游地区年平均气温呈现波动上升趋势［0.03～0.82 ℃?（10a）^-1］。其中,低辐射强迫情景下（SSP1-1.9、SSP1-2.6及SSP4-3.4）气温先呈现增加趋势,21世纪中期到达增幅峰值,之后增温呈现放缓趋势;而中、高辐射强迫情景下（SSP2-4.5、SSP3-7.0、SSP4-6.0及SSP5-8.5）气温表现为持续上升态势。空间上,未来气温增幅显著的区域位于黄河上游西部地区;时间上,呈现夏季增温快,春季增温慢。四季增温的空间分布呈现出一致特征,表现为西部增温强于东部,北部增温强于南部。研究结果可为黄河流域水资源管理及气候变化的适应性研究提供科学依据。     

电感耦合等离子体发射光谱法测定铅锌矿中15个主次量元素

铅锌矿多元素同时测定中,样品预处理过程受多种不同因素的影响,单因素试验往往无法同时为多因素实验提供合理的优化方案。本文采用盐酸-硝酸溶矿,考察了样品预处理过程中混合酸配比、混合酸体积、消解时间及消解温度对消解效果的影响,利用正交试验设计及极差分析确定了最优消解条件:盐酸-硝酸混合酸配比3∶1,混合酸体积10 m L,消解时间120 min,消解温度90℃,用电感耦合等离子体发射光谱法同时测定铅锌矿中15个主次量元素(铅锌铜锰砷银铋镉钴镍镓铟钼锗锑)的含量。在最优条件下,用富铅锌矿石国家标准物质(GBW 07165)进行试验,大多数元素的精密度(RSD,n=12)和准确度小于5%,方法检出限为0.0019~0.048μg/g。基于铅锌矿主要是以硫化物形式存在,采用一系列硫化物国家标准物质验证方法的准确度及可行性,检测结果基本都在标准值的误差范围内。应用该方法分析湖南某矿区中大批量铅锌矿样品,标样质量统计合格率为100%,密码质量统计合格率大于97%。