Shale oil recovery by CO_(2)injection in Jiyang Depression,Bohai Bay Basin,East China

Laboratory experiments,numerical simulations and fracturing technology were combined to address the problems in shale oil recovery by CO_(2)injection.The laboratory experiments were conducted to investigate the displacement mechanisms of shale oil extraction by CO_(2)injection,and the influences of CO_(2)pre-pad on shale mechanical properties.Numerical simulations were performed about influences of CO_(2)pre-pad fracturing and puff-n-huff for energy replenishment on the recovery efficiency.The findings obtained were applied to the field tests of CO_(2)pre-pad fracturing and single well puff-n-huff.The results show that the efficiency of CO_(2)puff-n-huff is affected by micro-and nano-scale effect,kerogen,adsorbed oil and so on,and a longer soaking time in a reasonable range leads to a higher exploitation degree of shale oil.In the"injection+soaking"stage,the exploitation degree of heavy hydrocarbons is enhanced by CO_(2)through its effects of solubility-diffusion and mass-transfer.In the"huff"stage,crude oil in large pores is displaced by CO_(2)to surrounding larger pores or bedding fractures and finally flows to the production well.The injection of CO_(2)pre-pad is conducive to keeping the rock brittle and reducing the fracture breakdown pressure,and the CO_(2)is liable to filter along the bedding surface,thereby creating a more complex fracture.Increasing the volume of CO_(2)pre-pad can improve the energizing effect,and enhance the replenishment of formation energy.Moreover,the oil recovery is more enhanced by CO_(2)huff-n-puff with the lower shale matrix permeability,the lower formation pressure,and the larger heavy hydrocarbon content.The field tests demonstrate a good performance with the pressure maintained well after CO_(2)pre-pad fracturing,the formation energy replenished effectively after CO_(2)huff-n-puff in a single well,and the well productivity improved.

Groundwater Potential Assessment and Characterization of Genale-Dawa River Basin

Genale-Dawa River Basin is one of the largest and most drought prone re-gions in Ethiopia. As a result a search for alternative source of water has always been a major issue in the region. This study therefore, aims at characterizing and evaluating the ground water potential resource of the river basin. The delineation and numerical discretization of the river basin were primarily done for the proper development of the problem domain. And the groundwater governing equation was solved at individual nodal points of the discretized region using a 3-D numerical ground water model called TAGSAC. This resulted in monthly average groundwater table fluctuation and replenishable ground water potential estimation. Hydrogeological classification was also done based on hydraulic conductivity values obtained from model calibration.

Quantity-based shipment consolidation and stock replenishment policy with lead time in vendor-managed inventory system

In order to let the supplier make more reasonable supply decisions,an integrated continuous replenishment policy for the vendor-managed inventory system is presented,which considers the quantity-based shipment consolidation and stock replenishment with lead time.Then the system cost is analyzed and a mathematical model is built.Since the model is rather complex,the bounds of the optimal policy are first attained,then the problem is solved by a heuristic algorithm.Through experiments the relationship between the order lead time and the corresponding integrated policy is discussed,and the influence on the system cost is also analyzed.The results reveal that the lead time's influence on the system is more serious with the increase of the order lead time,the integrated policy with the order lead time is more reasonable and the optimal policy can minimize the total system cost.Finally,the parameter sensitivity of the model is analyzed.

Optimal decision of retailer for replenishment cycle under a deteriorating product supply chain

In order to minimize the total cost of the retailer, an optimal replenishment cycle is studied by considering the deteriorating product, two-level trade credits, the limited storage capacity of their own warehouse and credit-linked order quantity simultaneously. A two-echelon supply chain model, which consists of a supplier and a retailer, is established. Then, the retailer＇s optimal replenishment cycle under all the cases are derived by using the optimization theory and method. On the basis of these, the effects of system parameters on the optimal replenishment cycle are examined by using the numerical studies. The results show that, when the retailer＇s trade credit period is longer （shorter） than the customer＇s trade credit period, the optimal replenishment cycle should he increased （decreased） as the retailer＇s trade credit period increases; if the minimum order quantity is high （low）, the optimal replenishment cycle should be increased （not changed） as the minimum order quantity increases.

