Round-off error in long-term orbital integrations using multistep methods

Techniques for reducing roundoff error are compared by testing them on high-order Störmer and summetric multistep methods. The best technique for most applications is to write the equation in summed, function-evaluation form and to store the coefficients as rational numbers. A larger error reduction can be achieved by writing the equation in backward-difference form and performing some of the additions in extended precision, but this entails a larger cpu cost.     

Advances in Numerical Modeling of Astrophysical and Space Plasmas

Advances in the simulation of astrophysical and cosmic plasmas are the direct result of advances in computational capabilities, today consisting of new techniques such as multilevel concurrent simulation, multi-teraflop computational platforms and experimental facilities for producing and diagnosing plasmas under extreme conditions for the benchmarking of simulations. Examples of these are the treatment of mesoscalic plasma and the scaling to astrophysical and cosmic dimensions and the Accelerated Strategic Computing Initiative whose goal is to construct petaflop (10¹⁵ floating operations per second) computers, and pulsed power and laser inertial confinement plasmas where megajoules of energy are delivered to highly-diagnosed plasmas. This paper concentrates on the achievements to date in simulating and experimentally producing plasmas scaled to both astrophysical and cosmic plasma dimensions. A previous paper (Part I, Peratt, 1997) outlines the algorithms and computational growth.     

On the Possibility of Directional Analysis for Geo-neutrinos

The possibility of terrestrial antineutrino directionality studies is considered for future unloaded liquid scintillator detectors. Monte-Carlo simulations suggest that the measurable displacement between prompt and delayed antineutrino signals makes such studies possible. However, it is estimated that on the order of 1000 terrestrial antineutrino events are required to test the simplest models, demanding detectors of 100 kt size to collect sufficient data in a reasonable period of time.     

矿物在油气形成过程中的作用

模拟实验结果表明,矿物组成对煤岩样品的模拟实验结果具有很大的影响,不同的矿物以及同一种矿物在不同的有机质热演化阶段,对有机质形成石油和天然气的过程均具有非常明显的影响:高岭土矿物在高温热演化阶段对有机质形成气态烃的过程具有一定的催化能力;碳酸钙（CaCO3）在中—低热演化阶段能抑制气态物质的形成,而在较高的热演化阶段则对煤生成气态物质的过程具有一定的催化作用;四氧化三铁（Fe3O4）在主要的模拟实验阶段对煤生成气态物质的过程均具有一定的催化作用;硫化亚铁（FeS）在低温热演化阶段对天然气的形成具有一定的抑制作用,而在高温热演化阶段则对天然气的形成具有较明显的催化作用;元素硫（S）在不同的热演化阶段对天然气的形成过程均具有非常明显的催化作用。矿物介质不仅对气态产物的产率有比较明显的影响,而且对液态产物的产率也有较明显的影响。模拟实验结果表明:高岭土和四氧化三铁（Fe3O4）在中、低温模拟实验阶段对“氯仿沥青A”的形成具有较明显的催化作用;硫化亚铁（FeS）在整个模拟实验的全过程中对“氯仿沥青A”的形成均具有明显的催化作用;碳酸钙（CaCO3）矿物对“氯仿沥青A”的形成具有明显的抑制作用;元素硫（S）对“氯仿沥青A”的形成具有非常明显的催化作用。     

基于GIS技术的河流污染动态模拟系统

本文论述了稳态均匀流场中的一维河流水质扩散模型,并在此模型的基础上,基于GIS二次开发技术开发了河流污染动态模拟系统。该系统以乐山市涌斯江支流段的地理底图及基本数据为例,实现了对河流污染情况的动态模拟。通过对该系统的一系列操作,可以对河段任意位置污染物浓度进行自动查询和修改,从而为河流水质污染的监测和管理提供可视化参考依据。     

