P—Ⅲ型分布Φ值数值计算方法比较

皮尔逊Ⅲ型频率曲线的数值计算方法是水文频率计算中的难点之一。本文对近年来提出的几种数值计算方法进行了分析比较，并推导出了各数值计算方法试算离均系数Φ值的迭代公式。最后推荐计算精度高、运算速度快、收敛范围广的级数展开法迭代公式。     

盾构隧道下穿管道施工引起的管-土相互作用研究

采用Loganathan公式研究了盾构隧道下穿管道施工引起的地下管道处土体竖向位移,利用考虑土中剪力传递的Pasternak模型模拟管-土相互作用,运用修正Vlasov模型中的迭代流程计算出Pasternak模型的关键参数——弹性系数k与剪切系数gs。将计算结果与已有文献结果及工程监测数据进行对比,深入分析了迭代求解k、gs值的Pasternak模型与传统模型的计算差异,并进一步研究了土中剪力、管道与隧道的夹角、土体弹性模量及隧道半径的变化对管-土相互作用的影响。研究结果表明：迭代求解的k、gs值能提升Pasternak模型的精确度;土中剪力对管道竖向位移计算值的影响可达15.3%;随着管道与隧道夹角的减小,管道的竖向位移增大、弯矩减小;土体弹性模量与隧道半径的增大均会增加管道的竖向位移和弯矩。     

基于高性能有限单元的单桩稳定性分析法

现行的桩基设计方法主要基于线弹性理论及采用半经验假定,难以准确地检验长桩在土体非线性条件下的稳定性。基于非线性有限单元分析理论,提出了高性能桩单元分析法用于非线性分析,可直接检验单桩稳定性且无需假定桩的计算长度系数。在桩单元推导过程中,通过整合在单元内部的连续弹簧以考虑土-结构相互作用(SSI),能够大幅提升计算效率。使用牛顿-拉夫逊迭代法进行迭代运算,利用推导的相应单元切线刚度矩阵预测位移,并通过割线关系减少每一步迭代中产生的误差,在桩的大变形条件下采用更新拉格朗日法确定平衡条件。算例验证表明,桩单元模型在考虑土体非线性条件下,能高效、可靠地对单桩进行分析和设计。     

考虑桩-岩接触面流变的桩岩联合受力性能研究

桩-岩接触面的性质对桩岩联合受力的状况有重要影响。大量的工程实践表明,桩-岩接触面存在显著的剪切流变现象。在前人研究的基础上,介绍了以位移协调法为基础的双曲传递函数迭代计算方法,并据此建立了一种可考虑桩-岩接触流变的桩岩联合受力迭代计算方法。与现场模型桩试验实测结果对比分析表明,双曲传递函数迭代计算方法能够较好地描述桩岩联合受力性能。当考虑了桩-岩接触面的剪切流变后,桩身各处轴力随时间的增长逐渐增加,直至趋于稳定,桩的轴力比不考虑流变时桩的轴力大。迭代计算结果与实测结果吻合较好,对桩岩联合受力性能的研究有一定参考价值。     

大型地下洞室群围岩应力-损伤-渗流耦合分析

以渗透体力来考虑渗流场的力学效应,建立了应力-损伤-渗透系数关系方程来考虑应力和损伤对渗流场的影响,结合岩体结构的三维弹塑性损伤有限元分析,建立了大型地下洞室开挖围岩应力-损伤-渗流耦合的计算模型。该计算模型求解的难度主要体现在岩体材料弹塑性、损伤、渗流自由面边界、渗流溢出边界、应力-损伤-渗流相互影响关系等。提出分步迭代法对以上因素进行归纳后按一定顺序分别进行迭代求解,取得了良好的计算效果。将该方法应用于某水电站大型地下洞室群分析,得出了一系列有意义的结论。     

-种确定土体固结系数C_V的新方法

计算固结系数的方法很多, 它们各有千秋, 但无论哪种都不能计算土体在任意固结时刻的固结系数。本文把几种计算方法的优点结合起来, 找到了-种计算任意固结时刻固结系数的新方法, 并用室内系列固结实验数据进行了计算。现场观测资料分析计算结果与之对比结果表明, 该方法行之有效。     

应用遗传算法对镭氡平衡曲线的研究

