海洋大地电磁揭示了西太平洋苏达海山的电阻率结构

为了研究海洋板块内火山的深部电阻率结构及其形成过程,在西太平洋苏达海山进行了海洋大地电磁测量。采用主流数据处理方法SSMT2000进行数据处理,分别对旋转后的实测数据以及根据实测数据得到的两组旋转不变量数据进行一维结构假设下的大地电磁响应测试,对响应最好的YX方向数据进行一维反演,并结合一维正演、地质资料对反演结果进行综合解释。反演结果表明,苏达海山区域地壳厚度为21.5 km左右;较厚的火山碎屑岩层表明苏达海山的构建过程以喷发性的岩浆活动为主,其侵入性的岩浆活动较弱。     

南黄海磁性基底特征分析和综合解释

本文利用南通幅区域地质调查的重磁异常数据,在重磁异常特征分析的基础上,采用一维频谱分析法对南黄海进行磁性基底反演。结果表明:南黄海属于扬子块体向海域的延伸,具有元古代结晶基底和古生代褶皱基底的双基底结构,其磁性基底深度在2—10km之间变化,是褶皱基底与磁性层最小埋深的结合,在坳陷和隆起上表现不同。通过利用综合地质和地球物理综合方法,对南黄海区域的四大构造单元华北—狼林地块、扬子—京畿地块、华南—岭南地块、东海大陆架进行重磁异常变化特征和磁性基底特征分析,为研究区断裂等区域构造特征解释提供了依据。     

基于GIS技术的河流污染动态模拟系统

本文论述了稳态均匀流场中的一维河流水质扩散模型,并在此模型的基础上,基于GIS二次开发技术开发了河流污染动态模拟系统。该系统以乐山市涌斯江支流段的地理底图及基本数据为例,实现了对河流污染情况的动态模拟。通过对该系统的一系列操作,可以对河段任意位置污染物浓度进行自动查询和修改,从而为河流水质污染的监测和管理提供可视化参考依据。     

基于动土压力影响的坝基动剪切模量反演研究

对于土石坝的坝基中土体应力-应变关系的求解,由于非自由地表以及土体所受水平正应力变化的影响,理想一维剪切梁模型方法并不适用。利用下覆饱和砂土地基的土石坝离心振动台试验,分析了不同输入地震峰值下非坝体覆盖区地基、坝中、坝趾处土体加速度的响应规律,得到了坝基中上部土层动土压力在水平向的变化特征,提出了考虑水平向动土压力影响的坝基土体剪应力-剪应变反演计算方法,并与利用理想一维剪切梁模型所绘制的应力-应变滞回圈对比,分析了动土压力对坝基土动剪切模量的影响。结果表明：由坝基底部至坝顶,加速度的放大效应基本呈线性增大;土压力在水平向的变化情况在剪应力-剪应变反演计算中不能忽略;考虑水平向动土压力的影响后,土体剪切模量的计算值减小。随着深度增加,动土压力对坝基中部土体剪切模量的影响比上部土体更大,且计算相对误差随着深度和输入地震波峰值加速度的增加而变大。     

考虑实测锤击数影响的动力触探杆长修正系数研究

现行规范中动力触探锤击数修正系数的确定需考虑实测锤击数影响,而业内对此少有研究。在一次人工做功击打时,动力触探落锤和探杆通过碰撞传递机械能,基于一维碰撞理论,若恢复系数不为0,则二者碰撞后速度不等,在重力和土阻力作用下会继续碰撞。推导了仅计算首次碰撞后探杆能量、计算首次碰撞后落锤与探杆能量之和及计算首次碰撞后有限次碰撞3种计算方案下的杆长修正系数计算方法。研究结果表明：在碰撞过程中和在一次人工做功击打的两次碰撞间是否考虑土体阻力作用,对杆长修正系数的计算都有影响。使用计算首次碰撞后落锤与探杆能量和的计算函数对现行规范中杆长修正系数拟合,得到重型和超重型决定系数分别为0.9935、0.9751。仅考虑两次碰撞间土体对探杆阻力的影响,使用计算首次碰撞后有限次碰撞的计算方法,当恢复系数≤0.3时,计算得到的杆长修正系数随锤击数增大而减小。     

基于Lévy飞行的粒子群算法在大地电磁反演中的应用

粒子群优化算法在大地电磁测深反演中相较于一般的线性反演算法具有多种优点。然而标准粒子群算法在多维优化问题中存在早熟问题,为此,采用基于Lévy飞行随机游走策略的优化粒子群算法来处理局部最优解,增加寻优能力。通过对地电模型的反演对比表明,改进后的粒子群算法相较于标准粒子群算法适应度值下降速度更快、寻优能力更好。最后将该算法应用于已知钻孔旁实测数据,结果较好,表明该算法具有较好的实用性。     

连续排水边界下考虑变荷载的软土一维非线性固结分析

