Linear Elastic and Cohesive Fracture Analysis to Model Hydraulic Fracture in Brittle and Ductile Rocks

Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a Perkins-Kern–Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone effect on the ductile rock fracture.     

Spatial patterns and storage composition of woody debris in a natural secondary forest dominated by <Emphasis Type="Italic">Pinus tabulaeformis</Emphasis> on Loess Plateau,China

Woody debris (WD) is an important part of natural Pinus tabulaeformis mixed stands, and it affects the forest ecosystem stability and development. The WD spatial patterns are especially important structural characteristics that can provide insights into forest dynamics. In this paper, the WD storage, WD spatial patterns and WD associations among the main species were examined in the natural secondary forest on Loess Plateau in northwest China. Data were collected in a 1 ha (100 m × 100 m) permanent plot, and all the trees with a diameter at breast height of more than 3 cm were measured and stem-mapped. Ripley’s K functions from the spatial-point-pattern-analysis method were used to analyze the spatial distribution and associations. The results showed that: (1) The total storage of WD was 10.73 t/ha, fallen wood was the main source of WD, and the majority diameters were greater than 20 cm, and in intermediate levels of decay; (2) The overall spatial pattern was closely related to the spatial scale, which exhibited an aggregated pattern on a small scale, and a random pattern on a large scale. The spatial patterns of coarse woody debris also gradually transitioned from an aggregated pattern in fine scales to a random pattern in broader spatial scales, which matched the overall spatial pattern. The spatial intensity was gradually decreased with the increasing diameters, and increased with the decomposition classes; (3) The WD of Pinus tabulaeformis species was negatively associated with Betula platyphylla and Populus davidiana on a small scale but positively associated with these species on a large scale. The spatial pattern and interspecies relations were the results of long-term interactions between the natural secondary forest community and the surrounding natural environment. These findings would provide a scientific basis for the sustainable management and protection of natural secondary forest ecosystems on Loess Plateau.     

The origin and biogeochemistry of organic matter in surface sediments of Lake Shihwa and Lake Hwaong

To understand the origin and biogeochemistry of the organic matter in surface sediments of Lake Shihwa and Lake Hwaong, organic nitrogen, inorganic nitrogen, labile organic carbon, and residual organic carbon contents as well as stable isotope ratios for carbon and nitrogen were determined by KOBr-KOH treatment. Ratios of organic carbon to organic nitrogen (C_org/N_org) (mean = 24) were much higher than ratios of organic carbon to total nitrogen (C_org/N_tot) (mean=12), indicating the presence of significant amounts of inorganic nitrogen in the surface sediments of both lakes. Stable isotope ratios for organic nitrogen were, on average, 5.2‰ heavier than ratios of inorganic nitrogen in Lake Shihwa, but those same ratios were comparable in Lake Hwaong. This might be due to differences in the origin or the degree of degradation of sedimentary organic matter between the two lakes. In addition, stable isotope ratios for labile organic carbon were, on average, 1.4‰ heavier than those for residual organic carbon, reflecting the preferential oxidation of¹³C-enriched organic matter. The present study demonstrates that KOBr-KOH treatment of sedimentary organic matter can provide valuable information for understanding the origin and degradation state of organic matter in marine and brackish sediments. This also suggests that the ratio of C_org/N_org and stable isotope ratios for organic nitrogen can be used as indexes of the degree of degradation of organic matter.     

Vertical distributions of aerosol optical properties during haze and floating dust weather in Shanghai

A comparative study on the vertical distributions of aerosol optical properties during haze and floating dust weather in Shanghai was conducted based on the data obtained from a micro pulse lidar.There was a distinct difference in layer thickness and extinction coefficient under the two types of weather conditions.Aerosols were concentrated below 1 km and the aerosol extinction coefficients ranged from 0.25 to 1.50km^-1 on haze days.In contrast,aerosols with smaller extinction coefficients(0.20 0.35 km^-1) accumulated mainly from the surface to 2 km on floating dust days.The seasonal variations of extinction and aerosol optical depth（AOD） for both haze and floating dust cases were similar greatest in winter,smaller in spring,and smallest in autumn.More than 85%of the aerosols appeared in the atmosphere below 1 km during severe haze and floating dust weather.The diurnal variation of the extinction coefficient of haze exhibited a bimodal shape with two peaks in the morning or at noon,and at nightfall,respectively.The aerosol extinction coefficient gradually increased throughout the day during floating dust weather.Case studies showed that haze aerosols were generated from the surface and then lifted up,but floating dust aerosols were transported vertically from higher altitude to the surface.The AOD during floating dust weather was higher than that during haze.The boundary layer was more stable during haze than during floating dust weather.     

