Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

煤岩对压裂裂缝长期导流能力影响的实验研究

煤层中一般都伴有煤层气,通过压裂可以将煤层气释放出来,一方面能够获得清洁能源,另一方面可以提高煤田的安全开采。煤岩层自身特点决定了其在压裂中应用的实验评价方法及结论应具有特殊性,本文主要针对煤层气井的煤岩层压裂裂缝长期导流能力进行实验评价研究。实验表明,与砂岩地层不同,煤岩的硬度较小,压裂中支撑剂嵌入情况较严重;煤岩易破碎,碎屑颗粒充填到支撑剂中;同时支撑剂存在破碎及吸附伤害等,多种作用导致支撑裂缝的导流能力降低。实验中,以煤岩层闭合压力为一个重要的参数,实验对比分析了在不同闭合压力下,铺砂浓度、支撑剂粒径、支撑剂种类及组合的选择、支撑剂颗粒的破碎与嵌入、压裂液残渣等因素对煤岩层导流能力的影响。本文的实验结论为,嵌入对煤岩层支撑裂缝的导流能力伤害很大,加大铺砂浓度、采用低破碎率支撑剂、降低压裂液伤害能在一定程度上提高煤层裂缝的导流能力,但影响程度随闭合压力的增加而逐渐降低。本文所得出的结论对今后煤层气的开发及设计施工具有一定的指导意义。     

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高     

利用Jason-1卫星确定全球平均海面高

介绍了卫星测高数据处理的基本方法,分析研究了卫星测高的主要误差,提出了Jason-1的数据编辑准则,并用其测高数据确定了全球平均海面高,与CLS01模型进行了比较分析.     

基于权变理论的滑坡灾害监测预报思路探讨

文章简要评述了地质灾害基础理论与应用技术发展现状、滑坡灾害多种监测预报判据的利弊。利用综合信息处理决策方法,提出了基于权变理论的滑坡灾害监测预报新思路。分析了滑坡成灾的权变特征、环境因素和决策因素,建立了滑坡灾害预报决策概念模型。进一步探讨了在预报决策中应遵循的动态性及满意性原则,为提高地质灾害监测预报理论的科学性提供了新的理论依据与技术途径。     

中国海相探矿权区块定量评价方法

在研究了我国海相油气探矿权区块基本油气地质条件、勘探程度、技术经济条件等要素和特点的基础上,给出了48个评价参数和归一化取值参照表,建立了探矿权区块评价工作流程、评分参数及计算公式。提出对我国海相勘探层探矿权区块评价应该注意:与国外海相碳酸盐岩评价的差别;与国内陆相碎屑岩评价的差别;成藏组合的划分;评价内容的有效性、不确定性和评价参数的可信度;成败经验的总结和勘探程度的研究。     

基于ArcGIS的动态符号化设计与实现

地图符号是地图的语言单位,通过对地图符号的解读,可以直观的了解地图所表达的地理信息。ArcGis的符号化功能对于测绘符号的绘制是有一定局限的。而其他大部分测图软件的符号系统基本上是针对各个软件的,不同的软件都有独立的符号系统,离开了软件平台,如不打散块状符号或不添加辅助线就不能为其它软件所调用,本文中提出的ArcSymbol插件实现了在ArcGis操作环境下无需打散或添加辅助线即可方便、快速、无损的查看和打印满足地形图图式要求的图形,减少空间数据库中的数据冗余。并且ArcSymbol支持Shape、Coverage、Geodatabase格式,不需要数据的转换就能轻松地实现图形符号化,实现符号化的无损性。