超高压气藏天然气拟对比压缩系数的数值解法

快速准确地获取超高压气藏天然气藏流体压缩系数是一个备受关注的问题。塔里木油田M超高压气藏的测试实例证实了DAK、DPR及HY 3种方法计算偏差因子的可靠性,在此基础上使用一阶中心差分算法计算了该气藏的天然气压缩系数。与Trube图版的对比表明,DPR方法的计算结果具有最高的精度,推荐数值计算时采用。数值实验结果表明,Tpr(拟对比温度)大于1.2时计算精度较高且与差分步长无关,相对误差小于5%;Tpr小于1.2时不建议使用数值方法。在借鉴部分超高压天然气偏差因子研究基础上,拓展了Trube图版的使用范围,将拟对比压力Ppr应用范围从15扩大到30,Tpr应用范围从2.0扩大到3.0,并重绘了图版。研究结果对超高压气藏的相关研究具有一定的理论参考和实际应用意义。     

NAMs柯石英中结构水的红外光谱和第一性原理计算

从微观尺度研究结构水的分布状态可以为超高压变质岩的形成环境、构造演化动力学过程提供重要的依据.为探讨大别山地区超高压变质岩中"名义上无水矿物"（nominal anhydrous minerals,NAMs）结构水的分布特征、赋存状态与超微结构缺陷的关系,对大别山石马地区榴辉岩中的柯石英进行了傅立叶变换红外光谱（FTIR）分析和第一性原理计算.FTIR研究表明柯石英主要吸收峰为（Ⅰ）3 561~3 580 cm-1、（Ⅱ）3 433~3 462 cm-1和（Ⅲ）3 412~3 425 cm-1;柯石英颗粒结构水含量为15×10-6~52×10-6,平均值是32×10-6.第一性原理理论计算得到了柯石英（4H）_Si和（AlH）_Si复合缺陷超晶胞模型（2×1×1）的形成能分别是-4.92 eV和-3.10 eV;含氢缺陷模型计算结果得到3 526 cm-1和3 198 cm-1的拉曼峰与柯石英的合成实验结果基本符合.FTIR分析表明石马地区柯石英结构水含量具有不均一性;模拟计算得到（4H）_Si复合缺陷模型比（AlH）_Si有更低的复合缺陷形成能,有更加稳定的结构,柯石英结构水中（OH）₄$ \Leftrightarrow $Si氢结合机制是优先模式,为实验研究提供理论依据.      

沁水盆地柿庄北深部煤层水平井CO2注入参数研究

基于CO₂的吸附能力大于CH₄的原理,向深部煤层中注入并封存CO₂同时置换驱替CH₄增加煤层气井产量,具有环境和经济双重效益。研究深部煤储层特征,通过气、水两相多组分,单孔和双孔模型的三维储层数值模拟技术,模拟不同井底压力下,完成注入量所需要的时间。优化后的注入施工参数为井口注入压力低于10MPa,注入速率低于47 L/min。根据数值模拟结果选定注入压力,注入结果与模拟连续注入时间一致,说明该方法可以有效地模拟注入压力和注入量之间的关系,保障CO₂注入施工安全有效。      

基于梯度提升算法的岩性识别方法

传统的岩性识别方法如岩屑录井、钻井取心及测井资料解释等技术,对录井质量的依赖程度较高,识别精度与效率低,泛化能力差。随着计算机技术的迅速发展,将测井资料与计算机技术相结合开展岩性研究已成为岩性识别的有效手段。本文提出了一种基于梯度提升算法XGBoost和LightGBM的岩性识别方法。以苏里格气田苏东41-33区块下碳酸盐岩储层为例进行测试验证,采用该方法结合测井资料中的声波时差、自然伽马、光电吸收截面指数、密度、深侧向电阻率和补偿中子等6种参数进行岩性识别,并与KNN （K近邻分类器）、朴素贝叶斯和支持向量机等传统算法进行对比,结果表明,3种传统算法的岩性识别准确率分别为78.45%、74.43%和78.72%,基于梯度提升算法XGBoost和LightGBM的识别准确率分别达到了98.90%和98.72%,远高于传统算法。     

单多普勒雷达反演涡旋风场方法分析

