Beyond dual-porosity modeling for the simulation of complex flow mechanisms in shale reservoirs

The state of the art of modeling fluid flow in shale reservoirs is dominated by dual-porosity models which divide the reservoirs into matrix blocks that significantly contribute to fluid storage and fracture networks which principally control flow capacity. However, recent extensive microscopic studies reveal that there exist massive micro- and nano-pore systems in shale matrices. Because of this, the actual flow mechanisms in shale reservoirs are considerably more complex than can be simulated by the conventional dual-porosity models and Darcy’s law. Therefore, a model capturing multiple pore scales and flow can provide a better understanding of the complex flow mechanisms occurring in these reservoirs. This paper presents a micro-scale multiple-porosity model for fluid flow in shale reservoirs by capturing the dynamics occurring in three porosity systems: inorganic matter, organic matter (mainly kerogen), and natural fractures. Inorganic and organic portions of shale matrix are treated as sub-blocks with different attributes, such as wettability and pore structures. In kerogen, gas desorption and diffusion are the dominant physics. Since the flow regimes are sensitive to pore size, the effects of nano-pores and micro-pores in kerogen are incorporated into the simulator. The multiple-porosity model is built upon a unique tool for simulating general multiple-porosity systems in which several porosity systems may be tied to each other through arbitrary connectivities. This new model allows us to better understand complex flow mechanisms and eventually is extended into the reservoir scale through upscaling techniques. Sensitivity studies on the contributions of the different flow mechanisms and kerogen properties give some insight as to their importance. Results also include a comparison of the conventional dual-porosity treatment and show that significant differences in fluid distributions and dynamics are obtained with the improved multiple-porosity simulation.     

Adaptive particle swarm optimization for optimal orbital elements of binary stars

The paper presents an adaptive particle swarm optimization (APSO) as an alternative method to determine the optimal orbital elements of the star η Bootis of MK type G0 IV. The proposed algorithm transforms the problem of finding periodic orbits into the problem of detecting global minimizers as a function, to get a best fit of Keplerian and Phase curves. The experimental results demonstrate that the proposed approach of APSO generally more accurate than the standard particle swarm optimization (PSO) and other published optimization algorithms, in terms of solution accuracy, convergence speed and algorithm reliability.     

Micrometeorological Observations of a Microburst in Southern Finland

On the afternoon of 3 July 2004 in Hyytiälä (Juupajoki, Finland), convective cells produced a strong downburst causing forest damage. The SMEAR II field station, situated near the damage site, enabled a unique micrometeorological analysis of a microburst with differences above and inside the canopy. At the time of the event, a squall line associated with a cold front was crossing Hyytiälä with a reflectivity maximum in the middle of the squall line. A bow echo, rear-inflow notch, and probable mesovortex were observed in radar data. The bow echo moved west-north-west, and its apex travelled just north of Hyytiälä. The turbulence data were analysed at two locations above the forest canopy and at one location at sub-canopy. At 1412 EET (Eastern European Time, UTC+2), the horizontal and vertical wind speed increased and the wind veered, reflecting the arrival of a gust front. At the same time, the carbon dioxide concentration increased due to turbulent mixing, the temperature decreased due to cold air flow from aloft and aerosol particle concentration decreased due to rain scavenging. An increase in the number concentration of ultra-fine particles (< 10 nm) was detected, supporting the new particle formation either from cloud outflow or due to rain. Five minutes after the gust front (1417 EET), strong horizontal and downward vertical wind speed gusts occurred with maxima of 22 and 15 m s^?1, respectively, reflecting the microburst. The turbulence spectra before, during and after the event were consistent with traditional turbulence spectral theory.     

Phytoplankton productivity in the western subarctic gyre of the North Pacific in early summer 2006

We deployed a profiling buoy system incorporating a fast repetition rate fluorometer in the western subarctic Pacific and carried out time-series observations of phytoplankton productivity from 9 June to 15 July 2006. The chlorophyll a (Chl a) biomass integrated over the euphotic layer was as high as 45–50 mg Chl a m⁻² in the middle of June and remained in the 30–40 mg Chl a m⁻² range during the rest of observation period; day-to-day variation in Chl a biomass was relatively small. The daily net primary productivity integrated over the euphotic layer ranged from 144 to 919 mg C m⁻² day⁻¹ and varied greatly, depending more on insolation rather than Chl a biomass. In addition, we found that part of primary production was exported to a 150-m depth within 2 days, indicating that the variations in primary productivity quickly influenced the organic carbon flux from the upper ocean. Our results suggest that the short-term variability in primary productivity is one of the key factors controlling the carbon cycle in the surface ocean in the western subarctic Pacific.     

