微型抗滑桩单桩设计计算模型及算法研究

考虑微型桩与周围岩土体间摩擦力的作用,提出了一种新的微型抗滑桩单桩设计计算模型,并给出了具体算法。由微型桩加固滑坡体的变形特点,分析微型桩与岩土体之间的相互作用机制,将微型桩与周围岩土体的摩擦作用引入其受力分析中;根据微型桩上各部分受力特点的不同,将微型桩分成上部摩擦受拉段、中部滑坡推力作用段和下部锚固段3段进行分析,推导了微型桩总的变形控制方程及各分段的变形控制方程;采用初参数法对控制方程进行求解,得到了微型抗滑桩上的内力分布及变形规律。计算结果表明,在滑坡推力作用下,微型桩的变形主要发生在滑面附近及以上桩段,滑面附近桩段将产生较大的内力和弯曲变形,受拉段弯曲变形较小,近似水平移动;桩与岩土体间的摩擦力是微型桩与周围岩土体相互作用的重要组成部分,摩擦力的作用能显著减小微型桩的弯曲变形,有效控制滑坡体的位移。     

西安地区地裂缝带Q3原状黄土流变特性试验研究

汾渭盆地是我国乃至世界上地裂缝最为发育、灾害最为严重的地区。考虑在漫长的地质历史过程中黄土的流变特性对地裂缝破裂扩展模式的影响,采用分级循环加卸载方式,对西安地区长安地裂缝（fc1）带Q3原状黄土进行了不同围压下的固结不排水三轴流变特性试验,得到了地裂缝带Q3原状黄土的蠕变加、卸载曲线和应变速率与时间关系曲线。试验结果表明,不同围压下地裂缝带Q3原状黄土表现为低应力下为减速蠕变,高应力下为减速蠕变和等速蠕变的特征。并基于统一流变力学模型理论,建立了适用于西安地区地裂缝带Q3原状黄土的三维蠕变本构模型,获得了相应的本构模型参数。该成果为研究隐伏地裂缝破裂扩展模式及其力学机制奠定了基础。     

松散承压含水层水文地质参数分区及水流场数值模拟

为了揭示华北型煤田松散承压含水层水文地质参数及其对地下水流数值模拟的意义,以安徽淮北煤田宿南矿区祁东煤矿松散层承压第四含水层(简称“四含”)为研究示范,对多个影响因素综合分析,采用层次分析-模糊综合评价法,对研究区四含水文地质参数进行分区,合理确定各分区的水文地质参数,并用于采煤情景下地下水流动态数值模拟。数值模拟结果表明：祁东煤矿2008-2012年浅部煤层开采过程中四含地下水流场没有明显改变,但南北分区水头差逐年增大,地下水流向始终从南向北,而且在井田北部水力梯度逐渐减小,在井田南部水力梯度逐渐增大,南北区以倾角为60°~70°、断距为10~320 m的魏庙断层为分界线,基岩中大倾角与大落差断层的出现是引起上覆松散承压含水层水力梯度异常的根本原因。     

土壤重金属连续提取方法的优化

重金属在污染土壤中的形态分布决定着重金属的迁移性和危害的程度。土壤重金属形态分析应用最多的是Tessier和BCR连续提取法。Tessier和BCR连续提取法没有考虑土壤样品的特征。美国环保署危险废物浸出毒性鉴别标准法虽然考虑土壤pH值,但没有划分形态。采集不同pH土样,结合Tessier、BCR和毒性浸出鉴别方法的特征,优化出针对不同pH值土壤的连续提取方法,将土壤中重金属划分为活性态、次生碳酸盐结合态、次稳定态和稳定态。用优化的连续提取方法对土样进行连续提取,重金属回收率为85%~115%。优化连续提取方法的结果与Tessier和BCR结果对比显示优化连续提取方法克服了Tessier连续提取法对非石灰质土壤提取过量BCR连续提取法对石灰质土壤提取不足的缺点。优化后的连续提取方法数据稳定可靠,可作为重金属形态分析方法使用。     

基于选权拟合法的电离层电子密度层析重构

附加约束的电离层层析算法是解决电离层电子密度反演中不适定问题的主要方法,为避免此类方法中约束权阵的选取不当对电子密度分布重构产生的不良影响,本文将选权拟合法应用到电离层层析成像技术中.该方法特别设计了依据电子密度空间分布特性构造参数权矩阵的方案.新方法有明确的物理意义,挖掘了隐含的信息量,为解决电离层电子密度反演中由于观测数据的不足等因素引起的不适定问题提供了一种新途径,可以得到符合客观实际的结果.数值模拟实验和实测数据的反演结果证实了该算法有效性、可靠性和优越性.     

Two-dimensional physical-based inversion of confined and unconfined aquifers under unknown boundary conditions