Integrated continuous replenishment planning in VMI system of retail industry

In order to let suppliers monitor customers＇ inventory levels and make supply decisions regarding order quantities and delivery time, an integrated continuous replenishment planning （CRP） strategy of fresh food oriented to retail industry is presented, which integrates the consolidation of shipment and inventory replenishment considering the deterioration of items to guarantee the quality safety of fresh food during the replenishment cycle. A vendor managed inventory （VMI） mathematical model to compute upper-level inventory and delivery cycle is built. Based on the real-time sales data exchange, it enables suppliers to make the decision of the optimal time and the quantities of replenishment for retailers during a delivery cycle, in order to replenish the consumers＇ stock initiatively and minimize the long-run average cost.

Response of Soil Moisture to Rainfall Event in Black Locust Plantations at Different Stages of Restoration in Hilly-gully Area of the Loess Plateau, China

Precipitation plays an important role in the water supplies that support ecological restoration by sustaining large-scale artificial plantations in arid and semiarid regions, especially black locust(Robinia pseudoacacia) plantations(RP plantations), which are widely planted due to R. pseudoacacia being an excellent pioneer species. Characterizing the response of soil moisture to rainfall events at different stages of restoration is important for assessing the sustainability of restoration in RP plantations. In this study, we quantified the response of soil moisture to rainfall events at different years of restoration(15, 20 and 30 yr) representing different restoration stages in RP plantations in a typical hilly-gully area, i.e., the Yangjuangou Catchment, of the Loess Plateau, China. Over the growing season(June to September) of 2017, smart probes were placed at nine depths(10, 20, 40, 60, 80, 100, 120, 150, and 180 cm below the soil surface) to obtain volumetric soil water information at 30-min intervals in the three RP plantations. The advance of the wetting front was depicted, and the total cumulative water infiltration was measured. Soil moisture was mainly replenished by eight heavy rainfall events(mean rainfall amount = 46.3 mm), accounting for 88.7% of the rainfall during the growing season. The mean soil moisture content profiles of RP plantations at the three restoration stages were ordered as 30-yr(14.07%) > 20-yr(10.12%) > 15-yr(8.03%), and this relationship displayed temporal stability. Soil moisture was primarily replenished by rainfall at the 0-60 cm soil depth, and soil moisture remained stable below the 100-cm soil depth. The rainfall regime influenced the advancement of the wetting front. Here, a single rainfall event of 30 mm was the rainfall threshold for infiltration into the 60-cm soil layer. The total infiltration time ranged from 310.5-322.0 h, but no significant differences were found among RP plantations at different restoration stages. Young and old RP plantations had more total infiltration(more than 228.2 mm) and deeper infiltration depths(80-100 cm) than middle-aged plantations. The RP plantation at the intermediate restoration stage exhibited minimal total infiltration(174.2 mm) and a shallow infiltration depth(60 cm) due to the soil physical structure of the plot, which may have limited rain infiltration. More stand conditions that may affect infiltration should be considered for priority afforestation areas.

MILP formulations for highway petrol station replenishment in initiative distribution mode

To investigate highway petrol station replenishment in initiative distribution mode,this paper develops a mixed-integer linear programming(MILP)model with minimal operational costs that includes loading costs,unloading costs,transport costs and the costs caused by unpunctual distribution.Based on discrete representation,the working day is divided into equal time intervals,and the truck distribution process is decomposed into a pair of tasks including driving,standby,rest,loading and unloading.Each truck must execute one task during a single interval,and the currently executing task is closely related to the preceding and subsequent tasks.By accounting for predictive time-varying sales at petrol stations,real-time road congestion and a series of operational constraints,the proposed model produces the optimal truck dispatch,namely,a detailed task assignment for all trucks during each time interval.The model is tested on a real-world case of a replenishment system comprising eight highway petrol stations,one depot,one garage and eight trucks to demonstrate its applicability and accuracy.