Hydrothermal plume mapping as a prospecting tool for seafloor sulfide deposits: a case study at the Zouyu-1 and Zouyu-2 hydrothermal fields in the southern Mid-Atlantic Ridge

Seafloor hydrothermal polymetallic sulfide deposits are a new type of resource, with great potential economic value and good prospect development. This paper discusses turbidity, oxidation–reduction potential, and temperature anomalies of hydrothermal plumes from the Zouyu-1 and Zouyu-2 hydrothermal fields on the southern Mid-Atlantic Ridge. We use the known location of these vent fields and plume data collected in multiple years (2009, 2011, 2013) to demonstrate how real-time plume exploration can be used to locate active vent fields, and thus associated sulfide deposits. Turbidity anomalies can be detected 10 s of km from an active source, but the location precision is no better than a few kilometers because fine-grained particles are quasi-conservative over periods of many days. Temperature and oxidation–reduction potential anomalies provide location precision of a few hundred meters. Temperature anomalies are generally weak and difficult to reliably detect, except by chance encounters of a buoyant plume. Oxidation–reduction potential is highly sensitive (nmol concentrations of reduced hydrothermal chemicals) to discharges of all temperatures and responds immediately to a plume encounter. Real-time surveys using continuous tows of turbidity and oxidation–reduction potential sensors offer the most efficient and precise surface ship exploration presently possible.     

Timing of the initial collision between the Indian and Asian continents

There exist three mainstream opinions regarding the timing of the initial collision between the Indian and Eurasian continents,namely,65±5,45±5,and 30±5 Ma.Five criteria are proposed for determining which tectonic event was related to the initial collision between India and Asia:the rapid decrease in the rate of plate motion,the cessation of magmatic activity originating from the subduction of oceanic crust,the end of sedimentation of oceanic facies,the occurrence of intracontinental deformation,and the exchange of sediments sourced from two continents.These criteria are used to constrain the nature of these tectonic events.It is proposed that the 65±5 Ma tectonic event is consistent with some of the criteria,but the upshot of this model is that the magmatic activity originating from the Tethyan subduction since the Mesozoic restarted along the southern margin of the Asian continent in this time after a brief calm,implying that the subduction of the Neotethys slab was still taking place.The magmatic activity that occurred along the southern margin of the Asian continent had a 7-Myr break during 72-65 Ma,which in this study is interpreted as having resulted from tectonic transformation from subduction to transform faulting,indicating that the convergence between the Indian and Asian continents was once dominated by strike-slip motion.The 30±5 Ma tectonic event resulted in the uplift of the Tibetan Plateau,which was related to the late stage of the convergence between these two continents,namely,a hard collision.The 45±5 Ma tectonic event is in accordance with most of the criteria,corresponding to the initial collision between these two continents.     

Effects of cold eddy on phytoplankton production and assemblages in Luzon strait bordering the South China Sea

The biochemical effects of a cold-core eddy that was shed from the Kuroshio Current at the Luzon Strait bordering the South China Sea (SCS) were studied in late spring, a relatively unproductive season in the SCS. The extent of the eddy was determined by time-series images of SeaWiFS ocean color, AVHRR sea surface temperature, and TOPEX/Jason-1 sea surface height anomaly. Nutrient budgets, nitrate-based new production, primary production, and phytoplankton assemblages were compared between the eddy and its surrounding Kuroshio and SCS waters. The enhanced productivity in the eddy was comparable to wintertime productivity in the SCS basin, which is supported by upwelled subsurface nitrate under the prevailing Northeastern Monsoon. There were more Synechococcus, pico-eucaryotes, and diatoms, but less Trichodesmium in the surface water inside the eddy than outside. Prochlorococcus and Richelia intracellularis showed no spatial differences. Water column-integrated primary production (IPP) inside the eddy was 2–3 times that outside the eddy in the SCS (1.09 vs. 0.59 g C m⁻²d⁻¹), as was nitrate-based new production (INP) (0.67 vs. 0.25 g C m⁻²d⁻¹). INP in the eddy was 6 times that in the Kuroshio (0.12 g C m⁻²d⁻¹). IPP and INP in the eddy were higher than the maximum production values ever measured in the SCS basin. Surface chlorophyll a concentration (0.40 mg m⁻³) in the eddy equaled the maximum concentration registered for the SCS basin and was higher than the wintertime average (0.29 ± 0.04 mg m⁻³). INP was 3.5 times as great and IPP was doubled in the eddy compared to the wintertime SCS basin. As cold core eddies form intermittently all year round as the Kuroshio invades the SCS, their effects on phytoplankton productivity and assemblages are likely to have important influences on the biogeochemical cycle of the region.     