An inverse method is developed to simultaneously estimate multiple hydraulic conductivities, source/sink strengths, and boundary conditions, for two-dimensional confined and unconfined aquifers under non-pumping or pumping conditions. The method incorporates noisy observed data (hydraulic heads, groundwater fluxes, or well rates) at measurement locations. With a set of hybrid formulations, given sufficient measurement data, the method yields well-posed systems of equations that can be solved efficiently via nonlinear optimization. The solution is stable when measurement errors are increased. The method is successfully tested on problems with regular and irregular geometries, different heterogeneity patterns and variances (maximum K_max/K_min tested is 10,000), and error magnitudes. Under non-pumping conditions, when error-free observed data are used, the estimated conductivities and recharge rates are accurate within 8% of the true values. When data contain increasing errors, the estimated parameters become less accurate, as expected. For problems where the underlying parameter variation is unknown, equivalent conductivities and average recharge rates can be estimated. Under pumping (and/or injection) conditions, a hybrid formulation is developed to address these local source/sink effects, while different types of boundary conditions can also exert significant influences on drawdowns. Local grid refinement near wells is not needed to obtain accurate results, thus inversion is successful with coarse inverse grids, leading to high computation efficiency. Furthermore, flux measurements are not needed for the inversion to succeed; data requirement of the method is thus not much different from that of interpreting classic well tests. Finally, inversion accuracy is not sensitive to the degree of nonlinearity of the flow equations. Performance of the inverse method for confined and unconfined aquifer problems is similar in terms of the accuracy of the estimated parameters, the recovered head fields, and the solver speed.     

A parallel dynamic programming algorithm for multi-reservoir system optimization

This paper develops a parallel dynamic programming algorithm to optimize the joint operation of a multi-reservoir system. First, a multi-dimensional dynamic programming (DP) model is formulated for a multi-reservoir system. Second, the DP algorithm is parallelized using a peer-to-peer parallel paradigm. The parallelization is based on the distributed memory architecture and the message passing interface (MPI) protocol. We consider both the distributed computing and distributed computer memory in the parallelization. The parallel paradigm aims at reducing the computation time as well as alleviating the computer memory requirement associated with running a multi-dimensional DP model. Next, we test the parallel DP algorithm on the classic, benchmark four-reservoir problem on a high-performance computing (HPC) system with up to 350 cores. Results indicate that the parallel DP algorithm exhibits good performance in parallel efficiency; the parallel DP algorithm is scalable and will not be restricted by the number of cores. Finally, the parallel DP algorithm is applied to a real-world, five-reservoir system in China. The results demonstrate the parallel efficiency and practical utility of the proposed methodology.     

鲁东地区视应力时空变化特征分析

利用山东数字地震台网记录到的山东东部地区2008年以来的地震波形资料,在去除地震序列的余震之后,采用Brune模型,并结合遗传算法反演计算了225 次Mr≥2. 0地震的视应力,讨论了视应力与震级的关系,分析研究了视应力的时空变化特征,并通过视应力的时空变化来研究该区应力状态的特征。结果表明:鲁东地区中小地震的视应力范围在0.06~3.91 MPa,平均应力水平在0.31 MPa,可作为该区的背景应力水平;视应力与震级呈正相关关系,随震级的增加而升高,利用稳态回归拟合初步得到适合于该区的定量统计关系;计算的视应力值减去由震级得到的理论视应力,去除震级的相关性干扰,得到差视应力,能更真实的反应该地区的应力变化特征;视应力的时间变化显示鲁东地区目前正处于应力的积累增强阶段,可能与近期鲁东地区震群频发有关;空间分布显示目前高视应力集中区主要分布在沂沭断裂带安丘段附近和烟台附近地区,对于山东地区年度危险区的判定有一定意义。     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.     

Structural modal parameter identification and damage diagnosis based on Hilbert-Huang transform

Traditional modal parameter identifi cation methods have many disadvantages,especially when used for processing nonlinear and non-stationary signals.In addition,they are usually not able to accurately identify the damping ratio and damage.In this study,methods based on the Hilbert-Huang transform(HHT) are investigated for structural modal parameter identifi cation and damage diagnosis.First,mirror extension and prediction via a radial basis function(RBF) neural network are used to restrain the troublesome end-effect issue in empirical mode decomposition(EMD),which is a crucial part of HHT.Then,the approaches based on HHT combined with other techniques,such as the random decrement technique(RDT),natural excitation technique(NExT) and stochastic subspace identifi cation(SSI),are proposed to identify modal parameters of structures.Furthermore,a damage diagnosis method based on the HHT is also proposed.Time-varying instantaneous frequency and instantaneous energy are used to identify the damage evolution of the structure.The relative amplitude of the Hilbert marginal spectrum is used to identify the damage location of the structure.Finally,acceleration records at gauge points from shaking table testing of a 12-story reinforced concrete frame model are taken to validate the proposed approaches.The results show that the proposed approaches based on HHT for modal parameter identifi cation and damage diagnosis are reliable and practical.