Replenishment policy and inventory optimization for supply-hub with liability period consideration

A replenishment decision-making model for supply-hub is firstly established from the angle of supplier, and optimal replenishment decision of the supplier is analyzed. Then, inventory optimization model for supply-hub is formulated from the angle of the manufacturer, and the optimization algorithm for obtaining optimal inventory levels is given. The result shows that liability period decides the share of the inventory cost between two sides in supply chain. With the increase of liability period, the service level has been quickly reduced even though the manufacturer's cost has been cut down by transferring the inventory cost to the supplier. As to the safety inventory, if the lower bound of components safety inventory increases, the supplier's cost will rise up more slowly than the liability period does, while the service levels increases as the safety inventory's lower bound is raised.

Magma Replenishment as Revealed by Textural and Geochemical Features of Plagioclase Phenocrysts in Subduction-Related Lavas

Various petrographic features and geochemical characteristics indicative of disequilibrium are preserved in plagioclase phenocrysts from basaltic to andesitic lavas in East Junggar, northwest China. These characteristics indicate that they crystallized in a magma chamber, which was replenished by less differentiated and high-temperature magmas. The petrographic and geochemical features of the plagioclase phenocrysts are interpreted to record responses to changes in temperature, composition and mechanical effect during magma replenishment. Distinct rare earth element(REE) patterns between cores and rims of the same plagioclase crystal suggest derivation from two end-member magmas. From core to rim, plagioclase phenocrysts commonly display sharp fluctuations of anorthite(An) content up to 20, which either correspond to reverse zoning associated with ovoidal cores and resorption surface(PI), or normal zoning with euhedral form and no resorption surface(P2). Plagioclase crystals with diverse textures and remarkably different An content coexist on the scale of a thin-section. Cores of these plagioclases in each sample display a bimodal distribution of An content. From core to rim in PI, concentrations o f FeOT and Sr increase remarkably as An content increases. During magma replenishment, pre-existing plagioclase phenocrysts in the andesitic magma, which were immersed into hotter and less differentiated magmas, were heated and resorbed to form ovoidal cores, and then were overgrown by a thin rim with much higher contents of An, FeO^T and Sr. However, pre-existing plagioclase phenocrysts in the basaltic magma were injected into cooler and more evolved magmas, and were remained as euhedral cores, which were later enclosed by oscillatory zoned rims with much lower contents of An, Sr and Ba.

Key theoretical and technical issues and countermeasures for effective development of Gulong shale oil, Daqing Oilfield, NE China

Aiming at the four issues of underground storage state,exploitation mechanism,crude oil flow and efficient recovery,the key theoretical and technical issues and countermeasures for effective development of Gulong shale oil are put forward.Through key exploration and research on fluid occurrence,fluid phase change,exploitation mechanism,oil start-up mechanism,flow regime/pattern,exploitation mode and enhanced oil recovery(EOR)of shale reservoirs with different storage spaces,multi-scale occurrence states of shale oil and phase behavior of fluid in nano confined space were provided,the multi-phase,multi-scale flow mode and production mechanism with hydraulic fracture-shale bedding fracture-matrix infiltration as the core was clarified,and a multi-scale flow mathematical model and recoverable reserves evaluation method were preliminarily established.The feasibility of development mode with early energy replenishment and recovery factor of 3o%was discussed.Based on these,the researches of key theories and technologies for effective development of Gulong shale oil are proposed to focus on:(1)in-situ sampling and non-destructive testing of core and fluid;(2)high-temperature,high-pressure,nano-scale laboratory simulation experiment;(3)fusion of multi-scale multi-flow regime numerical simulation technology and large-scale application software;(4)waterless(CO_(2))fracturing technique and the fracturing technique for increasing the vertical fracture height;(5)early energy replenishment to enhance oil recovery;(6)lifecycle technical and economic evaluation.Moreover,a series of exploitation tests should be performed on site as soon as possible to verify the theoretical understanding,optimize the exploitation mode,form supporting technologies,and provide a generalizable development model,thereby supporting and guiding the effective development and production of Gulong shale oil.

