水流损失和热补偿共同作用对增强型地热系统(EGS)产能影响的研究

基于青海共和盆地-3705m地热田实测数据,结合流固耦合传热理论并运用Comsol软件,建立了离散型裂隙岩体流体传热模型。考虑水流损失和热补偿共同作用,模拟得到了开采过程中上、下岩层(盖层和垫层)为绝热不渗透、传热不渗透、渗透传热时,储层(上、下岩层和压裂层)温度场的变化特征,分析了产出流量、水流损失、产出温度、产热速率的变化规律。研究结果表明:采热过程中产出流量始终小于注入流量;产出流量增幅速率先增大后减小,最后趋于稳定,前3a产出流量增幅超过总增幅量的3/4;忽略水流损失,将高估产热速率,采热初期甚至达到考虑水流损失时产热速率的3倍以上;考虑水流损失,产热速率呈先快速上升再趋于稳定后逐渐下降的趋势,最优开采时间为3a^11a;研究上、下岩层对产出温度的影响,仅考虑传热,采热寿命延长5.43%,同时考虑渗流传热时,采热寿命延长2.71%;采热前9a,水流损失占主导作用,即流入上、下岩层水流损失对产热速率的影响高于热补偿效应,开采10a后,热补偿效应占主导作用;同时考虑水流损失和热补偿效应得到的产热速率变化规律与实际工程更为符合,建议选择低渗透能力的上、下岩层延长增强型地热系统(EGS)运行时间。     

离子浓度及表面结构对岩石孔隙内水流动特性的影响

酸性环境引发的岩石孔隙表面溶解增加了孔隙内水溶液的盐离子浓度,破坏了孔隙的表面结构.本文采用分子动力学模拟的方法研究了纳米级岩石孔隙内水溶液的流动特性,分析了盐离子浓度和孔隙表面结构对水流速度分布的影响及原因.研究结果表明:纳米级岩石孔隙内的水溶液流动符合泊肃叶流动特性,流速呈"抛物线"分布;随盐离子浓度增加,水溶液内部氢键网络变得更为致密,水黏度随其呈线性增长;水溶液中离子浓度越大,孔隙表面对水流动的阻力越大,最大流速越小,速度分布的"抛物线"曲率半径越大;岩石孔隙表面结构的破坏改变了流动表面的粗糙程度,增加了孔隙表面对H2O分子的吸引力.随表面结构破坏程度的增大,水溶液在近壁区域的密度增大,流速降低;当表面破坏程度达到50%时,水溶液在近壁区域出现了明显的负边界滑移现象.     

苯并三氮唑对变压器绕组硫腐蚀抑制机理的分子模拟研究

硫腐蚀现象对变压器绝缘性能的破坏不容忽视,其中腐蚀性硫化物与铜绕组作用产生的硫化亚铜是导致绝缘性能降低的一个重要原因.在变压器内部无氧和有氧的情况下,硫化亚铜的生成主要是因为腐蚀性硫化物与铜或者与氧化亚铜发生了一系列反应.为深入研究铜缓蚀剂苯并三氮唑(BTA)对铜绕组的缓蚀机理,采用分子模拟技术分析了BTA分子与铜晶体的电荷密度分布和态密度分布;同时利用过渡态搜索、红外光谱研究了BTA分子与氧化亚铜的反应过程.结果表明:加入BTA分子前后的铜表面态密度分布发生了变化,并且利用电荷密度分布图说明BTA分子与铜晶体形成了配位键;同时分析了BTA与氧化亚铜反应前后红外光谱的差异,发现BTA分子的结构已发生改变,与氧化亚铜形成了共价键.以上模拟结果说明BTA分子对铜和氧化亚铜有保护作用,能够抑制腐蚀性硫化物对铜绕组的腐蚀.     

硫原子在镍基合金825(001)面吸附的第一性原理研究

元素硫在镍基合金表面吸附产生严重的电化学腐蚀,为从原子尺度研究硫腐蚀机理,采用第一性原理方法,构建并优化了镍基合金825的晶胞结构模型,计算分析了S原子在镍基合金825耐蚀性较差面(001)晶面的吸附及电子转移情况.结果表明:Ni原子占据顶角, Cr原子和Fe原子对称占据面心是镍基合金825稳定的晶胞结构;原子S在镍基合金825(001)面上最稳定的吸附位为四重穴位,吸附能为-6.51 eV; S吸附前后的态密度(DOS)和二维电荷差分密度图(DCD)对比发现,镍基合金825中Fe与S之间电荷偏移明显,形成离子键,易生成腐蚀产物Fe_xS_y. S的吸附对镍基合金825中Cr原子的电子分布影响不大,且合金中Cr和Ni抑制了合金中Fe与S之间的相互作用,从而提高了合金耐蚀性.     

Continuous catalytic pyrolysis of oily sludge using U-shape reactor for producing saturates-enriched light oil

