南极东方站深冰层及冰下湖钻探技术

自1970年至今,前苏联和俄罗斯在南极东方站持续进行了近50年的冰层钻探活动,先后攻克了包含粒雪层、冰层、冰岩夹层和湖水冻结冰的复杂冰层钻进难题,逐渐形成了一套集热融取芯钻探、电动机械取芯钻探和分支孔钻探等为一体的深冰芯钻探技术。创造了冰层最深干孔钻进深度记录（952.4 m）、最深热融取芯钻进记录（2755 m）、最深冰芯钻探记录（3769.3 m）,累计进尺达13000 m,并获取了总长超46 m的含湖水冻结冰样品的冰芯。东方站的钻探活动对极地冰层钻探技术的发展起到了巨大的推动和引领作用,同时积累了宝贵的深冰钻探经验。通过对东方站深冰钻探技术的系统梳理,将为我国正在实施的深冰芯钻探和即将开启的冰下湖科学钻探提供重要的借鉴。     

寒区渠基黏土热参数最优概率分布

为分析寒区渠基黏土热参数的随机分布特征及概率分布模型,以寒区渠基黏土的导热系数为样本,结合经典分布拟合法、多项式逼近法、最大熵法和正态信息扩散法,分别对寒区渠基黏土热参数的概率分布规律进行了研究。首先通过分析热参数的离散性,并比较概率分布曲线、拟合检验值和累计概率分布值,对不同方法描述热参数随机性的优劣进行了评价;然后,基于寒区渠基黏土热学参数对温度的敏感性,提出了一个可以达到理想拟合精度的寒区渠基黏土热参数概率推断的区间取值标准。研究结果表明：寒区渠基黏土的热参数具有随机变量的特征;正态信息扩散法可以描述热参数样本的随机波动性;在4种方法中,正态信息扩散法的拟合精度最高。使用3.5σ法,将[μ-3.5σ,μ+3.5σ]（μ为随机变量的均值,σ为标准差）作为概率函数推断时的取值区间,同时考虑偏度的影响,可使得累计概率值达到1.000 0的精度,能够较准确地推断热参数的概率分布函数。     

活动断裂带岩石的变形机制和摩擦性质—以圣·安德列斯断层为例

岩石的变形机制和力学性质对认识活动断裂带的孕震机制和发震过程都具有重要意义。圣·安德列斯（San Andreas）转换断层是太平洋板块与北美板块的构造边界,也是重要的地震活动带,总结了美国南加州地震中心对圣·安德列斯断层钻孔断层岩的研究进展。对断层泥的高速摩擦实验发现：动力弱化现象只在高速滑移时才发生,伴随着摩擦生热,岩石的摩擦系数显著降低。断层泥中富镁黏土矿物的富集可以显著降低断层的摩擦系数,促进蠕滑,黏土矿物的不均一分布可导致断裂带变形行为和地震分布的突变。与断层泥和碎裂岩不同,粉碎岩形成于快速传播地震的断裂端部,具高孔隙率和高渗透率,流体活动性强。此外,断层岩的有机质成熟度可用来追溯古地震造成的热异常。在川滇中国地震科学实验场开展科学钻探项目,研究断层岩样品的显微结构、物理性质和变形历史,以期为理解大陆汇聚与转换挤压过程中复杂的地震活动提供重要信息。     

热红外遥感技术在江西寻乌地区地热勘查中的应用

热红外遥感是地热资源勘查中的常用技术,对江西寻乌地区的Landsat 8热红外影像采用单窗算法,根据MODIS影像计算大气水汽含量,并逐步获取反演参数,得到地表温度反演结果;利用MOD11_L2地表温度产品对结果进行对比分析,验证反演有效性。通过对已有地热资料和地质构造背景的研究,总结地热形成的有利构造条件,分析地表反演温度与已有地热及断裂构造之间的空间关系,建立地热异常区圈定依据,并剔除假异常,最终在研究区圈定了11个潜在地热异常区,为该区地热勘查提供依据。     

An Operational Precipitable Water Vapor Retrieval Algorithm for Fengyun-2F/VISSR Using a Modified Three-Band Physical Split-Window Method

The Visible and Infrared Spin-Scan Radiometer(VISSR) onboard the Fengyun-2(FY-2) satellite can provide valuable thermal infrared observations to help create a precipitable water vapor(PWV) product with high spatial and temporal resolutions. The current FY-2/VISSR PWV product in operation is produced by using a traditional two-band physical split-window(PSW) method, which produces low quality results under dry atmospheric conditions. Based on the sensitivity characteristics of FY-2 F/VISSR water vapor channel and two split-window channels to atmospheric water vapor, this study developed a new, robust operational PWV retrieval algorithm for FY-2 F to improve the operational precision of the current PWV product. The algorithm uses a modified three-band PSW method, which adds a scale for the water vapor channel in the improved three-band PSW method. Integrated PWV products from the radiosonde data in 2016 are used here to validate the precision of the PWV retrieved by the modified three-band and traditional two-band PSW methods. The mean bias, root mean square error(RMSE), and correlation coefficient of the PWV retrieved by the modified three-band PSW method are 0.28 mm, 4.53 mm, and 0.969, respectively. The accuracy is much better than the PWV retrieved by the two-band method, whose mean bias, RMSE, and correlation coefficient are 12.67 mm, 29.35 mm, and 0.23. Especially, in mid-or high-latitude regions, the RMSE of the PWV is improved from 10 to 2 mm by changing the inversion in the two-band method to the modified three-band PSW method. Furthermore, the modified three-band PSW results show a better consistency with the radiosonde PWV at any zonal belt and season than the two-band PSW results. This new algorithm could significantly improve the quality of the current FY-2 F/VISSR PWV product, especially at sites where the actual PWV are lower than 15 mm.     