土体往往具有非线性压缩特性,不同压缩特性下土体的固结规律存在差异。综合考虑土体非线性特性、变荷载作用以及连续排水边界条件建立了一维固结方程。采用无条件稳定的有限差分法和半解析法对固结方程进行了求解,两种解答方法的可靠性通过连续排水边界条件退化以及两种解答结果对比得到验证。在有限差分法解答的基础上,详细分析了界面参数、荷载参数以及非线性参数对土体固结的影响。结果表明：连续排水边界的界面参数取值越大,则超静孔隙水消散速率及土体沉降速率越大,但界面参数取值对最终沉降量没有影响;超静孔隙水压力在加载阶段逐渐增大,在恒载阶段逐渐消散;加载速率越大,则超静孔隙水压力峰值越大,并且土体固结速率越快,说明延长施工周期有利于降低超静孔隙水压力的影响;工程中要准确预测土体固结速率绝非易事,采用固结理论预测时需保证土体模型、边界条件以及土体计算参数等因素的准确性。     

大地电磁测深反演技术有效性对比试验

设计了一个简单的三维模型:三维低阻异常体和高阻异常体位于一维层状介质模型中,以"十"字型和"米"字型观测剖面方式作为模拟方式,开展模型一维、二维、三维反演技术有效性对比试验,其中一维反演计算采用自适应正则化(ARIA)反演,二维反演计算采用非线性共轭梯度(NLCG)反演,三维反演计算采用REBOCC三维反演,将剖面下方不同的反演结果与原始模型进行对比分析,发现:在该模型条件下,一维和二维反演都能得到反映模型真实信息的结果,REBOCC三维反演更偏向重建原始模型的宏观电性特征;在一维反演结果中,对于深部的电性结构特征,通常TE极化模式的反演结果好于TM极化模式,TE/TM几何平均值反演结果介于前两者之间;在二维反演结果中,通常TM极化模式的反演结果好于TE极化模式,TE/TM联合模式反演结果与TM极化模型相当,甚至更好;模型REBOCC三维反演相比一维、二维反演更易受反演结果多解性影响,REBOCC三维反演结果偏向重建原始模型的宏观电性特征。     

城市小型河流水动力弥散和潜流交换过程

为研究城市小型河流中污染物的物理迁移过程规律,分析基流条件下流动水体与暂态存储区之间的滞留交互作用,采用溴化锂(LiBr)作为保守性示踪剂进行野外现场示踪试验,结合一维溶质运移存储模型(One-dimensional Transport with Inflow and Storage model, OTIS)定量解析潜流交换特性,估算纵向弥散系数(D)、潜流交换面积(A_s)、主河道断面面积(A)和潜流交换系数(α).模型度量指标D_aI值和均方根误差值结果表征参数模拟结果可靠性高,拟合效果理想.由泵入点O至下游1 300 m设置的A、B、C、D 4处监测点的模拟结果表明,水文参数D、A_s、A和α均随水文条件而变,OB河段(0~600 m)潜流交换能力较弱,主要以对流弥散过程为主;BD河段(600~1 300 m)具有较强的暂态存储能力,对溶质的滞留时间长;BC(600~1 000 m)和CD(1 000~1 300 m)河段交换系数分别为(3.42×10-6±0.65×10-6)s-1和(2.87×10-6±0.81×10-6 )s-1;河段BC存在2.2×10-5m3/(s·m)的侧向补给流量.4个河段对比发现,城市河流渠道化、河床沉积物贫瘠等特征导致潜流交换能力弱化.     

一维热扩散湖模式在太湖的应用研究

利用在太湖获得的2010年8月11-28日的观测资料研究了一维热扩散湖模式在太湖的适用性,通过对比模拟进一步研究了影响太湖湖表温度模拟的主要因子。该湖模式对太湖最初的模拟结果并不理想,模拟的湖表面温度与观测有较大的系统性偏差,温度的日变化幅度与观测相比也偏小。通过分析该模式对太湖的模拟效果不理想的可能原因,针对太湖的生态环境和污染情况,设计了18个测试参数的敏感性试验,从敏感性试验的结果分析得到,适用于太湖的、依赖于湖泊类型的3个参数应做如下修改:消光系数（η₀）应放大3倍,湖泊表层吸收的太阳辐射系数（β）应取0.8,粗糙长度（z₀）采用公式计算得到。用新得到的适用于太湖的3个参数,模拟得到的结果与最初的模拟结果和观测资料对比,发现采用新的参数后,模拟结果比最初的模拟结果有了很大的改进,模拟的湖表温度基本接近观测,模拟的湖水垂直剖面时间序列图也跟观测吻合得较好,随之的感热、潜热通量的模拟也都与观测接近。最后,对输入湖模式的主要大气参量（太阳辐射、2 m气温及风场）±10%的误差引起的模式模拟的湖表面温度误差进行分析,结果表明该湖模式对大气强迫场的误差敏感度不高;相比之下,模拟结果对风场敏感性最小,对辐射和气温的敏感度相当。