本文通过分析常规计算镭氡平衡系数的误差因素及钻孔中镭氡平衡的动力学过程,建立了求解镭氡平衡系数的数学模型.应用遗传算法来设计、求解相应模型的待定参数,编写相关程序段,进行迭代求解.计算出了单矿段及铀矿床总镭氡平衡系数.     

黏弹-双曲线Drucker-Prager塑性模型应力更新隐式算法

黏弹-塑性模型是由黏弹性模型与塑性元件串联而成,可视为黏弹-黏塑性模型中黏塑性黏度参数趋于0的一种极限情况,为黏弹性材料的结构破坏分析提供了一个途径。把黏弹-塑性模型的应变增量分解为黏弹和塑性增量两部分,考虑黏弹性应变历史,把黏弹性积分型本构关系在一个时步内线性化,定义与时间增量相关的剪切模量和体积模量,导出应力递推公式,把黏弹-塑性本构积分转化为与弹塑性相似的形式。针对由黏弹性和双曲线Drucker-Prager塑性、各向同性硬化的黏弹-塑性模型,通过黏弹性预估和塑性校正"二步"算法实现对应力的更新,给出完全隐式算法和最终的计算公式。算例比较分析表明,由于迭代过程中仅需要简单的函数计算,该算法具有很好的收敛性。一般经过2次迭代运算后,屈服函数值已达到10-10的量级,应力点便返回到屈服面上。     

TEM中心回线法计算考虑关断时间的全区视电阻率

详细推导了TEM法中心回线装置发射电流为斜阶跃波形条件下全区视电阻率迭代反演计算方法.对正演过程中对关断时间的积分采用了解析表示.对反演过程中求导运算,通过变量变换,将对电阻率的求导变换为对时间的求导,从而也给出简单的解析关系.避开了数值积分运算和差分替代导数的近似,提高了计算精度和速度.为克服迭代反演过程中要多次计算概率积分的耗时问题,采用标准数据库的方法一劳永逸地加以解决,也大大提高了计算精度和速度.对实测资料进行了计算,并通过算例对一些问题进行了讨论.     

乌库尔齐铀矿床镭-氡放射性平衡系数计算及讨论

从分析镭-氡放射性平衡破坏的原因出发,以新疆乌库尔齐可地浸砂岩型铀矿床为例,探讨了镭-氡放射性平衡系数的计算方法,并对不同的计算方法进行了讨论。该系数的取得,为乌库尔齐铀矿床测井数据的修正、矿体边界的划分、铀矿资源量的计算提供了科学依据。     

The significance of unusual zoning in olivines from FAMOUS area basalt 527-1-1

Major element, Ni, Mn and Ca electron microprobe analyses of olivine phenocrysts in one of the most primitive basalts from the FAMOUS area of the Mid-Atlantic Ridge, 527-1-1, reveal two different olivine populations, distinguished by their zoning characteristics. The often skeletal Group I olivines have zoning profiles with high forsterite, high Ni, low Mn, and low Ca cores. These profiles can be explained by low pressure crystallization from the 527-1-1 magma on cooling. The equant, often megacrystic Group II olivines also have high forsterite, low Mn cores, but the cores have low Ni compared to rims. Thus they are normally zoned with respect to Mg, Fe, and Mn but reversely zoned with respect to Ni. For example, Ni ranges from 1,700 ppm at Fo_90.5 in the core to 2,100 ppm at Fo_89.5 at the rim. In view of the published whole rock and mineral data of le Roex et al. (1981), the most likely explanation for these data is that the Group II olivines are xenocrysts assimilated from solidified plagioclase-pyroxene basalts through which the 527-1-1 basalt ascended. The diffusion rate of Ni and size of the xenocrystic olivines are used to calculate the residence time of the assimilated olivines in the magma. An alternative hypothesis would be a high pressure origin for the Group II olivines. This would be possible if future experiments show that with increasing pressure the partition coefficient for Ni decreases relative to the partition coefficient for Mg for a given bulk composition. Current evidence suggests this is unlikely. The data from 527-1-1 and other samples from the FAMOUS area require magmas with distinct Ni-MgO-FeO characteristics. In general, MORB from different ridge segments fall on distinct trends on plots of MgOvs. Ni as well as MgOvs. FeO. Calculation of the MgO and Ni contents of primary magmas suggests a mechanism by which such distinct trends could come about.     

既有房屋饱和土地基极限承载力增长的探讨