Numerical study on wave loads and motions of two ships advancing in waves by using three-dimensional translating-pulsating source

A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.

Effects of improved fallow with Sesbania sesban on maize productivity and Striga hermonthica infestation in Western Kenya

Striga hermonthica is a major constraint to smallholder sub-sistence agriculture production in the sub-Saharan African region. Low soil fertility and overall environmental degradation has contributed to the build-up of the parasitic weed infestation. Improved cropping systems have to be introduced to address the interrelated problems of S. her-monthica and soil fertility decline. Thus, the effects of improved fallow with leguminous shrub Sesbania sesban on maize yields and levels of S. hermonthica infestation on farm land in the bimodal highlands of western Kenya were investigated. The experimental treatments were arranged in a phased entry, and randomized complete block scheme were six months Sesbania fallow, 18 months Sesbania fallow, six months natural fallow consisting of regrowth of natural vegetation without cultivation, 18 months natural fallow, continuous maize cropping without fertilizer application, and continuous maize cropping with P and N fertilization. Results show that Sesbania fallows significantly （p0.05） increase maize yield relative to continuous unfertilized maize. S. hermonthica plant populations decrease in continuous maize between the first season （mean = 428 000 ± 63 000 ha-1） and second season （mean=51 000 ± 15 000 ha-1）, presumably in response to good weed management. S. hermonthica seed populations in the soil decrease throughout the duration of the experi-ment in the continuous maize treatments. Short-duration Sesbania fal-lows can provide modest yield improvements relative to continuous unfertilized maize, but short-duration weedy fallows are ineffective. Continuous maize cultivation with good weed control may provide more effective S. hermonthica control than fallowing.

Data replenishment of five moderate earthquake sequences in Japan, with semi-automatic cluster selection

Missing early aftershocks following relatively large or moderate earthquakes can cause significant bias in the analysis of seismic catalogs.In this paper,we systematically address the aftershock missing problem for five earthquake sequences associated with moderate-size events that occurred inland Japan,by using a stochastic replenishing method.The method is based on the notion that if a point process(e.g.,earthquake sequence)with timeindependent marks(e.g.,magnitudes)is completely observed,it can be transformed into a homogeneous Poisson process by a bi-scale empirical transformation.We use the Japan Meteorological Agency(JMA)earthquake catalog to select the aftershock data and replenish the missing early events using the later complete part of each aftershock sequence.The time windows for each sequence span from 6 months before the mainshock to three months after.The semi-automatic spatial selection uses a clustering method for the epicentral selection of earthquakes.The results obtained for the original JMA catalog and replenished datasets are compared to get insight into the biases that the missing early aftershocks may cause on the Omori-Utsu law parameters’estimation,characterizing the aftershock decay with time from the mainshock.We have also compared the Omori-Utsu law parameter estimates for two datasets following the same mainshock;the first dataset is the replenished sequence,while the second dataset has been obtained by waveform-based analysis to detect early aftershocks that are not recorded in the JMA catalog.Our results demonstrate that the Omori-Utsu law parameters estimated for the replenished datasets are robust with respect to the threshold magnitude used for the analyzed datasets.Even when using aftershock time windows as short as three days,the replenished datasets provide stable Omori-Utsu law parameter estimations.The p-values for all the analyzed sequences are about 1.1 and c-values are significantly smaller compared to those of original datasets.Our findings prove that the replenishment method is a fast,reliable approach to address the missing aftershock problem.

Hydrodynamic interactions between two bodies in waves in 3D time domain