夏季亚洲高空急流纬向非对称变异与北大西洋海温和欧亚陆面热力异常的可能联系

基于1979~2015年ERA-Interim再分析资料,分析了夏季亚洲高空急流纬向非对称变异特征及其可能的外强迫因子。研究发现夏季亚洲200 hPa纬向风异常EOF第二模态（方差贡献为16.4%）主要表现出了急流纬向非对称的空间异常形态,反映了西亚和东亚区域急流南北偏移的反位相变化。通过进一步的诊断分析,我们发现急流纬向非对称变异与北大西洋海表温度（简称海温）和欧亚陆面热力异常可能存在一定的联系。北大西洋三极型海温异常会激发出向下游传播的异常波列,夏季该波列在欧亚大陆上空的异常环流中心与急流纬向非对称相关的异常环流中心对应一致,其中东欧平原的异常反气旋和巴尔喀什湖附近的异常气旋对西亚急流变化存在影响,东亚地区急流的变化与贝加尔湖北部异常气旋和贝加尔湖南部的异常反气旋有关。对比欧亚土壤湿度关键区内垂直环流,陆面热力异常可能会改变局地环流进而影响急流变异,且这种影响存在区域差异。     

由非季风区向季风区过渡过程中大气边界层结构的变化分析

利用中国西北中部具有代表性的非季风区、夏季风影响过渡区和季风区的7个高空站的2013年夏季晴天07时、13时、19时（北京时）的大气边界层资料,通过分析大气边界层位温、比湿、风速的垂直结构,发现大气边界层结构及厚度在不同区域的分布特征:稳定边界层厚度、残余层顶高度和对流边界层厚度从非季风区、夏季风影响过渡区至季风区出现阶梯性大幅降低,从非季风区至夏季风影响过渡区,以及从夏季风影响过渡区至季风区,对流边界层厚度降幅依次为25.6%和81.8%,稳定边界层厚度降幅依次为58.3%和41.8%;在稳定边界层条件下,可观察到低空急流的存在,非季风区低空急流出现高度明显高于夏季风影响过渡区和季风区,且非季风区的低空急流风速也明显大于夏季风影响过渡区和季风区。通过分析与大气边界层发展最为密切的陆面热力因素在不同气候区的分布,净辐射值、日地-气温差最大值以及感热通量值在非季风区大于夏季风影响过渡区和季风区,从陆面热力过程为非季风区大气边界层厚度大于夏季风影响过渡区和季风区提供了理论依据。     

定量研究雅安地形坡向坡度对降水分布的影响

利用四川省雅安市30 m分辨率基础高程数据,提取栅格的坡向和坡度参数,将雅安307个区域自动站在2017年汛期(6—9月)共50次的降水天气个例,分为16次大尺度降水和34次中小尺度降水,使用对应时次的欧洲中心细网格0.25°×0.25°再分析风场资料,根据不同的站点地形高度将风场合成平均风场,和各站点地形的坡向和坡度计算出其动力抬升作用,同时使用当天日照和天文太阳辐射值来计算地形的热力抬升作用,与对应降水过程的降水分布进行多元线性回归,根据回归的标准系数的大小确定各自变量对降水分布的影响,得出以下结论:(1)中小尺度降水中,地形的热力抬升作用对降水分布的影响作用最大,其次是海拔高度,地形的动力抬升作用在三者中对降水分布的影响最小;(2)在大尺度降水中刚好相反,地形的动力抬升对降水的分布影响作用最大,海拔次之,热力抬升作用在三者中影响作用最小;(3)日降水量最大值的站点海拔高度基本位于1 000 m左右,与抬升凝结高度对应较好;(4)从长期的统计来看,地形的动力作用和地表的植被情况对降水分布的影响最大。在实际预报工作用,根据不同的降水类型,关注不同的动力和热力作用对于判断降水分布大值区的位置有较好的参考作用。     

A 3D discrete FEM iterative algorithm for solving the water pipe cooling problems of massive concrete structures

Water pipe cooling has been widely used for the temperature control and crack prevention of massive concrete structures such as high dams. Because both under‐cooling and over‐cooling may reduce the efficiency of crack prevention, or even lead to great harm to structures, we need an accurate and robust numerical tool for the prediction of cooling effect. Here, a 3D discrete FEM Iterative Algorithm is introduced, which can simulate the concrete temperature gradient near the pipes, as well as the water temperature rising along the pipes. On the basis of the heat balance between water and concrete, the whole temperature field of the problem can be computed exactly within a few iteration steps. Providing the pipe meshing tool for building the FE model, this algorithm can take account of the water pipe distribution, the variation of water flow, water temperature, and other factors, while the traditional equivalent algorithm based on semi‐theoretical solutions can only solve problems with constant water flow and water temperature. The validation and convergence are proved by comparing the simulated results and analytical solutions of two standard second‐stage cooling problems. Then, a practical concrete block with different cooling schemes is analyzed and the influences of cooling factors are investigated. In the end, detailed guidance for pipe system optimization is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.