Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.     

中国古代牧业社会骨制品的初步考察:以新疆巴里坤石人子沟遗址为例

考古遗址出土骨制品的研究对于揭示古代社会的动物资源获取和利用、手工业生产、社会组织结构等问题具有重要意义。国内骨制品研究目前多集中在农业区域,这些研究为探讨新石器时代至青铜时代的动物使用及其与社会发展、早期国家形成的关系等问题作出了重要贡献。然而,针对牧业社会骨制品的考察十分缺乏。新疆哈密地区巴里坤草原分布有大量古代牧业文化遗址,石人子沟遗址（43°31'12.8"~43°34'28.9"N,93°13'44.8"~93°16'49.1"E）是其中一处青铜时代晚期至铁器时代早期的大型聚落。本文从动物考古学视角研究该遗址2006年至2011年发掘出土的426件骨制品。结果显示,石人子沟遗址骨制品的原料主要为以羊（Ovis aries/Capra hircus）、马（Equus caballus）为主的家养动物和以鹿（Cervidae）为主的野生动物,其中羊的比例最高（69.7%）,鹿（13.8%）、马（2.3%）次之,这与中原地区青铜时代农业文化遗址的骨制品多以牛为原料的情况明显不同。羊在骨制品原料中占绝对多数且大量使用羊距骨制品的现象是对石人子沟遗址以牧业为主的生业经济方式的直接反映。石人子沟遗址不同类型骨制品的制作各具特点,但整体表现出"省时省力"的特点,即对使用部位细致打磨,对非使用部位仅做简单处理。与中原地区商周时期大型制骨作坊规范化、规模化、产业化的骨器生产不同,石人子沟遗址未见专门的制骨场所,骨制品生产的操作链条也并不完整。遗址的骨制品生产可能是以家庭为单元进行的,产品的专业化、精细化程度也相对较低。该研究填补了我国古代牧业文化遗址出土骨制品动物考古学专门研究的空白,为进一步探究我国古代不同区域、不同生业经济基础下的制骨手工业面貌提供了重要资料。

     

Comparison of two drought indices in studying regional meteorological drought events in China

The composite-drought index (CI), improved weighted average of precipitation index (IWAP), and the objective identification technique for regional extreme events (OITREE) were employed to detect China's regional meteorological drought events (CRMDEs) during 1961–2010. Compared with existing references, CI and IWAP both showed strong ability in identifying CRMDEs. Generally, the results of CI and IWAP were consistent, especially for extreme and severe CRMDEs. During 1961–2010, although the frequencies of extreme and severe CRMDEs based on CI and IWAP both showed weak decreasing trends, the two mean-integrated indices both showed increasing but not significant trends. However, the results of IWAP were more reasonable than CI’s in two aspects. Firstly, the monthly frequency of extreme and severe CRMDEs based on IWAP showed a clear seasonal variation, which coincided with the seasonal variation of the East Asian monsoon over central–eastern China, whereas the frequency based on CI presented a much weaker seasonal variation. Secondly, the two sets of results were sometimes inconsistent with respect to the start and end times of a CRMDE, and CRMDEs based on CI generally showed two unreasonable phenomena: (1) under non-drought conditions, a severe drought stage could suddenly occur in a large area; and (2) during the following period, drought could alleviate gradually in cases of non-precipitation. Comparative analysis suggested that the IWAP drought index possesses obvious advantages in detecting and monitoring regional drought events.     

近50年中国地面气候变化基本特征

采用国家基准气候站和基本气象站的地面资料,系统地分析了中国大陆地区1951年以来近地表主要气候要素演化的时间和空间特征。结果表明,中国近50 a来年平均地表气温变暖幅度约为1.1℃,增温速率接近0.22℃/(10 a),比全球或半球同期平均增温速率明显偏高。地表气温增暖主要发生在最近的20余年,其季节和空间特征与前人分析结论基本一致。降水量变化趋势对所取时间段和区域范围敏感。1951年以来全国平均降水量变化趋势不明显,但1956年以来略有增加。降水变化的空间特征明显而相对稳定,东北北部、包括长江中下游的东南部地区和西部广大地区降水增加,而华北地区以及东北东南部和西北东部地区降水明显减少。分析还发现,近50a来全国平均的日照时数、平均风速、水面蒸发等气候要素均呈显著下降趋势,但积雪地带的最大积雪深度却有所增加。中国日照时间和水面蒸发量变化的空间特征很相似,减少最明显的地区均发生在华北和华东,新疆次之。影响中国年代以上尺度气候变化的因子错综复杂,人类活动引起的大气中温室气体浓度增高可能在一定程度上影响了中国近50 a来的气候,但考虑到尚存的不确定性,目前仍不能给出明确结论。中国东部大部分地区日照时间和水面蒸发量减少可能均起源于人为排放的气溶胶影响,平均风速减弱也有利于水面蒸发量下降,而在西部地区云量和降水量的变化可能更重要。     

Geological Characteristics and Stability Evaluation of Wanjia Middle School Slope in Wenchuan Earthquake Area

