Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system

Simulations of flow for a discrete-fracture model in fractured porous rocks have gradually become more practical, as a consequence of increased computer power and improved simulation and characterization techniques. Discrete-fracture models can be formulated in a lower-dimensional framework, where the fractures are modeled in a lower dimension than the matrix, or in an equi-dimensional form, where the fractures and the matrix have the same dimension.     

Variational iteration solving method for El Nino phenomenon atmospheric physics of nonlinear model

A class of El Nino atmospheric physics oscillation model is considered. The El Nino atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conceptual oscillator model should consider the variations of both the eastern and westem Pacific anomaly patterns. An El Nino atmospheric physics model is proposed using a method for the variational iteration theory. Using the variational iteration method, the approximate expansions of the solution of corresponding problem are constructed. That is, firstly, introducing a set of functional and accounting their variationals, the Lagrange multiplicators are counted, and then the variational iteration is defined, finally, the approximate solution is obtained. From approximate expansions of the solution, the zonal sea surface temperature anomaly in the equatorial eastern Pacific and the thermocline depth anomaly of the seaair oscillation for El Nino atmospheric physics model can be analyzed. El Nino is a very complicated natural phenomenon. Hence basic models need to be reduced for the sea-air oscillator and are solved. The variational iteration is a simple and valid approximate method.     

汶川地震前的b值变化

Gutenberg-Richter定理给出了地震频度随震级的分布特征.大量研究表明频度-震级关系的斜率（b值）在大地震孕育的过程中会出现减小.为了考察b值在汶川地震孕育过程中的时间演化特征,本文尝试基于破裂断层选取研究区域,考察了2000年1月至2008年4月间,汶川地震（M_S8.0）破裂区的地震活动性,并对该区域b值的变化进行了探讨.结果表明,在2005年中至2006年底,地震月频度及季度频度有一个较明显的下降.b值从2002年始至地震前呈现出一个长期趋势性减小;在地震前约半年,出现快速、显著的下降.b值的这一时间变化特征与其他研究者报道的日本东北M_W9.0级地震前的b值变化特征具有很高相似性,可能反映了大地震准备过程中的应力变化.以上结果有益于认识和理解大地震孕育演化过程,同时也表明b值在中长期地震灾害评估中具有潜在价值.

     

自由站立式立管总体运动疲劳分析

由于良好的疲劳特性,自由站立式立管(FSHR)正广泛应用于深海油气田的开发当中。首先根据设计参数建立FSHR详细的有限元模型,然后利用谱分析方法对FSHR进行总体运动疲劳分析,以确定FSHR总体设计参数是否满足运动疲劳要求。最后对工作于相同环境下的钢悬链立管(SCR)进行总体运动疲劳分析,探讨哪种立管系统的疲劳性能更为优良。计算结果表明FSHR系统中刚性主管最大运动疲劳损伤出现在刚性主管顶部,FSHR总体设计参数满足规范要求;FSHR总体运动疲劳寿命要远大于SCR,体现了FSHR系统对浮体运动与刚性主管之间具有良好的解耦作用。     

基于多体分析的浅水FPSO和水下软钢臂系泊系统运动特性研究

基于多体动力分析方法进行FPSO和水下软钢臂系泊系统的运动特性研究。相较于非线性弹簧模拟软钢臂系泊系统或者其他近似模拟方法,多体分析方法可以充分考虑系泊系统具体结构形式及其动力项对FPSO运动性能的影响,更好的预报系统运动响应和系泊力。本文将FPSO和水下钢臂结构模拟成2个具有6自由度的独立结构,两者用系泊链组进行连接。基于三维势流理论应用汇源分布法,首先在频域内进行FPSO的水动力参数分析,进而在时域内对系统进行耦合动力分析。本文重点讨论系泊系统黏性力和二阶波浪力对系统响应的影响,计算结果发现系泊系统黏性力对系泊力有一定影响,而在浅水条件下二阶波浪力的计算对准确预报系统运动及系泊力非常重要。     

