基于分形几何理论的DFN模型构建方法研究

以重庆梁-忠（梁平县-忠县）高速公路礼让隧道为工程研究背景,通过测线法调查现场节理,获得了节理产状分布概率密度函数,并从分形几何学的角度分析了节理间距及迹长的分形分布规律,推导出能反映节理间距及迹长分布状态的分形维数D及分形分布概率密度函数。在该基础上采用Matlab软件以及Monte-Carlo随机分析方法,产生节理参数随机数,结合3DEC中最新模块离散裂隙网络(DFN)技术,建立了能反映节理裂隙分布特征的离散裂隙网络模型并验证了模型的有效性,结果表明,分形分布比负指数分布包含更多的间距、迹长分布信息,更接近于实际分布;分形维数D反映了节理间距、迹长在其变化范围内的分布特征,分形维数的大小取决于小间距、小迹长部分数量在总节理数量中的比例,为节理裂隙岩体网络模型构建提供了一种新方法。     

Laboratory Estimation of Rock Joint Stiffness and Frictional Parameters

Numerical modeling of complex rock engineering problems involves the use of various input parameters which control usefulness of the output results. Hence, it is of utmost importance to select the right range of input physical and mechanical parameters based on laboratory or field estimation, and engineering judgment. Joint normal and shear stiffnesses are two popular input parameters to describe discontinuities in rock, which do not have specific guidelines for their estimation in literature. This study attempts to provide simple methods to estimate joint normal and shear stiffnesses in the laboratory using the uniaxial compression and small-scale direct shear tests. Samples have been prepared using rocks procured from different depths, geographical locations and formations. The study uses a mixture of relatively smooth natural joints and saw-cut joints in the various rock samples tested. The results indicate acceptable levels of uncertainty in the calculation of the stiffness parameters and provide a database of good first estimates and empirical relations which can be used for calculating values for joint stiffnesses when laboratory estimation is not possible. Joint basic friction angles have also been estimated as by-products in the small scale direct shear tests.     

淄博市张店南部地热异常成因

通过研究井温曲线类型,确定研究区纯传导地温梯度为2.5℃/100m;结合各单井地温梯度,圈定了地热异常区域;揭示了该区地热异常具有下段以对流为主、上段以传导为主的特点.通过热水水化学和同位素数据的综合分析表明,该区地下热水补给来源为大气降水,热水赋存于断裂裂隙系统中,且沿断裂裂隙垂向运动十分强烈.     

内蒙古西拉木伦断裂带韧性变形特征

西拉木伦断裂带在早二叠世末发生了韧性右行走滑,形成了包括韧性右行走滑带、韧性压扁带及弱韧性变形带在内的韧性变形带。据动力和运动分析,表明该韧性变形带为一普通剪切带,其中韧性右行走滑带在平面上以简单剪切为主,韧性压扁带在平面上以纯剪切为主。     

Isotopic studies of leaf water. Part 1: A physically based two-dimensional model for pine needles

Understanding of isotopic variations in leaf water is important for reconstruction of paleoclimate and assessment of global biochemical processes. We report here a study of isotopic distributions within a single needle of two pine species, Pinus resinosa Ait and Pinus strobes L., with the objective of understanding how isotopic compositions of leaf water are controlled by environmental and physiological variables. A 2D model was developed to simulate along-leaf isotopic variations and bulk leaf water isotopic compositions. In addition to variables common to all leaf water isotopic models, this 2D model also takes into account the specific geometry and dimensions of pine needles and the isotopic transport in xylem and mesophyll. The model can successfully simulate oxygen isotopic variations along a single needle and averaged over a leaf (bulk leaf water). The simulations suggest that isotopic composition of the bulk leaf water does not always depend only upon the average transpiration rate, which in turn raises questions about using leaf water isotopic values to estimate transpiration rates. An unsuccessful attempt to simulate along-needle hydrogen isotopic variations suggests that certain unknown biological process(es) may not have been incorporated into our 2D model, and if so, it calls for a reevaluation of all other models for hydrogen isotopic simulations of leaf water since they too lack these processes.Existing leaf water isotopic models are reviewed in this work. In particular, we evaluate the most frequently used model, the stomatal boundary layer model (also referred to as the Craig-Gordon model). We point out that discrepancy between the boundary layer model and the measured bulk leaf water seems to depend upon relative humidity. Using our 2D model, we show that this humidity dependency is a result of an interplay between environmental and physiological conditions: if the transpiration rate of plant leaves decreases with increasing relative humidity, our 2D model can reproduce the pattern of isotopic discrepancy between boundary layer model predictions and observations, enabling us to understand better the reason behind this discrepancy.     