饱和土地基在既有房屋长期荷载作用下,地基土强度参数值随着地基土的排水固结和压密而提高,地基极限承载力也同时增大,如何估算此提高值而挖掘其潜力,传统方法只是借助于统计回归讨论了设计承载力。本文从土力学基本原理出发,借助Prandtl假定,对既有房屋饱和土地基极限承载力的增长进行了探讨。     

Preparation and Certification of High-Grade Gold Ore Reference Materials (GAu 19-22)

Two types of gold deposit with both good homogeneity and a high-grade of gold were selected to prepare four gold ore reference materials (GAu-19,GAu-20,GAu-21 and GAu-22) by the Institute of Geophysical and Geochemical Exploration (IGGE),China.Seven laboratories participated in the certification programme. Volumetric methods for GAu19-21 and atomic absorption spectrometry for GAu-22 were used for the homogeneity testing,the coefficient of variation being found to be less than 3%.After sample digestion and preconcentration,the samples were analysed by flame atomic absorption spectrometry (AAS),colorimetry,neutron activation analysis (NAA)and volumetric analysis. The certified values for the gold concentration in GAu19-22 are 18.3 μg g -1, 32.3 μg g -1, 53.0 μg g -1 and 5.72 μg g -1,respectively.     

Spectrophotometric Determination of Vanadium in Carbonaceous Shales (Stone Coal Ores) Using 2-(5-bromo-2-pyridylazo)-5-diethylaminophenol

A selective and sensitive spectrophotometric method has been developed for determining the total amount of vanadium in carbonaceous shales (stone coal ores). The method is based on the reaction of vanadium(V) with the chromophore reagent 2-(5-bromo-2-pyridylazo)-5-diethylaminophenol (5-Br-PADAP) in the presence of hydrogen peroxide. In a 0.072 mol l⁻¹ sulfuric acid medium, 5-Br-PADAP reacts with vanadium(V) to form a red-violet complex with a maximum absorption peak at 596 nm with an apparent molar absorption coefficient of the complex of 8.45 × 10⁴ l mol⁻¹ cm⁻¹. Beer's law was obeyed in the range 0–25 μg vanadium in 25 ml of solution, with a correlation coefficient of 0.9995. Interferences due to various non-target ions were also investigated and high quantities of other common inorganic ions were tolerable. The method involved the dissolution of the ore sample by Na₂O₂ fusion, followed by filtering of the alkali solution after which Fe(III), Cu(II), Ni(II) and Co(II) etc. were effectively separated from the solution by precipitation in a NaOH solution. Selectivity was increased with the use of EDTA as a masking agent. The vanadium in ore samples was determined with a relative deviation (RSD) between 0.20 and 0.76%, and has been successfully applied to the determination of vanadium-bearing stone coal ores. The results indicated that the accuracy of 5-Br-PADAP spectrophotometry is comparable with the ICP-AES method. The characteristics of the method, i.e., simplicity, selectivity, sensitivity and rapid calibration, make it specially suitable for routine analysis.     

Calculation of the contribution of single-zone production by ultraviolet spectrum technique： A case study in well QHD32-6-3 field

Based on the established mathematic model and graphic interpretation, a new method, which is used to calculate the contribution of single-zone production in a commingled producing well by the ultraviolet spectrum technique, has been established. The standard plate was drawn using the extinction coefficient E of sample oils formulated artificially as y-axis and the wavelength as x-axis. The curve resulting from the UV analysis of sample oils in the commingled well was inserted into the standard plate and compared. The proportion of each single zone in the commingled producing well was identical with the proportion of the curve which is closest to the curve of sample oils formulated artificially. In the well QHD32-6-3 field, taking well A22 for example and using this method, the production contribution of a single zone was calculated. The result showed that the Nm4 zone is a major "contributor", the proportion of the Nm4 zone is 70%, and that of the Nm1 zone is 30%. The ultraviolet spectrum technique provided a new reservoir geochemical technique of monitoring production contribution, especially for biodegraded heavy oil, but it has some limitation, just depending on the GC fingerprint technique.     

1965-2015年新疆夏季不同等级降水的空间分布特征