In this paper, a 3D time domain technique is adopted to calculate the coupled hydrodynamic interaction between two bodies without flare in waves. For verifying the code, two same cylinders are selected to calculate coupled hydrodynamic effects by comparison with the results obtained by 3D frequency method which has been proved to be efficient for solving such problems. In order to improve efficiency of calculation, the effect of history time has been discussed, and an improved method is presented. Moreover, the effect of lateral separation distance is also discussed in detail. The technique developed here may serve as a more rigorous tool to analyze the related transient problems of two ships doing underway replenishment in waves.

Development of a Portfolio Simulation System for Apparel Supply Chain

In this paper,a generic retailer-oriented portfolio simulation system linking manufacturer,retailer,and customer in apparel supply chain is presented.The purpose of the simulator is to generate a portfolio consisting of replenishment strategy and performance index under different sales forecasting errors to satisfy the retailer-defined customer service level in apparel industry.After analyzing the main parameters in the portfolio simulator,the procedure and detailed structure of the simulator are then described.With the use of data from the industry,one case study of the portfolio simulator is achieved and the process of the simulation is validated.

Maintaining eco-health of urban waterscapes with imbedded integrating ecological entity： Experimental approach

An imbedded integrating ecological entity(IIEE) was designed to combine landscaping, replenishing-water purifying and ecosystem maintaining simultaneously. With this IIEE, within 15 d experiment, simulated replenish water(SRW) with high(SRW-Ⅰ) or low(SRW-Ⅱ) nutrients concentration was well purified. Relative removal rates of CODCr, TP, TN, Chl-a and turbidity reached 84.87%, 84.05%, 94.76%, 188.17%, 110.93% when dealing SRW-Ⅰ, and 52.62%, 90.05%, 82.44%, 166.15%, 202.99%, respectively, when dealing SRW-Ⅱ. The well grew flora and fauna of IIEE benefit eco-maintaining and landscaping. Separately, the maximal root and stem length-increments of Cyperus alternifolius Linn. were 26.1 mm and 28.4 mm, while for Potamogeton crispus Linn. 18.3 mm and 25.7 mm. Mortality for both Bellamya aeruginosa and Misgurnus anguillicaudatus was both under 2.96%. The analysis of variance(ANOVA) indicated that most experimental indexes in each group performed more significantly better than those in their control. All results indicated that the IIEE is a promising technology for future urban waterscapes construction.

Textural and compositional zoning in plagioclase phenocrysts:implications for magma chamber processes in the Emeishan large Igneous Province,SW China

Textural and compositional zoning within plagioclase phenocrysts records the magma chamber processes,such as magma differentiation,magma recharge and mixing,and crustal contamination.The plagioclase phenocrysts in the Daqiao and Qiaojia plagioclase-phyric basalts from the Emeishan Large Igneous Province(LIP)show complex textural and compositional zoning patterns,e.g.,normal,reverse,oscillatory,and patchy zoning patterns.Most plagioclase phenocrysts exhibit a core–rim normal zoning pattern(Pl-A)with euhedral high-An cores(An=76–78%,in mole fraction)and low-An rims(An=68–72%),indicative of the crystal regrowth processes caused by recharge of relatively evolved magmas after the formation of high-An cores.Some phenocrysts have a core–rim reverse zoning pattern(Pl-B)with irregular ovaloid cores,characterized by extremely low An(60–61 mol%)and Ba(84–88 ppm)contents and extremely high87Sr/86Sr ratios(0.7120–0.7130).The rims of the Pl-B have relatively high An(69–72%),Ba(~160 ppm)contents,and low87Sr/86Sri(~0.7056).These Pl-B plagioclase phenocrysts preserve the information about the interaction between the crustal xenocrysts and the transporting magmas.Some plagioclase phenocrysts show a core–mantle–rim oscillatory zoning pattern(Pl-C)with multiple oscillations of An(70–80%),Ba(88–147ppm)from core to rim,revealing replenishment and mixing of multiple batches of basaltic melts with diverse compositions.87Sr/86Sr ratios of the Pl-C do not vary significantly(0.7050–0.7054).A small portion of phenocrysts has patchy patterns in the cores(Pl-D),where the low-An patches(72–75%)in form of elliptical or irregular elongated shapes were enclosed by the high-An domains(80–87%).These features can be attributed to crystal dissolution and regrowth processes during the reaction between earlyformed low-Cumulates and recharged hot primitive melts.The cores,mantles,and rims of different types of plagioclase phenocrysts(except the core of Pl-B)commonly display nearly constant Sr isotopic compositions,implying insignificant wall-rock assimilation at shallow-level magma reservoir(s)during the growth of these plagioclase phenocrysts.In conclusion,the massive crystallization of plagioclase in the late stage was an important controlling factor for the formation of iron-rich basalts in the Emeishan LIP.

