盐穴储油（气）库水溶造腔管柱动力特性初探

针对盐穴储气库单井水溶造腔过程中管柱动力特性开展了初步探讨,认为造腔工程中,管柱在受限空间的自激振动和动力失稳是导致管柱发生严重弯曲或者破损的主要原因。初步分析了无空间约束条件下内外管柱的振动特性,得到了理想条件下管柱临界流速的计算式。计算表明,无约束条件下造腔管柱一般会由于液-固耦联振动而发生动力失稳,而且造腔内管一般会先于造腔外管发生失稳破坏。鉴于问题的重要性和复杂性,还对造腔管柱临界流速全面计算力学模型应考虑的诸因素进行了探讨,初步分析了管道轴力、液体压力、造腔循环方式、管道水平支撑,以及轴向水击和管道横向振动耦合的影响,为进一步结合工程实际开展空间约束条件下造腔管柱动力特性全面分析提出了建议。     

基于最小二乘法的边坡稳定性分析

利用经典土压力理论设定合理土条推力线范围,根据静力平衡及力矩平衡条件建立以条块界面安全系数 为变量的线性超定方程组,应用MATLAB软件基于最小二乘方法对此方程组求解,依据极大值原理及合理性条件： ,最终获得边坡整体安全系数 。经算例验证得出如下结论：当条块界面与条块整体安全系数相等时,边坡抗剪能力发挥最大;安全系数随推力线位置的升高而减小,但影响并不大;由于考虑了极大值原理,该法求得的安全系数略微偏大于其他方法。     

Integration rules from the equivalent differential equation method applied to equations of motion

The numerical integration of equations of motion necessarily implies the presence of errors that depend on initial conditions as well as the different physical parameters under consideration. More particularly, dumping or dissipative terms can appear and it is especially interesting to determine its causes. The equivalent differential equation method may allow the errors from a certain numerical scheme to be analyzed and, together with other considerations, can help us to eliminate or reduce them.     

不规则区域海浪嵌套模式设计及其应用实例

本研究针对第三代MASNUM海浪模式的并行设计特点,设计实现了并行化不规则嵌套计算,从而提高了海浪模式嵌套模拟的灵活性和时效性.选取浙江舟山、宁波附近海域,利用201 1年的NCEP再分析风场作为强迫,采用大、中、小网格并利用该嵌套模式进行了三重嵌套应用.其中大、中区域实验均采用规则矩形方式设置开边界,大区域的模拟结果为中区域嵌套实验提供开边界条件,同时中区域嵌套输出小区域嵌套实验在不规则开边界上的边界条件.中区域的实验结果与卫星观测资料的对比表明:嵌套实验比无嵌套实验的误差小30％左右,二者在开边界附近的有效波高差异较大,达到0.5m以上.8月份中区域的无嵌套实验所得有效波高平均为0.7m,最大值为3.6m,嵌套实验结果则分别为1.0、7.1m,后者更接近以往观测资料所显示的平均波高1.5～2.0 m和最大波高5.0 ～10.0 m的结果,表明该常规嵌套方式的可行性.小区域实验也包含嵌套和无嵌套2个实验,开边界均采用不规则形状.实验结果表明:在5～8月期间,嵌套实验的波向与无嵌套的波向差异较明显,主要表现为外海传入;在6～7月份,嵌套实验与无嵌套实验有效波高的差值达1.0m以上;嵌套实验的周期普遍大于无嵌套的模拟结果,其最大差异为8 s.上述实验结果表明,该嵌套模式可实现规则开边界区域和不规则区域的嵌套模拟,可应用于复杂近岸区域的海浪数值模拟和预报,有效提高对关注区域海浪的模拟能力.     

The Stability of Tide Gauges in the South Pacific Determined from Multiepoch Geodetic Levelling, 1992 to 2010

Tide gauge data is important for determining global or local sea level rise with respect to a global geocentric reference frame. Data from repeated precise levelling connections between the tide gauges and a series of coastal and inland benchmarks, including Continuous GPS (CGPS) benchmarks, are used to determine the stability of tide gauges at 12 locations in the South Pacific. The method for determining this stability is based on a constant velocity model which minimises the net movement amongst a set of datum benchmarks surveyed since the installation of the tide gauges. When assessed at a 95% confidence interval, and with the exception of the Solomon Islands, none of the tide gauges were found to be in motion relative to their CGPS benchmarks. The Solomon Islands estimate is considered to be unreliable since the CGPS benchmark was recently established and has been surveyed fewer than three times. In Tonga and Cook Islands, the tide gauges were found to be disturbed or affected by survey errors whereas the Vanuatu results were affected by earthquakes.     

基于最多叶子生成树的中国航空网络轴辐结构构建

航空网络的轴辐 (Hub-Spoke) 结构是实现规模经济发展的重要交通运输网络结构,本文为此提出了一种全新的航空网络轴辐结构构建方法。该方法从图论和地理学的角度出发,引入地理距离约束,改进了传统的最多叶子生成树 (Maximum Leaf Spanning Tree) 算法,直接从现有的中国航空网络中抽取树形轴辐结构形成航空支线网络,然后选取支线网络中度前10的节点作为航空枢纽点,并将枢纽点之间在原图中的航线抽取为航空干线网络,最后将支线网络和干线网络合并形成中国航空网络的轴辐结构。在与相关研究的对比分析中,本文方法虽是从图论角度出发,但构建的中国航空轴辐结构符合实际地理环境,划分支线网络距离阈值的选择更加客观合理,所选的航空枢纽点地理意义更为明显,干支线网络的覆盖度更为全面。     

Comparison of Ensemble Kalman Filter groundwater-data assimilation methods based on stochastic moment equations and Monte Carlo simulation

Traditional Ensemble Kalman Filter (EnKF) data assimilation requires computationally intensive Monte Carlo (MC) sampling, which suffers from filter inbreeding unless the number of simulations is large. Recently we proposed an alternative EnKF groundwater-data assimilation method that obviates the need for sampling and is free of inbreeding issues. In our new approach, theoretical ensemble moments are approximated directly by solving a system of corresponding stochastic groundwater flow equations. Like MC-based EnKF, our moment equations (ME) approach allows Bayesian updating of system states and parameters in real-time as new data become available. Here we compare the performances and accuracies of the two approaches on two-dimensional transient groundwater flow toward a well pumping water in a synthetic, randomly heterogeneous confined aquifer subject to prescribed head and flux boundary conditions.     

储层弹性与物性参数地震叠前同步反演的确定性优化方法

储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性.     

附加线性不等式约束的条件平差模型未知参数的解算

讨论附加线性不等式约束的条件平差模型的解算思路。根据K-T定理,对拉格朗日乘子λ进行约束,通过迭代算法求解满足K-T条件的拉格朗日乘子,进而求得未知参数的最佳估值,并用算例验证该方法的可行性。     

半拉格朗日模式风的半解析解及与差分解的比较试验

质点的轨迹计算是半拉格朗日模式的重要基础,传统的数值计算方法由于采用时间差分代替微分,只能得到质点运动轨迹终点的速度,因此质点的移动轨迹（位移）只能靠风速外推的方法计算,导致了模式计算不稳定等问题.借鉴精细积分法中使用半解析解的思路,利用正压原始方程研究了用运动方程的半解析解构建数值模式的可能性.求解了运动方程的一阶和二阶微分方程组的半解析解,通过时间积分半解析解计算质点运动轨迹.数值试验表明,一阶微分方程组的半解析解比差分解略有优势.二阶微分方程组的半解析解在时间步长增大时优势非常明显,而且在保证计算精度的前提下,节省计算时间,这对提高模式性能有重要作用.