A Numerical Simulation Study of Temperature and Pressure Effects on the Breakthrough Pressure of CO2 in Unsaturated Low-permeability Rock Core

Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock, and it also plays important roles in safety and capacity of CO2 geological storage. Based on the published experimental results, we present numerical simulations on CO2 breakthrough pressure in unsaturated low-permeability rock under 9 multiple P-T conditions (which can keep CO2 in gaseous, liquid and supercritical states) and thus, a numerical method which can be used to accurately predict CO2 breakthrough pressure on rock-core scale is proposed. The simulation results show that CO2 breakthrough pressure and breakthrough time are exponential correlated with P-T conditions. Meanwhile, pressure has stronger effects on experimental results than that of temperature. Moreover, we performed sensitivity studies on the pore distribution index λ (0.6, 0.7, 0.8, and 0.9) in van Genuchten-Muale model. Results show that with the increase of λ, CO2 breakthrough pressure and breakthrough time both show decreasing trends. In other words, the larger the value of λ is, the better the permeability of the caprock is, and the worse the CO2 sealing capacity is. The numerical method established in this study can provide an important reference for the prediction of gas breakthrough pressure on rock-core scale and for related numerical studies.     

碧口地块东-北缘中、新生代构造隆升及演化：来自磷灰石和锆石裂变径迹的证据

青藏高原具有复杂的构造演化特征,该地区自中、新生代以来的构造隆升和构造演化机制一直是地质研究的热点。为精细刻画青藏高原板块、华北板块和华南板块之间的拼合关系及差异性隆升特征,对位于青藏高原东北端的碧口地块进行了磷灰石和锆石裂变径迹测试,以及热史模拟和岩石冷却速率计算。结果锆石和磷灰石裂变径迹年龄分别在（118±5～265±29）Ma和（29.0±2.7～54.0±7.0）Ma之间;碧口地块东北缘及北缘冷却速率接近,在3.125～3.448 ℃/Ma之间,东缘冷却速率相对较低,为2.041～2.273 ℃/Ma。结果表明,中、新生代以来,碧口地块及其周缘总体上经历了持续隆升过程,但不同地区隆升特征具有差异性：碧口地块北侧在早、中侏罗世（151±7）Ma经历了构造挤压和隆升过程;东部相对较晚,在晚侏罗世（143±11）Ma经历了构造隆升阶段;东北端在早白垩世才与华北板块拼接并进入持续构造隆升阶段。进入古近纪（54.0±7.0）Ma隆升阶段,即始新世早期后,碧口地块东缘在始新世中后期（44 Ma）开始发生构造隆升,北缘自渐新世中晚期（29～32 Ma）开始发生显著的构造隆升。上述区域在10 Ma（中新世晚期）共同进入快速隆升阶段。     

河间潜山地热资源开发方案数值模拟

古潜山地热资源具备岩溶孔隙发育程度高、热储面积厚度大、地热水储量大的优点。冀中坳陷内古潜山分布密集且地热资源丰富, 河间潜山位于冀中坳陷饶阳凹陷中东部, 具有良好的地热地质条件, 开发潜力巨大。本文基于河间潜山及其周缘地区测井资料、岩石热物性并进行了计算, 发现其地温梯度为29.8 ℃/km到44.5 ℃/km之间, 平均值为40.7 ℃/km。大地热流值介于64.8~80.6 mW/m2之间, 平均值为 73.4 mW/m2。通过水热耦合模拟方法模拟选定的地热资源有利区的温度变化, 结果发现河间潜山合理的开采井距为800 m, 合理开采量为60 L/s, 回灌温度为35 ℃, 总可开采量为6.32×1016 J, 单年可开采量为 6.32×1014 J, 可供暖面积为1.22×106 m2, 对于冀中坳陷潜山地热资源的开发利用具有一定的指导意义。     

页岩储层物性参数对CO2不同封存机制及封存量的影响

CO₂增强页岩气开采技术(CO₂-ESGR)一方面可提高CH4产量,另一方面又可实现CO₂地质封存。为了分析页岩储层物性参数对CO₂封存机制的影响,文章以鄂尔多斯盆地延长组页岩为研究对象,采用CMG-GEM软件建立双孔双渗均质模型,分析了CO₂-ESGR中页岩储层垂直渗透率与水平渗透率之比(Kv/Kh)、含水饱和度和孔隙度对不同CO₂封存机制封存量的影响;并设计了27组正交试验采用极差分析法比较了三种因素的影响程度。研究表明,Kv/Kh在0.1～1范围增大会增加不同CO₂封存机制的封存量,封存总量最大可增加69.96%,其中吸附封存量最大可增加97.96%,受到影响最大;含水饱和度在0～0.9范围增大引起CO₂封存总量先增加后减小,封存总量最大可减少67.12%,其中溶解封存量最大可减少83.35%,范围波动最大;页岩储层孔隙度在0.1～0.99范围增大会导致CO₂封存总量减少,封存总量最大...