Numerical Study on Motion Responses of Two Parallel Ships During Underway Replenishment by the Semi-Analytical HOTP Method

Motion responses of two ships advancing parallel in waves with hydrodynamic interactions are investigated in this paper. Within the framework of the frequency-domain potential flow theory, a semi-analytical higher-order translating-pulsating source(HOTP) method is presented to solve the problems of coupled radiation and diffraction potential. The method employs nine-node bi-quadratic curvilinear elements to discretize the boundary integral equations(BIEs) constructed over the mean wetted surface of the two ship hulls. In order to eliminate the numerical oscillation, analytical quadrature formulas are derived and adopted to evaluate the integrals related to the Froudedependent part of the Green’s function along the horizontal direction in the BIEs. Based on the method, a numerical program is originally coded. Through the calculations of hydrodynamic responses of single ships, the numerical implementation is proved successful. Then the validated program is applied in the investigations on the hydrodynamic interactions of two identical Wigley Ⅲ hulls and the underway replenishment of a frigate and a supply ship in waves with and without stagger, respectively. The comparison between the present computed results with experimental data and numerical solutions of other methods shows that the semi-analytical HOTP method is of higher accuracy than the pulsating source Green’s function method with speed correction and better stability than the traditional HOTP method based on Gauss quadrature. In addition, for two ships with obviously different dimensions,the influence of hydrodynamic interactions on the smaller ship is found to be more noticeable than that on the larger ship, which leads to the differences between the motions of frigate with and without the presence of supply ship.

Simulation-Based Optimization for the Fast Fashion Replenishment

Fast fashion is a commercial pattern which provides fashionable clothes at affordable price.This mode needs rapid response supply chain to respond to varying fashion trends.New styles are introduced in every sale period to cover fashion trends.In order to maximize profits,replenishment quantity of each style should be decided in every period.The purchasing and replenishing process over multiple periods based on uncertainty customer demand is modeled,which is formulated by a stochastic choice process.Heterogeneous consumers visit a store in a stochastic sequence and choosing dynamically from the available fashion styles(buy or not buy) according to a utility maximization criterion.The purchase process in a retail shop for multi-period is simulated.An algorithm which combines simulated anneal(SA) with gradient estimation is proposed to find the optimal replenishing strategy from the simulation program.

Design of cost allocation rule for joint replenishment with controllable lead time

To encourage retailers to form cooperative alliances to jointly replenish inventory,considering that the supplier provides a flexible lead time and quantity discount to retailers,a model of average total cost per unit time of periodic joint replenishment is constructed,and an approximate algorithm,which can satisfy the requirement of any given precision,is given.The cost allocation rule in the core of the joint replenishment game is designed based on the cooperative game theory.The numerical experiment results show that the proposed algorithm can quickly solve the joint replenishment problem when the item number is not greater than 640.The retailer's cost saving rate is always greater than 0,and it increases with the increase in quantity discount and fixed cost after adopting the given cost allocation rule.With the increase in the safety stock level,the retailer's cost saving rate increases first and then decreases;and the retailer's cost saving rate increases with the increase in the size of the alliance,but it decreases as the number of product category increases.The proposed cost allocation rule can reduce the retailer's cost up to 20%,which is conducive to forming a cooperative coalition.