传统填图法在半覆盖区的改进——第四系岩性填图法

传统地质填图方法一般按成因划分第四系填图单位,但在半覆盖区,往往存在第四系不同成因物质混杂堆积的现象,若按成因划分填图单位,则会造成混乱。笔者在填图实践中,尝试按岩性划分第四系填图单位,这一做法能避免这种混乱,且有利于第四系中基岩信息的发掘,并能为勘查地球化学工作提供景观依据,这种新式填图方法可称为"第四系岩性填图法"。用岩性划分第四系填图单位需要使用统一的分类命名方案,本文使用了按碎屑特征进行分类的方案,主要考虑第四系松散堆积中碎屑的粒径、形态和含量变化。第四系岩性填图法中填图单位的划分需要参照这一分类命名方案,并可选用颜色、构造、次要碎屑类型等修饰性要素做为辅助划分依据。填图单位划分一般在踏勘和剖面实测中完成,第四系岩性填图法剖面测制较为复杂,除采用传统垂向测制方法外,还需运用横向测法和横向-垂向联合测法,通过剖面工作,可细化填图单位划分方案,并可评价不同填图单元基岩恢复的可靠性。第四系岩性填图法除需要编制传统地质图外,还需要编制基岩推断地质图或是在报告中制作插图以阐述此成果。     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

τ-P域地震资料三维叠前深度偏移及其实用算法

本文提出了在τ-P域实现三维叠前深度偏移的计算方法，该方法把地震道集由时-空域变换到τ-P域，然后对共P数据体进行相移或相移加插值偏移，将偏移后的各共P数据体叠加得到叠前深度偏移结果．该算法能够实现陡倾地层的正确归位，降低偏移处理维数，减少计算机内存需求，易于并行处理，但是计算量和所花费的计算机时间仍然很大．为此我们又提出了针对目的层的分时（层）偏移方法，由于只考虑目的层，可节省大量的计算时间．文中以二维实际地震资料及三维数值模型资料的处理说明了上述方法的可行性与有效性．     

基于Neuman井流理论确定玄武岩潜水含水层渗透系数的应用——以吉林靖宇煤矿为例

全面分析研究区玄武岩潜水含水层水文地质条件的基础上,对稳定流和非稳定理论确定潜水含水层渗透系数的优缺点进行了对比,认为非稳定井流公式确定玄武岩潜水含水层渗透系数是合理的。Aquifer Test软件基于Neuman井流理论配线法求参原理,运用计算机完成了自动求参,避免了人工配线的不足。该文根据吉林靖宇煤矿研究区玄武岩潜水含水层非稳定流抽水试验,分别用Aquifer Test软件和稳定井流公式进行了求参,并对比分析得出:玄武岩含水层弹性释水作用影响时间较长,利用Neuman潜水井流公式求参可靠,且更逼近现实;用Aquifer Test软件实现了自动求参,具有操作简单规范、可比性好的优点。     

新型贴井壁式阵列方位侧向测井数值模拟与响应特性

在现有侧向测井仪器基础上,提出了一种新型贴井壁式阵列方位侧向测井电极系,该仪器能提供径向和周向方位电阻率测量.该电极系采用贴井壁测量方式,提供5种径向探测深度.利用三维有限元方法模拟了阵列方位侧向测井电极系的井眼影响特性、径向探测特性、纵向分层能力、方位分辨能力,并模拟其对井周地层、水平井非对称泥浆侵入和倾斜地层的响应.在导电泥浆中最大探测深度为1.23m,纵向分辨率为0.3m,可以识别出0.1m薄层,方位分辨率为20°.贴井壁测量时,纵向分辨率不受泥浆和围岩电阻率的影响,能够准确测量井周方位电阻率,较不贴井壁测量具有很大优势,同时利用12条方位电阻率曲线能够反映出水平井泥浆非对称侵入特性,倾斜地层倾角和倾斜方向.