To investigate the formation mechanism and the stability of Wanjia middle school slope in Wenchuan Earthquake Area, the macroscopic geological characteristics and the failure process of the landslide are researched by engineering geology analysis method, limit equilibrium method, and finit element method. The results show that after the Wenchuan Earthquake, retaining walls, houses and other infrastructure on the foot of Wanjia middle school slope were severely destroyed, 10 cm wide tension fracture appeared at the trailing edge of the slope. Wanjia middle school slope is a type of medium-sized soil landslide. The area of the deformation body is about 19,314 m², the total volume of the deformation body is about 23 × 10⁴ m³. There may be two potential sliding surfaces in the unstable slope: shallow and deep landslide. The analysis results of the limit equilibrium method and the finite element method show that: under dead weight, dead weight + rainstorm, dead weight + earthquake conditions, the plastic zone occurs mainly at the middle part or the trailing edge of the slope, and it doesn’t fully cut through the deep landslide body, so the deep landslide is stable. However, under rainstorm or earthquake conditions, the plastic zone almost completely cut thorough the shallow landslide body, it shows that the shallow landslide is in the understable–basic stable state. It is found that the results of finite element method is concordant with the results of the limit equilibrium method (F _s = 1.06–1.29, the shallow landslide is in the basic stable–stable state). The calculation results show that shallow landslides are likely to occur in Wanjia middle school slope during a rainstorm or an earthquake, so monitoring and control of the slope should be strengthened. The shallow landslide should be managed by some measures, such as anti slide pile retaining structures and drainage works, and the dangerous rock bodies on the slope surface should be cleaned up.     

Recurrent mesoproterozoic continental magmatism in South-Central Norway

We report U–Pb dates and Lu–Hf isotope data, obtained by LAM-ICPMS, for zircons from metamorphic rocks of the Setesdalen valley, situated in the Telemark block south of the classic Telemark region of southern Norway. The samples include infracrustal rocks from the metamorphic basement, metaigneous rocks and metasediments from the Byglandsfjorden supracrustal cover sequence, and metaigneous rocks which intruded the whole succession. The main crustal evolution took place from 1,550–1,020 Ma, beginning with the emplacement of juvenile tonalitic melts; the contribution of older crustal material increased with time. Around 1,320 Ma, further addition of juvenile material occurred, involving both mafic and felsic melts, metamorphism and deformation. Acid magmas with high FeO*/MgO were intruded at 1,215 Ma, coinciding with underplating elsewhere in South Norway. The period starting at 1,215 Ma is represented by supracrustal rocks, principally metarhyolites with minor mafic material and immature sediments of the Byglandsfjorden Group. The crust generation processes ended with the intrusion of diorites and granodiorites at 1,030 Ma, late in the Sveconorwegian orogeny. Regional processes of metamorphism and deformation (around 1,290 and 1,000 Ma) can be related to the assembly of Rodinia. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

山东泰安朱家庄自然硫矿床地质概要

山东省泰安朱家庄自然硫矿床,赋存于鲁西断隆古近纪汶东陆相盆地。自然硫矿层在汶口组二段上亚段,泥质灰岩-含泥灰岩-含云灰岩-石膏岩韵律层中的白云岩-石膏岩层段发育。自然硫以胶态、晶态、土态3种形态产出,可划分为顺层型、准顺层型、斑杂型及不顺层型4大类10余种矿石类型。从矿区去膏化作用普遍、自然硫与石膏-硬石膏分别富集轻-重同位素的特征、油气显示明显及地下水主要径流方向自然硫矿化最强等各种证据,表明该矿属生物后生成因的矿床。     

Broad search for unstable resonant orbits in the planar circular restricted three-body problem

Unstable resonant orbits in the circular restricted three-body problem have increasingly been used for trajectory design using optimization and invariant manifold techniques. In this study, several methods for computing these unstable resonant orbits are explored including grid searches, flyby maps, and continuation. Families of orbits are computed focusing on orbits with multiple loops near the secondary in the Jupiter–Europa system, and their characteristics are explored. Different parameters such as period and stability are examined for each set of resonant orbits, and the continuation of several specific orbits is explored in more detail.     

Spatial decision support applications based on three-dimensional city models

The basic mathematic models, such as the statistic model, the time-serial model, the spatial dynamic model etc., and some typical analysis methods based on 3DCM are proposed and discussed. A few typical spatial decision making methods integrating the spatial analysis and the basic mathematical models are also introduced, e.g. visual impact assessment, dispersion of noise immissions, base station plan for wireless communication. In addition, a new idea of expectation of further applications and add-in-value service of 3DCM is promoted. As an example, the sunshine analysis is studied and some helpful conclusions are drawn.     

华北强地震余震震群应力场“窗口”效应

本文研究了华北最近地震活动期的三次强地震(邢台、海城、唐山)的余震序列,从实际观测资料与国内外的一些实验结果,提出余震震群的出现是在区域应力场增强的背景下余震区应力加强的效应。强震余震震群的起伏,可以作为监视华北地区中强以上地震孕育的“窗口”。文中以“多点场”观点讨论了余震震群与中强以土地震关系,认为强震与中强震在区域应力场中是处于支配与被支配的不同地位。