根据新疆51个台站1965-2015年夏季逐日降水资料,将降水划分为小雨、中雨及大雨3个等级,分析了新疆近51 a夏季不同等级降水量、降水日数及降水强度的空间分布特征,并讨论了各等级降水日、降水量及降水强度与总降水量的空间相似程度以及各等级降水对夏季总降水的贡献。结果表明：新疆降水主要集中在夏季,并以小雨为主。以天山山脉为界,南北两疆降水空间分布存在明显差异,北疆夏季降水量（日）占年降水量（日）的36%~45%（36%~39%）,南疆夏季降水量（日）占年降水量（日）的51%~63%（48%~60%）;新疆夏季不同等级降水量、降水日及降水强度的空间分布不均匀。新疆夏季总降水量与各等级降水量的空间相似系数最为密切,与各等级降水强度的空间相似系数相对较小;新疆夏季小雨贡献率最大,中雨其次,大雨最小,夏季降水量和降水日的变化主要受小雨的影响。     

模糊模式判别在油气预测中的应用--以泌阳凹陷前杜楼地区北部断块区98-2井区为例

通过对泌阳凹陷前杜楼地区北部断块区98-2井区资料的分析，运用模糊模式判别的原理，选取主振幅、主频率、有效吸收系数、砂岩厚度和层间速度差等参数来反映信息。建立了模糊因素集、模糊评判集、隶属函数和权重向量。利用所建立的模糊模式，对98-2区的7个砂层进行了模糊判别，得到98-2井区位Ⅰ级异常，即最有利含油气区。     

垃圾土的渗透特性试验

用取自京郊某垃圾堆填场的垃圾土样进行渗透试验,比较不同密度情况下渗透系数的差别。试验中同时记录了渗透系数随时间的变化情况。试验结果表明,垃圾土的渗透系数随密度的不同,可以在(10-3~10-7)cm/s之间变动。试验中也发现垃圾土的渗透系数随时间而逐渐减小,大约6天左右趋于稳定。对垃圾土的渗透系数室内试验测定应以观测7天时间为宜。     

1956-2012年黄河源区流量演变的新特征及其成因

利用黄河源区唐乃亥水文站1956-2012年的逐月流量实测资料,通过趋势分析、小波分析以及不均匀系数和贡献的计算,分析了1956-2012年57 a来黄河源区流量的趋势、年内双峰型变化新特征及其主要成因.结果表明：近57 a来黄河源区流量演变包括2个上升时段（1956-1967年和2002-2012年）、1个下降时段（1989-2002年）和1个振荡调整时段（1967-1989年）.流量年内分配不均匀性最大的是1980年代,最小的是1990年代.流量年内分配主要表现为双峰型,峰值点主要出现在7月和9月,但出现月份存在年际、年代际变化,且峰值点的月份也随着丰、枯水年存在着变化.降水对流量的正贡献呈7、9月双峰型,温度的负贡献则使9月峰值消失.当夏季高原附近区域,尤其是高原北侧低层600 hPa位势高度场降低,高原夏季风偏强时,有利于黄河源区降水增加、流量增大,出现流量年内分配的峰值点,且夏季高原高度场的变化对流量年内分配第二个峰值点（秋季）的出现有一定的指示意义.     

城镇化发展背景下那曲地区气温变化的特征分析

根据西藏那曲气象站（NQ站）1955-2014年和其相邻野外站（BJ 站）2001-2014年的气象观测资料,在那曲县城镇化发展背景下,对比分析了NQ站和BJ站气温变化的差异,从自然要素和人类活动要素两方面探讨了两站气温变化差异的原因。结果表明：近60a NQ站气温呈上升趋势,且2001-2014 年的气温上升率大于1955-2014年和1987-2000年的气温变化趋势;NQ站2001-2014年的平均气温、最高和最低气温的升温率,分别为0.52℃·（10a）^-1,0.60℃·（10a）^-1和0.67℃·（10a）^-1;同时期BJ站的平均气温、最高和最低气温升温率,分别为0.76℃·（10a）^-1、2.72℃·（10a）^-1和0.32℃·（10a）^-1。NQ站和BJ 站的最高气温变化均主要与自然要素有关,NQ站的最高气温与日照（0.67）、降水（-0.62）、水汽压（-0.58）相关系数较高,BJ站的最高气温与日照（0.67）、NDVI（-0.63）、水汽压（-0.57）相关系数较高;NQ站的最低气温与人均GDP（0.53）、GDP（0.50）、总人口（0.47）相关系数较高,说明NQ站的最低气温变化主要是受人类活动因素的影响,人类活动对NQ站最低气温的热岛贡献率为52.2%;BJ站的最低气温与各因素相关系数均较小且接近,是在气候背景下的自然变化。