基于电偶极子近似的多边形回线源瞬变电磁响应

针对大回线源中心装置在回线内场非中心位置处场点响应表达式复杂,视电阻率计算困难,只能套用中心回线装置计算公式,导致解释精度下降问题,给出了圆回线的响应公式。提出对发射回线沿边框进行电偶极子源分解,通过对各个边框的偶极子源产生的响应进行叠加;讨论了圆回线与圆内接多边形发射回线的等效性问题;完成了大回线内任意点的场值计算,并对基于电偶极子近似的大回线内场的分布特征进行了分析,得出了大回线源内非中心点探测的可行性结论。研究成果具有理论价值和实际指导意义。     

藏东三江构造带云南段α石英热活化ESR定年与新生代构造事件

运用α石英热活化ESR定年法(HAESRDQ)对三江构造带云南段49件样品定年结果的分析发现,该方法所获测年数据与地质事实及其他测年方法定年结果吻合,表明是一种普遍适用于中—新生代地质事件定年的方法。同时,证实三江构造带云南段是新生代陆内褶皱造山带,大量具“层控”特征的有色贵金属矿床的成矿时代也是新生代。     

莺歌海盆地成因及其大地构造意义

莺歌海盆地位于南海西北部 ,属 NW走向红河断裂带向南海海域的延伸。本文通过对盆地结构、沉降特征和构造 -沉积迁移过程的研究 ,提出莺歌海盆地在始新世—早渐新世期间属左旋扭张性断陷盆地 ,晚渐新世—早中新世的盆地演化阶段受到红河断裂带的左行剪切运动影响。莺歌海盆地的形成和演化历史反映了印藏板块碰撞过程对南海形成演化的影响历史     

永安盆地晚古生代—中—新生代沉积构造环境演化特征

永安盆地发育在前泥盆纪基底之上,位于华夏块体南缘,沿NE向政和—大埔断裂带展布。利用层序地层学方法,辅之以野外调查,研究了永安盆地的地层序列和岩石组合：① 对晚古生代至早三叠世地层格架及柱状剖面的分析表明,伴随海平面的升降,海相沉积序列具有阶段性,不同沉积相在盆地东、中、西区段分布不均,盆地沉积沉降中心大致位于龙岩和梅县地区;② 对横贯盆区的地质剖面及盆内辅助剖面的研究表明,中—新生代,盆内以断块作用为主,改造了晚古生代地层,代之以断褶复合地层;③ 有机碳质层为盆地的海陆变迁作了时域界定,并且其赋存部位和形态也在空间上表征了盆地构造演化的力学机制。通过分析构造事件、古地理及物源区、深部构造、边界构造等盆地要素,研究了盆地的沉积构造环境演化及其应力机制,结果表明：① 晚古生代,海进海退的快慢受控于不同时期构造事件的强弱变化,由此产生的沿岸隆起区的剥蚀以及拉张背景的山体风化剥蚀为盆地提供了物源,存在北、南两个物源区;中—新生代,盆地总体为火山活动背景下的山前及河湖相沉积环境。② 盆地具有隆起—伸展构造发育的特点,主要经历了华力西期海陆交互相巨厚沉积阶段,印支期稳定的台地向活动大陆边缘转变阶段并伴随挤压隆升,以及印支期后中—新生代由挤压向拉张机制转换的构造改造阶段。以上认识为中国东南部盆地的薄弱基础研究提供了新的基础参考信息。     

Sobol’ sensitivity analysis of parameters in the common land model for simulation of water and energy fluxes

The water and energy transfer of land surface is complex due to its large spatial and temporal variability. The modeling and simulation is an important means to study land water and energy transfer, but most selection and analysis of model parameters are empirical and qualitative. This paper has proposed a method of quantitatively identifying the most influential parameters of Common Land Model through Sobol?? sensitivity analysis. Considering sensible heat flux as the model output, the first order and total sensitivity indices of 25 model input parameters are estimated using an improved Sobol?? method. The simulated results are resampled using a bootstrapping method and the corresponding sensitivity indices are calculated. Confidence intervals for the bootstrapping sensitivity indices are estimated by using a percentile method. The results show that the parameters phi0 and porsl are the most important parameters, followed by ref(2,1), tran(2,1) and bsw. Five out of 25 parameters need to have an accurate evaluation, while the other parameters are fixed to a certain value. The sensitivity indices of parameters phi0 and porsl are decreasing after precipitation, while the sensitivity indices of parameters tran(2, 1) and ref(2, 1) are increasing after precipitation.