一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

基于广州实证的后改革时代中国城市住房权问题

社会公正是当前地理学的一个热点也是传统研究课题。住房问题往往是社会公正研究和实践的重点案例。中国住房制度改革之后,中国住房拥有率及人均住房面积均大幅度提高。然而在住房分配的过程中,不平等和不公正现象时有发生。如何区别平均和公正,如何看待中国住房改革中的社会公平与不公共存的现象,是了解中国新城市现象的一个重要途径。本文从近期西方政治学和地理学关于住房权（housing rights）、社会公正以及城市权利等文献的讨论出发,通过在广州的实证研究,运用结构方程模型和Pratt指数对涉及住房权主要方面的住房分异的发生途径及其公正性进行综合分析。实证研究说明在住房分配领域,存在着制度因素延续下来所导致的分异。一些制度因素已经转化为市场因素,直接或间接影响个人居住条件和住房所有权的获得,并通过后者进一步影响住房分异。与此同时,纯粹市场因素的影响则相对薄弱。上述结论有助于了解当前中国城市住房权的核心问题,认识社会公正的发展前景,以及探寻减少社会不公正的制度利器。     

卫星跟踪经纬仪的旋转建模及转向自决策算法研究

在对卫星跟踪经纬仪转动轨迹建模基础上,提出一种基于预报轨迹及历史轨迹的经纬仪转向自决策算法。此算法具有事先规划经纬仪旋转轨迹的特点,从而在保护硬件设备的同时避免了卫星跟踪经纬仪目前普遍采用的“断电制动/限位反转”策略所引起的观测弧段的损失。该算法可以应用于其他类似的依赖于卫星预报的卫星跟踪经纬仪。     

特超稠油油藏水平井蒸汽驱油物理模拟实验

以胜利油田A区块特超稠油油藏为目标区块,自主研发了高温高压二维比例物理模型,研究了不同开发阶段水平井蒸汽驱油机理,开展了油藏压力、井底干 度、注汽强度对水平井蒸汽驱的效果影响的研究。物模实验结果显示,蒸汽驱整个过程分为三个阶段：存水回采期、汽驱受效期、蒸汽突破期。存水回采期蒸汽腔向水平方向扩展,以水平驱替方式为主 ;汽驱受效期蒸汽腔向垂向和水平方向扩展,以水平驱替方式为主,纵向泄油为辅;后期蒸汽突破,蒸汽腔继续向上方扩展,以纵向泄油为主,实现油层整体动用。5 MPa时转蒸汽驱的采出程度高于7 MPa时转蒸汽驱的采出程度;蒸汽驱的井底蒸汽干度应不低于0.4,蒸汽干度越高,水平井蒸汽驱的开发效果越好;注汽强度为1.9 t/(d·ha·m)时温度场最为发育。同时,在高压下提高干度可实现水 平井蒸汽驱的有效开发。以上的实验结果较好地指导了胜利油田A区块特超稠油油藏水平井蒸汽驱开发,取得了明显的开发效果,该研究对于类似稠油和超稠油油藏水平井蒸汽驱开发具有重要的指导意义 。     

致密砂岩气藏可动流体赋存特征的微观地质因素: 以苏里格气田东部盒8段储层为例

以苏里格气田东部盒8段典型的致密砂岩气藏为例,运用核磁共振、恒速压汞、铸体薄片、物性、X-衍射等实验资料,探讨了影响可动流体赋存差异的微观地质 因素。结果表明,盒8段储层可动流体饱和度低,T₂谱分布均为左高峰右低峰的双峰态;黏土矿物的充填与孔隙类型是孔隙结构复杂的重要因素,孔隙结构是影响可动流体赋存特征的关键;面孔率、喉道 半径、孔喉半径比是影响可动流体饱和度的主要因素,有效孔隙体积、分选系数对可动流体饱和度影响明显,储层物性、黏土矿物、有效喉道体积、孔隙半径对可动流体饱和度影响较弱。     

南海北部神狐海域天然气水合物成藏模式研究