由于单雷达在探测台风临近、登陆时仍具有双部或多部雷达难以具备的优势,目前利用单多普勒雷达反演涡旋风场的方法并不多见,本文对能够反演涡旋风场的3种方法NIVAP(自然坐标系的积分VAP方法)、EVAPTC(反演热带气旋的EVAP方法)、GBVTD(地基速度追踪法)进行了比较,分别从每种方法的假设条件和数学方法以及结果精度来比较3种方法对涡旋风场的适用性。利用中尺度气旋Rankine模式模拟了纯涡旋性气旋、只有环境风场时、环境风场和辐合(辐散)同时存在时的径向速度场及风场,比较了各种情况下的反演风场和模拟风场的相似系数。在对模拟风场比较之后,选取了2014年"威马逊"台风登陆前1h、登陆时及登陆后1h的3个实测雷达资料,对比分析了3种方法反演实测资料的风场特征,3种方法中EVAPTC方法最好,GBVTD方法在应用中有一定的限制,NIVAP方法较差。     

一次极端寒潮天气过程等熵位涡分析

利用2016年1月1日至31日的FNL资料,对一次极端寒潮天气过程进行了等熵位涡分析。结果表明:高位涡主体由极涡分裂而来,前面低位涡区的阻挡与后侧低位涡大气的北上加强了位涡的经向交换,高位涡空气不断由极地向南输送,使得高位涡主体不断加强维持。高位涡在由北向南移动的同时,也由对流层顶向下输送。此次寒潮过程主要有3股冷空气由上而下发展,位置均在高空急流轴的北侧,最南端的一股下沉气流最旺盛,这是其与高空急流相互作用的结果。强盛的冷空气下沉使得寒潮影响范围触及我国华南地区。随着高位涡的向南向下传输,一方面引起对流层中高层低涡系统迅速发展,当它移到中国东部地区时,东亚大槽迅速加深,使槽后强冷空气迅速向南爆发;另一方面,在高位涡输送的过程中,其后侧有强烈的下沉运动,使得地面冷高压快速发展,导致强寒潮天气的爆发。     

1961-2010年青藏高原气候变化特征分析

利用1961-2010年青藏高原及其周边地区158个气象站温度(包括平均温度、最低和最高温度)、降水和风速资料,对青藏高原的气候变化特征进行了分析。结果表明:(1) 1961-2010年青藏高原主体正在变暖变湿,但是高原东侧部分地区正在变暖变干,同时高原整体风速都在减小。(2)升温主要是夜间的最低温度贡献的。不同地区升温速率有差异,中部地区高于东部地区;平均温度和最高温度分别在1994年和1997年发生突变,突变后升温速率明显加快;三种温度都存在准8年周期震荡,其他短周期及更长周期震荡表现不一致。(3)降水量空间分布上表现为从东南向西北逐级减少,并且出现过多次突变,突变时间分别为1965年、1977年和1995年,突变前后降水的变化速率明显不同,降水存在准4年和准10年周期震荡。风速存在18～20年周期震荡。(4)青藏高原平均温度、最低温度及最高温度EOF分解的第一载荷向量均表现出全区一致的正值,中心区位于94°E-97°E一带,说明青藏高原腹地是平均温度、最低温度及最高温度变化最敏感的地区。(5)平均温度、最低温度及最高温度EOF分解的第二载荷向量大体表现出高原主体与东部以及北部边缘地带变化趋势相反,即高原主体升温(降温)时,东部及北部边缘地带是降温(升温)的。     

钙稀土氟碳酸盐矿物中新多型体的HREM研究

利用高分辨电子显微术（ＨＲＥＭ）研究了钙稀土氟碳酸盐矿物系列中氟碳铈矿（Ｂ）／直氟碳钙铈矿（Ｓ）的混层结构。发现并确定了Ｂ＿５Ｓ＿２（ＩＨ，３Ｒ）、Ｂ＿５Ｓ＿４（３Ｒ，６Ｒ）和Ｂ＿（１８）Ｓ＿（１０）（ＩＨ）规则混层矿物中５种新多型的晶体结构类型、晶胞参数以及理论晶体化学式等。ＨＲＥＭ研究揭示出，在Ｂ＿ｍＳ＿ｎ（ｍ＞ｎ）型规则混层矿物新多型体的结晶过程中，成矿地质环境是复杂的。阐明了这些多型体是由该系列矿物中的氟碳铈矿（Ｂ）和直氟碳钙铈矿（Ｓ）结构单元层以不同比例沿ｃ轴方向的有序堆垛构成。