Temporal variations in the current structure and volume transport of the Coastal Oyashio revealed by direct current measurement

Direct current measurements by a shipboard and bottom-mounted acoustic Doppler current profiler and concurrent hydrographic observations with a CTD were conducted off southeastern Hokkaido, Japan, between January and May 2005 to reveal temporal variations in the current structure and volume transport of the Coastal Oyashio (CO). The CO, which has a baroclinic jet structure with southwestward speeds exceeding 90 cm s^?1 and a width of 7–8 km, was associated with a surface-to-bottom density front and was formed on the offshore side of the shelf break. The volume transport of CO (T _CO) was estimated by integrating the fluxes of lower-density water that was trapped against the coast along the density front represented by the 26.2 σ _θ isopycnal line. This transport decreased monotonously from 0.79 Sv (1 Sv = 10⁶ m³ s^?1) in January to 0.21 Sv in March and subsequently to 0.12 Sv in May, possibly due to the decay of the East Sakhalin Current Water in the Okhotsk Sea. Accompanied by a decrease in T _CO, the location of the jet structure associated with the density front moved toward the coast while the maximum speed of the jet decreased and the tilt of the front became more horizontal. Consequently, more saline offshore Oyashio water flowed into the deep part of the shelf area, and the current structure altered from relatively barotropic in winter to baroclinic in spring. This study is the first to estimate the observed volume transport of the CO from direct current measurements.     

Applications of a GIS-based geotechnical tool to assess spatial earthquake hazards in an urban area

A geotechnical information system (GTIS) was constructed within a spatial geographic information system (GIS) framework to reliably predict geotechnical information and accurately estimate site effects at Gyeongju, an urban area in South Korea. The system was built based on both collected and performed site investigation data in addition to acquired geo-knowledge data. Seismic zoning maps were constructed using the site period (T _G) and mean shear wave velocity to a depth of 30 m (V _S30), and these maps were presented as a regional strategy to mitigate earthquake-induced risks in the study area. In particular, the T _G distribution map indicated the susceptibility to ground motion resonance in periods ranging from 0.2 to 0.5 s and the corresponding seismic vulnerability of buildings with two to five stories. Seismic zonation of site classification according to V _S30 values was also performed to determine the site amplification coefficients for seismic design and seismic performance evaluation at any site and administrative subunit in the study area. In addition, we investigated the site effects according to subsurface and surface ground irregularities at Gyeongju by seismic response analyses in time domains based on both two- and three-dimensional spatial finite element models, which were generated using spatial interface coordinates between geotechnical subsurface layers predicted by the GTIS. This practical study verified that spatial GIS-based geotechnical information can be a very useful resource in determining how to best mitigate seismic hazards, particularly in urban areas.     

ESTIMATE OF METHANE EMISSIONS FROM RICE FIELDS IN CHINA BY CLIMATE-BASED NET PRIMARY PRODUCTIVITY

1INTRODUCTIONIrrigated ricefieldsarecharacterizedbylargespatialandtemporalvariationsin CH4 emissiontotheatmo-sphere.Accordingly,thereisagreatuncertaintyintheestimate ofCH4 emissionsfromricefields.GreateffortshavebeenmadetoestimatetheCH4 emissionsfromricefieldsandseveralapproacheshavebeendeveloped.TherepresentativemethodsincludetheIPCC(Inter-governmentPanelofClimateChange)region-specificemissionfactormethodandthemodelcalculationmethod.Toimprovethecalculationaccuracy,theIPCCmethodreq…