南海北部陆坡神狐海域是我国海洋天然气水合物勘探开发研究的重点靶区,独特的水合物成藏特征,难以利用当前观测到的沉积速率和流体流动条件对其成藏机理进行解释和量化说明,对其形成演化模式和控制因素尚不明确.本文构建了海洋天然气水合物形成演化过程的动力学模型,模型的主控参量为海底沉积速率和水流速率,以此计算了神狐海域天然气水合物聚集演化过程,并与饱和度的盐度测试值进行对比.最后,在研究神狐海域地质构造活动和水合物成藏动力学基础上建立了神狐天然气水合物形成演化模式.认为神狐海域当前的天然气水合物是在上新世末—更新世早期断裂体系水合物基础上继承演化而来的,神狐海域天然气水合物形成演化具有典型的二元模式.第一阶段水合物形成发生在距今1.5 Ma之前构造活动形成的断裂体系中,高达50 m/ka的孔隙水流动携带了大量的甲烷进入水合物稳定带,导致了水合物的快速生成,在4万年内形成了饱和度达20%的甲烷水合物;第二阶段发生在1.5 Ma以来,泥质粉砂沉积使沉积体渗透率骤减,0.7 m/ka的低速率水流使甲烷供给不足,在海底浅层新沉积体中无法生成水合物,仅在水合物稳定带底部有缓慢的水合物继承增长,并因此形成了神狐海域当前观测到的水合物产出特征,而且水合物资源量仍在减少.     

基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

Delineation of the northern limit of the Congo Craton based on spectral analysis and 2.5D modeling of aeromagnetic data in the Akonolinga-Mbama area,Cameroon

Aeromagnetic data of the Akonolinga-Mbama region are analyzed in order to elucidate the subsurface geology of the area. The available data in the form of a residual aeromagnetic map is interpreted as a vast magnetically quiet region, and complex zones which do not correlate with the surface geology of the region.Within the magnetically quiet zone, a high negative circular elongated anomaly zone seems to represent an intrusion of a plutonic rock into the metamorphic formations of the region. Spectral analysis and two-and-a-half dimensional (2½-D) modeling are used to estimate the depth of the causative bodies and determine the source rocks along three profiles crossing the suspected areas. Models from various zones of granitic intrusions are obtained, thereby proposing some shallow fault lines along zones of contact. This permits us to mark out the northern margin of the Congo Craton, thus enabling us to distinguish the cratonic formations from the Pan African fold belt. Part of the belt has been thrust over the northern portion of the Congo Craton in Cameroon.     

Bidirectional behavior of unreinforced masonry walls

Most of the studies related to the modeling of masonry structures have by far investigated either the in‐plane (IP) or the out‐of‐plane (OP) behavior of walls. However, seismic loads mostly impose simultaneous IP and OP demands on load‐bearing or shear masonry walls. Thus, there is a need to reconsider design equations of unreinforced masonry walls by taking into account bidirectional effects. The intent of this study is to investigate the bidirectional behavior of an unreinforced masonry wall with a typical aspect ratio under different displacement‐controlled loading directions making use of finite element analysis. For this purpose, the numerical procedure is first validated against the results of the tests on walls with different failure modes conducted by the authors. Afterward, the response of the wall systems is evaluated with increasing top displacement having different orientations. A set of 19 monotonic and three cyclic loading analyses are performed, and the results are discussed in terms of the variation of failure modes and load–displacement diagrams. Moreover, the results of wall capacity in each loading condition are compared with those of the ASCE41‐06 formulations. The results indicate that the direction of the resultant force, vectorial summation of IP and OP forces, of the wall is initially proportional to the ratio of stiffness in the IP and the OP directions. However, with the increase of damage, the resultant force direction inclines towards the wall's longitudinal direction regardless of the direction of the imposed displacement. Finally, recommendations are made for applicability of ASCE41‐06 formulations under different bidirectional loading conditions. Copyright © 2014 John Wiley & Sons, Ltd.