Two and Half Dimensional Simulation of Ridge Effects on the Ground Motion Characteristics

— Seismic responses of weathered and non-weathered ridge models were simulated to study the ridge effects on the ground motion characteristics. The range of ridge slope from 19.98° to 45° was considered to produce a possible set of generalized results. 2.5-D modeling based on parsimonious staggered grid approximation of elastodynamic wave equations was adopted in simulations. Computed results reveal an increase of amplitude of incoming waves with both elevation and the slope of the ridge. Further, the characteristics of surface waves are highly ridge slope dependent. The analysis of responses of weathered and non-weathered ridge models reveals that ridge has caused a strong generation of surface waves near its top. The surface waves are not dominating on the top of the ridge but at some lower elevation. The increase of weathering of ridge further intensified the ridge effect. Analysis in frequency domain, based on spectral ratio method, does not indicate any pattern in the spectral amplification factor and is very much sensitive to slope, source focal mechanism and location. However, on an average there is a continuous decrease of amplification with slope in the vertical component and increase in the transverse component, and it is increasing in the radial component up to slope =38.0° and thereafter decreasing.     

Calibrating the GOCE accelerations with star sensor data and a global gravity field model

A reliable and accurate gradiometer calibration is essential for the scientific return of the gravity field and steady-state ocean circulation explorer (GOCE) mission. This paper describes a new method for external calibration of the GOCE gradiometer accelerations. A global gravity field model in combination with star sensor quaternions is used to compute reference differential accelerations, which may be used to estimate various combinations of gradiometer scale factors, internal gradiometer misalignments and misalignments between star sensor and gradiometer. In many aspects, the new method is complementary to the GOCE in-flight calibration. In contrast to the in-flight calibration, which requires a satellite-shaking phase, the new method uses data from the nominal measurement phases. The results of a simulation study show that gradiometer scale factors can be estimated on a weekly basis with accuracies better than 2 × 10⁻³ for the ultrasensitive and 10⁻² for the less sensitive axes, which is compatible with the requirements of the gravity gradient error. Based on a 58-day data set, scale factors are found that can reduce the errors of the in-flight-calibrated measurements. The elements of the complete inverse calibration matrix, representing both the internal gradiometer misalignments and scale factors, can be estimated with accuracies in general better than 10⁻³.     

Decadal evolution of tidal flats and channels in the Outer Weser estuary,Germany

This study focuses on the medium scale morphodynamics of the tidal flat and channel system Fedderwarder Priel, located in the Outer Weser estuary (Wadden Sea, Germany). Tidal channels and adjacent flats are highly dynamic systems whose morphologic evolution are driven by tidal, wind, and wave forcings. These coastal environments are an important ecosystem and react to changes in hydrodynamic conditions in various spatial and temporal scales. Based on annual medium-resolution digital elevation models from 1998 to 2016, we describe changes in the surface area over depth with hypsometries and use vertical dynamic trends in order to analyze and visualize the morphologic evolution of the Fedderwarder Priel and adjacent tidal channels. It is shown that several intertidal flats rise in the order of 1.3 to 5.6 cm/year. The findings indicate that the Outer Weser estuary was not in an equilibrium state for the investigated period, and tidal flats accreted with a rate exceeding mean sea level rise.     

The diagenetic geochemistry and contamination assessment of iron,cadmium, and lead in the sediments from the Shuangtaizi estuary,China

