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.     

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.     

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.     

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.     

Energy distributions of the large-scale horizontal currents caused by wind in the baroclinic ocean

Wind is the main energy source for the generation of the internal waves and the ocean mixing. Wunsch[1] estimated that about 1 TW (1 TW = 1012 W) energy was transported into the ocean from the winds by us-ing the altimeter data. Watanabe et al.[2] numerically calculated that the mixing processes obtained 0.7 TW energy from the global wind, which afforded most of the energy needed by the maintenance of the Merid-ional Overturning Circulation (MOC). During the past 50 years, in the Norther…     

Comparison of debris-flow volume and activity under different formation conditions

Debris flows frequently occurred in Wenchuan earthquake region from 2008 to 2010, resulting in great damage to localities and being a prolonged threat to reconstruction. Forty three events' data including debris-flow volume, sediment volume and watershed area are analyzed and compared with other debris-flow events in Eastern Italian Alps, burned areas in USA and in Taiwan. The analysis reveals that there is a strong empirical relationship between debris-flow volume and loose materials volume in the earthquake region. In addition, the relationship between debris-flow volume and watershed area in the earthquake region has a wider variation range than that in other three regions while the debris volume also appears to be larger than that in the other three regions, which implies the volume of debris flows with strong influence of earthquakes is larger than that with no such influence and it is hard to predict the post-quake volume only by the watershed area. The comparison of the maximal debris-flow erosion modulus in the Wenchuan region and in Taiwan indicates that debris flows will be very active in a short time after strong earthquake.     

湘中邵阳凹陷二叠系龙潭组页岩含气性特征与气体成因

针对湘中地区邵阳凹陷二叠系龙潭组页岩气资源评价,部署了页岩气调查井2015H-D3井,通过现场解析气等相关样品测试定量分析了页岩含气性特征及其影响因素,并借助气体稳定碳、氢同位素对气体成因进行了初步探讨.结果显示:钻深从150m处开始出现气显,随深度增加,解析气含量呈现先逐渐增大后减小的趋势,300～425m为最高含气层段,累计厚度达125m,现场解析气含量全部大于0.5m^3/t,最高为2.35m^3/t,平均1m^3/t,证实龙潭组具良好页岩气资源潜力.当埋深大于300m时,解析气含量受有机碳含量控制,而埋深小于300m时,解析气含量并不简单受控于有机碳含量,而是受到保存条件的严重制约.气体碳、氢同位素测试显示δ13C1介于-29.87‰～-36.82‰,平均为-34.52‰,δ¹³C_2介于-29.45‰～-31.02‰,平均为-30.09‰,δD1介于-131.20‰～-178.40‰,平均为-167.40‰.气体成因分析揭示龙潭组页岩气属热成因中的油型气类型.基于沉积现象判断龙潭组为海陆过渡相环境,与氢同位素判断结果基本吻合,但目前还无法确定具体判定区间.     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

Saddle-shaped reticulate <Emphasis Type="Italic">Nummulites</Emphasis> from Early Oligocene rocks of Khari area,SW Kutch,India

Saddle-shaped reticulate Nummulites from the Early Oligocene rocks of Khari area, SW Kutch, India is reported here for the first time. Unusual shape of this Nummulites is due to the curved nature of the coiling plane, indicating space constrained postembryonic test growth. With regular development of chambers, septa and septal filaments, the saddle-shaped Nummulites constitutes the third morphotype of N. cf. fichteli Michelotti form A. Other morphotypes of the species reported earlier include inflated lenticular and conical tests. Multiple morphotypes of N. cf. fichteli form A indicates varied test growth in response to substrate conditions. Morphological variability exhibited by N. cf. fichteli form A from Kutch and some Early Oligocene reticulate Nummulites from the Far East are comparable. This faunal suite is morphologically distinct from the contemporary reticulate Nummulites of the European localities.