格网DEM侵蚀学坡长并行计算方法

针对现有侵蚀学坡长串行算法在处理大区域海量数据时计算能力的不足,基于消息传递接口(MPI)并行化工具,提出了一种格网DEM的侵蚀学坡长并行计算方法,采用缓冲区更新计算策略,解决了并行计算过程中的数据依赖性问题。选取陕北黄土高原的两组不同分辨率的DEM数据对并行算法进行了测试,结果表明,提出的并行计算方法可以有效降低侵蚀学坡长的计算时间,并取得了较好的并行效率。     

煤矿排水对枣庄市清凉泉水源地污染作用分析

近年来,随着枣庄市陶枣矿区煤矿资源的大规模开发,矿坑排水对矿区地下水资源及生态环境产生的负面影响愈发突出,尤其是矿坑水沿地表河流排放途中,通过灰岩“天窗”大量灌入岩溶地下水中,造成地处煤田下游的清凉泉水源地水质遭受明显污染,供水水质难以保证。根据矿坑排水量、水源地水质、开采现状等分析认为,对串层污染进行封堵和控制矿坑水排放渠道是防治水源地污染的有效措施。     

伊犁蒙其古尔铀矿床含矿层砂岩中黄铁矿形成机制及对铀成矿的指示意义

蒙其古尔铀矿床为伊犁盆地南缘大型层间氧化带砂岩型铀矿床,为查明该矿床含矿层中黄铁矿成因及其形成机制,探讨微生物参与铀成矿过程。文章对含矿层砂岩中黄铁矿与铀矿物矿物学特征、黄铁矿S同位素与碳酸盐胶结物的C-O同位素开展细致研究。研究表明:①蒙其古尔铀矿床中铀主要以铀矿物与吸附铀形式存在,吸附铀主要为有机质吸附铀,铀矿物以沥青铀矿为主,多与黄铁矿、炭屑共生;②蒙其古尔铀矿床含矿层砂岩中黄铁矿主要以自形晶、草莓状和不规则状集合体产出,多与沥青铀矿、碳酸盐胶结物共生,其中黄铁矿S同位素(δ~(34)S_(V-CDT)=-68.4‰～22.1‰)与碳酸盐胶结物的C-O同位素(δ~(13)C_(V-PDB)=-10.2‰～-7.4‰,δ~(18)O_(V-PDB)=-9.6‰～-5.8‰)分析表明黄铁矿具有细菌硫酸盐还原(BSR)与有机物热解2种成因,并探讨了这2种不同成因黄铁矿的形成机制。③结合前人研究成果,认为硫酸盐还原菌(SRB)参与蒙其古尔铀矿床铀成矿过程,以间接还原方式为主,在有机质、黏土矿物与颗粒表面吸附U(Ⅵ)的基础上,通过硫酸盐还原菌(SRB)还原SO_4~(2-)产生的H_2S将U(Ⅵ)被还原成U(Ⅳ),形成铀矿物。     

吐哈盆地巴喀气田八道湾组致密砂岩储层分析及孔隙度演化定量模拟

在吐哈盆地巴喀气田八道湾组砂岩储层特征及控制因素分析的基础上,以地史时间和埋藏深度为变量,以现今孔隙度为约束条件,紧密结合埋藏阶段分析和关键成岩阶段窗口分析,分阶段建立了巴喀气田八道湾组储层从埋藏初始直至现今的孔隙度演化定量模型,动态体现了整个孔隙度演化过程.定量模拟结果表明:研究区储层总孔隙度演化是一个七段式分段函数,各个阶段孔隙度演化或遵循减小模型或遵循减小与增大的叠加模型;机械压实阶段的孔隙度减小模型是以埋深为自变量的函数,压实和胶结综合作用阶段则是以埋深和埋藏时间为变量的函数;孔隙增大模型适用于R.值范围为0.6％～1.0％的溶蚀作用窗口内;早期的纯机械压实作用使八道湾组储层已接近致密,溶蚀作用之前的压实和胶结综合作用导致储层致密化,而后期溶蚀作用仅在一定程度上提高了孔隙空间.     

层状岩质边坡遍布节理模型的三维稳定性分析

采用遍布节理模型（ubiquitous-joint模型）描述层状岩体的各向异性特征,并探讨ubiquitous-joint准则中安全系数的计算方法,利用FLAC3D分析层理倾角、倾向与边坡稳定性之间的关系,结果表明：（1）对于顺层边坡,当岩层倾向与坡面倾向的夹角 较小时,边坡的安全系数随层理倾角 的增大呈现先减小后增大再减小的趋势; 20°～30°时,安全系数得到最小值, 60°时安全系数得到最大值;不同破坏型式导致安全系数变化规律之间的差异;当 较大时,曲线从先减小后增大的趋势转变为不断增大的趋势;根据 的大小将顺层边坡分为两类：当 45°时,为严格顺层边坡;当45° 90°时,为斜向层状边坡。（2）对于逆层边坡,当 45°时,曲线呈现先增大再减小然后又增大的趋势,各曲线随 的变化程度基本相同;但减小过程中,随着 的增大,曲线的斜率逐渐减小,边坡稳定性的各向异性程度减小;当 45°时,曲线随 的增大呈现非线性单调递增的趋势。（3）安全系数与 之间的关系表明,对于顺层坡,随着 的增大,边坡安全系数逐渐增大;对于逆层坡,当 较小时( 150°),边坡安全系数随 的增大而逐渐减小;随着 的增大( 120°),边坡安全系数与 呈现递增关系。     