Continuous catalytic pyrolysis of oily sludge was carried out in a special U-shape reactor for producing saturates-enriched light oil. The sludge underwent thermal pyrolysis first and then catalytic pyrolysis. During the thermal pyrolysis, chain hydrocarbons were first cracked and further polymerized into aromatics. The effect of temperatures (400–800 °C) on the products was investigated and the maximum liquid yield (67.7%) was obtained at 500 °C. High temperature promoted polymerization, thus the distribution of aromatics in the liquid product was increased and was more concentrated in polyaromatics at 800 °C. In the catalytic upgrading stage, dolomite was used as catalyst and aromatics were adsorbed on it, either aggregated or decomposed. As a result, a light oil product with 57.0% saturates was obtained at the residence time of 8.9?s due to the conversion of aromatics and heavy hydrocarbons into light aliphatic hydrocarbons such as straight chain hydrocarbons. Compared with the oil phase in the raw sludge sample, the content of saturates was increased by 45.0% and that of the asphaltenes was reduced by 88.5%. Meanwhile, the inherent moisture in the oily sludge could participate in the steam reforming reaction, promoting the decomposition of aromatics and leading to an increase in the H₂ generation. Moreover, the release of H₂S was reduced from 0.132 to 0.005?mol per kg sludge and the sulfur content of the oil product was also decreased in the presence of dolomite. The deactivation of dolomite can be attributed to the carbonization of CaO and deposition of polyaromatic coke on the catalyst surface.     

一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷，首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点，建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联，建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型，根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处，并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证，结果表明：拟合曲线与试验曲线吻合度较高，所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律，尤其对加速蠕变阶段的非线性特征描述效果很好，验证了模型的合理性．     

岩石微观结构CT扫描表征技术研究

岩石微观结构特征研究对非常规油气藏的演化规律、油气赋存状态、运移方式、渗流特征等基础地质问题研究具有重要的科学意义.利用CT 扫描技术对油砂、致密砂岩和页岩样品微观结构表征研究,并对比了常规测试方法与CT扫描表征技术的差异性. 对比CT扫描法和筛分法测试油砂矿物的粒度分布,两者结果十分接近,矿物颗粒大小分布总体趋势上差异小,但小于96μm的颗粒矿物分布差异性稍大.CT扫描结果显示致密砂岩样品的裂缝比较发育,还发育溶蚀孔隙. 由于测试方法和样品大小的差异性,CT扫描获得的孔隙度略大于氦气法孔隙度.微米CT 扫描可以表征页岩的层理发育情况,但无法表征内部的微观孔隙结构.与常规测试方法相比,纳米CT 扫描表征页岩中有机质和黄铁矿的含量等方面准确性好,但孔隙度测试结果偏小.纳米CT分辨率还不能完全满足表征页岩的微观孔隙结构的要求,同时有机质和孔隙的灰度值差异较小,两者区分难度大,因此纳米CT还无法完全准确表征页岩微观孔隙结构.     

地下工程爆破对危岩稳定性的影响

根据爆破峰值速度衰减规律及危岩体主控结构面的假定，通过欧拉方程确定作用于危岩体峰值爆破荷载，并对峰值荷载进行修正得到时程爆破惯性力，建立了3类危岩体动力物理模型及计算模型，结合极限平衡方法得到3类危岩体动力稳定系数计算方法，选取最小的稳定性系数作为整个危岩体的动力稳定性系数，最终给出了3类危岩动力稳定性评价方法。以重庆市万州区太白危岩体为例，选取5个代表性危岩体进行动力稳定性计算，结果表明，3类危岩体动力稳定系数计算方法能较好地评价和反映危岩体的动力性态。     

基于Oda裂隙结构张量的流纹岩各向异性试验研究

岩石内天然存在长度、倾角和形态不同的裂隙,造成岩石的各向异性特征。为揭示岩石内天然随机裂隙发育特征对岩石物理力学特性的影响规律,以泥巴山隧址区采集裂隙性流纹岩为研究对象,首先对试样裂隙进行素描统计分析;然后基于Oda裂隙结构张量,获得天然随机分布裂隙的几何统计参数;最后对裂隙性流纹岩试样分别进行单轴和常规三轴压缩试验,得到不同应力路径下流纹岩的应力-应变曲线及物理力学参数。分析Oda裂隙结构张量定义的各向异性参数与试验获得的力学参数之间的规律,研究结果表明:(1)Oda裂隙结构张量适用于天然随机分布裂隙的几何统计分析,各向异性参数A(F)越大,裂隙优势方向越明显;(2)单轴压缩下,随着各向异性参数I1和A(F)的增大,流纹岩各向异性程度增大,弹性模量减小,泊松比增大;(3)常规三轴压缩下,流纹岩弹性模量和泊松比随各向异性参数改变的规律较不明显,Oda裂隙结构张量不再适用。     

基于单元破裂的岩石裂纹扩展模拟方法

传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔-库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力-应变曲线各阶段之间的对应关系. 结果表明: 连续-非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.