Sediment and pore water samples have been collected from the coastal tidal flat in the Shuangtaizi estuary, China, in order to investigate the geochemical behavior of iron, cadmium, and lead during diagenesis and to assess the degree of contamination. The calculated enrichment factors and geoaccumulation indices for separate elements show that anthropogenic activities have had no significant influence on the distribution of Fe and Pb in the study area, whereas the distribution of Cd has been closely influenced in this way. The high percentage of exchangeable Cd (average of 56.34%) suggests that Cd represents a potential hazard to benthic organisms in the estuary. The calculated diffusive fluxes of metals show that the most mobilized metal is Fe (9.22 mg m^?2 a^?1), followed by Cd (0.54 mg m^?2 a^?1) and Pb (0.42 mg m^?2 a^?1). Low Fe²⁺ contents in surface pore water, alongside high chromium-reducible sulfur contents, and low acid-volatile sulfur, and elemental sulfur contents at 0–25 cm depth in sediments show that Fe²⁺ is formed by the reduction of Fe oxides and is transformed first to a solid phase of iron monosulfides (FeS) and eventually to pyrite (FeS₂). The release of adsorbed Pb due to reductive dissolution of Fe/Mn oxides during early diagenesis could be a source of Pb²⁺ in pore water. From the relatively low total organic carbon contents measured in sediments (0.46–1.28%, with an average of 0.94%) and the vertical variation of Cd²⁺ in pore water, sulfide or Fe/Mn oxides (instead of organic matter) are presumed to exert a significant influence on carrying or releasing Cd by the sediments.     

On the assessment of Argo float trajectory assimilation in the Mediterranean Forecasting System

The Mediterranean Forecasting System (MFS) has been operational for a decade, and is continuously providing forecasts and analyses for the region. These forecasts comprise local- and basin-scale information of the environmental state of the sea and can be useful for tracking oil spills and supporting search-and-rescue missions. Data assimilation is a widely used method to improve the forecast skill of operational models and, in this study, the three-dimensional variational (OceanVar) scheme has been extended to include Argo float trajectories, with the objective of constraining and ameliorating the numerical output primarily in terms of the intermediate velocity fields at 350 m depth. When adding new datasets, it is furthermore crucial to ensure that the extended OceanVar scheme does not decrease the performance of the assimilation of other observations, e.g., sea-level anomalies, temperature, and salinity. Numerical experiments were undertaken for a 3-year period (2005–2007), and it was concluded that the Argo float trajectory assimilation improves the quality of the forecasted trajectories with ~15%, thus, increasing the realism of the model. Furthermore, the MFS proved to maintain the forecast quality of the sea-surface height and mass fields after the extended assimilation scheme had been introduced. A comparison between the modeled velocity fields and independent surface drifter observations suggested that assimilating trajectories at intermediate depth could yield improved forecasts of the upper ocean currents.     

Transformation of Carbon and Nitrogen by Earthworms in the Decomposition Processes of Broad-leaved Litters

Earthworms are the important constituents in the decayed food web and the main ecological conditioners in the process of decomposition and nutrient mineralization. The transformation of organic carbon (C) and total nitrogen (N) in the broad-leaved litters ingested by earthworms was researched by means of a laboratory experiment. Experimental samples were collected from broad-leaved Korea Pine mixed forest in Liangshui National Natural Reserve (47°10′50″N, 128°53′20″E) in the northeastern Xiao Hinggan Mountains of Northeast China. The contents of organic C and total N in earthworms, leaf litters and earthworm faeces were analyzed. Results show that the organic C content was in the fol- lowing order: leaf litters>faeces>earthworms, while total N content was contrary to that of the organic C. The organic C contents in the different leaf litters were in the following order: Tilia amurensis>Betula costata>Acer mono, whereas the total N contents in the different leaf litters were: Betula costata>Tilia amurensis>Acer mono. The contents of organic C and total N in the faeces from the different leaf litters were almost consistent with the contents of the leaf litters. After the leaf litters were ingested by earthworms, the organic C, which was transformed to increase earthworms' weights, ac- counted for 3.90%-13.31% of the total ingestion by earthworms, while that in the earthworm faeces accounted for 6.14%-13.70%. The transformed organic C through the other metabolism (e.g., respiration) of earthworms accounted for 75.04%-89.92%. The ingested organic C by earthworms was mostly used for metabolic activities. The N ingested by earthworms was less than organic C. It is estimated that 37.08% of total N was transformed to increase the earthworm's weight, 19.97% into earthworm faeces and 47.86% for the consumption of the earthworm's activities. The earthworms not only increased the content of organic C and total N in the soil, but also decreased the values of C/N in the soil and leaf litters. Earthworms play a major role in the leaf litters' decomposition and transformation.