低剂量对比剂联合低监测阈值在头颈部CTA中的应用

目的：探讨Philips Brilliance iCT低剂量对比剂联合低监测阈值在头颈部CTA中的应用。资料与方法：选取我院需行头颈部CTA的患者60例,随机分为A组（观察组）、B组（实验组）。分别行不同的注射扫描方案：A组对比剂60 mL,生理盐水40 m L,监测阈值150 HU;B组对比剂40 mL,生理盐水40 mL,监测阈值90 HU;两组注射速率均以5.0 mL/s,均监测主动脉弓。分别对两组图像进行主观评价（目测法）和客观评价（分别测量上腔静脉、主动脉弓、左右颈总动脉、基底动脉及左右大脑中动脉的CT值）来进行图像质量分析。结果：①主观评价A、B两组图像均能满足诊断要求,B组明显减少上腔静脉对比剂残留;②A、B两组图像质量等级评分无统计学意义（P> 0.05）;③A、B两组各部位测量的CT值比较均有统计学意义（P <0.05）。结论：Philips Brilliance iCT在头颈部CTA检查中采用40 m L对比剂配合主动脉弓90 HU的监测阈值的方法,虽然图像CT值比较有差异,但图像质量完全能够满足临床诊断要求。且在保证图像质量的前提下,降低了对比剂用量,减少了上腔静脉对比剂的残留,降低了过敏反应的发生及对比剂对肾功能的损伤。     

多时相遥感影像提取湖泊边界信息的融合算法

为了有效改善遥感影像提取湖泊边界信息的可靠性和精度,减少人为误差,提出了一种利用多时相遥感影像提取边界信息的加权平均融合算法以及误差的域法修正处理方法.结果表明,该方法能有效融合各时相影像信息,提高湖泊边界信息提取的可靠性,并且对融合时相变化较大的湖泊边界都有一定的普适性.通过融合算法提取的呼伦湖面积为1928.35km2,修正后的面积为1929.85km2.通过利用地统计学理论对算法的验证及误差分析,得出相对误差空间变异拟合模型的块金方差与基台值之比都小于25%,具有很强的空间相关性,修正后的数据空间相关性要优于融合数据,空间变程也得到了有效地降低,显示修正后的数据对半方差函数理论模型的拟合程度更好.     

淮河流域夏季极端降水频次空间分布的客观分类及其形成机理

本文基于1961～2016年淮河流域四省（江苏、安徽、河南、山东）逐日降水观测资料及全球大气再分析资料,利用K均值聚类、旋转经验正交函数分解对淮河流域夏季极端降水频次分布进行了客观分类,利用统计诊断和数值模拟的手段讨论了其相关环流异常和形成机理。结果表明:（1）淮河流域夏季极端降水频次的空间分布可客观分为以极端降水主要发生在淮河流域33°N以南地区的南部型,主要发生在32°～36°N之间的中部型,和主要发生在34°N以北的北部型这三种分布类型;（2）南部型极端降水频次分布与西北太平洋副热带高压异常偏西偏南有关,西北太平洋异常反气旋北侧的异常气旋性环流使得水汽输送停留在淮河流域南部,导致南部极端降水频次偏多。中部型对应淮河流域受鞍型场环流结构控制,导致中部极端降水频次偏多。北部型极端降水频次分布时,淮河流域处于反气旋性环流异常西南侧,偏南风将水汽输送至淮河流域北部,导致北部极端降水频次偏多;（3）南部型和北部型的极端降水频次分布相关环流异常分别受厄尔尼诺和拉尼娜相关海表温度异常所影响,而中部型极端降水频次分布的相关环流异常是巴伦支海/喀拉海海冰异常在欧亚大陆上空激发的自西北向东南传播的准定常罗斯贝波所导致的。     

录井校正新方法在东海盆地西湖凹陷的应用

录井检测受井场环境影响较大,主要包括钻井环境与检测环境,钻井环境包括井径、钻时、钻井液排量、是否取心等,检测环境主要是井场自然环境与检测器内部环境。一般情况下检测器内部环境基本一致,而井场自然环境在一定时间内变化不大,故而钻井环境对录井检测的影响较大,钻井环境校正是进行录井解释前的重要工作。常规Tg(全烃含量)校正方法在取心、低异常、多互层等层段使用效果较差,经过分析,将校正方法在取心、钻时、排量、基值方面进行优化,实践验证,新校正方法效果较好,适应了西湖凹陷的地质研究需求。     

Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation, and is called groundwater contaminant source identifcation. This is a kind of typical groundwater inverse problem, and the solution is usually ill-posed. Especially considering the spatial variability of hydraulic conductivity field, the identification process is more challenging. In this paper, the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model (MT3DMS) and a data assimilation method (Iterative local update ensemble smoother, ILUES). In addition, Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction. In practical problems, the geostatistical method is usually used to characterize the hydraulic conductivity feld, and only the contaminant source information is inversely calculated in the identifcation process. In this study, the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field. The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site, and ILUES has good performance in solving high-dimensional parameter inversion problems.