Geomorphology, sedimentary processes and development of the Zenisu deep-sea channel, northern Philippine Sea

The Zenisu deep-sea channel originates on the Izu-Ogasawara island arc, and disappears in the Shikoku Basin of the Philippine Sea. The geomorphology, sedimentary processes, and the development of the Zenisu deep-sea channel were investigated on the basis of swath bathymetry, side-scan sonar imagery, submersible observations, and seismic data. The deep-sea channel can be divided into three segments according to the downslope gradient and channel orientation. They are the Zenisu Canyon, the E-W fan channel, and the trough-axis channel. The sediment fill is characterized by turbidite and debrite deposition and blocky-hummocky avalanche deposits on the flanks of the Zenisu Ridge. In the Zenisu Canyon and the Zenisu deep-sea channel, sediment transport by turbidity currents generates sediment waves (dunes) observed during the Shinkai 6500 dive 371. The development of the Zenisu Canyon is controlled by a N-S shear fault, whereas the trough-axis channel is controlled by basin subsidence associated with the Zenisu Ridge. The E-W fan channel was probably affected by the E-W fault and the basement morphology.

Sedimentary processes and development of the Zenisu deep-sea channel, Philippine Sea

Zenisu deep-sea channel originated from a volcanic arc region, Izu-Ogasawara Island Arc, and vanished in the Shikoku Basin of the Philippine Sea. According to the swath bathymetry, the deep-sea channel can be divided into three,segments. They are Zenisu canyon, E-W fan channel and trough-axis channel. A lot of volcanic detritus were deposited in the Zenisu Trough via the deep-sea channel because it originated from volcanic arc settings. On the basis of the swath bathymetry, submersible and seismic reflection data, the deposits are characterized by turbidite and debrite deposits as those in the other major deep-sea channels. Erosion or few sediments were observed in the Zenisu canyon, whereas a lot of turbidites and debrites occurred in the E-W channel and trough axis channel. Cold seep communities, active fault and fluid flow were discovered along the lower slope of the Zenisu Ridge. Vertical sedimentary sequences in the Zenisu Trough consist of the four post-rift sequence units of the Shikoku Basin, among which Units A and B are two turbidite units. The development of Zenisu canyon is controlled by the N-S shear fault, the E-W fan channel is related to the E-W shear fault, and the trough-axis channel is related to the subsidence of central basin.

Suspended particulate matter fluxes through the Straits of Dover, English Channel: observations and modelling

Suspended Particulate Matter (SPM) concentrations at various levels within the water column, together with salinity and temperature, were measured using water samples collected from six stations across the Straits of Dover. The sampling programme covered a 16-month period, undertaken during 23 cruises. On the basis of the spatial variability in the concentrations, the water bodies are divided by several boundaries, controlled by tidal and wind conditions. Within the water column, SPM concentrations were higher near the sea bed than in the surface waters. Throughout the cross-section, maximum concentrations occurred adjacent to the coastlines. Temporal variability in the SPM concentration exists on daily and seasonal scales within the coastal waters (4.2 to 74.5 mg L-1): resuspension processes, in response to semi-diurnal tidal cycles (with a period of around 12.4 h) and spring-neap cycles (with a period of 15 days) make significant contributions. Distinctive seasonal/annual concentration changes have also been observed. In the offshore waters, such variability is much less significant (0.9 to 6.0 mg L-1). In the summer the English Coastal Zone is associated with relatively high SPM concentrations: the Central Zone has a low and stable SPM concentration between these zones, there is a Transitional Zone, where there is a rapid response of SPM concentration to wind forcing. Finally, the French Coastal Zone is characterized by variable (sometimes high) SPM concentrations. Because of the zonation, SPM fluxes within the Dover Strait are controlled by different transport mechanisms. Within the Central Zone, the flux can be represented by the product of mean water discharges and SPM concentrations. However, within the coastal zones fluctuations in SPM concentrations on various time-scales must be considered. In order to calculate the maximum and minimum SPM fluxes, 10 cells were divided in the strait. A simple modelling calculation has been proposed for this complex area. The effect of spring-neap tidal cycles and seasonal changes can contribute significantly to the overall flux, which is of the order of 20 x 10(6) t.yr(-1) (through the Dover Strait, towards the North Sea). Such an estimate is higher than most obtained previously. (C) 2000 Ifremer/CNRS/IRD/Editions scientifiques et medicales Elsevier SAS.

Resuspension and advection processes affecting suspended particulate matter concentrations in the central English Channel

Suspended particulate matter (SPM) measurements obtained along a cross-section in the central English Channel (Wight-Cotentin transect) indicate that the area may be differentiated into: (1) an English coastal zone, associated with the highest concentrations; (2) a French coastal zone, with intermediate concentrations; and (3) the offshore waters of the Channel, characterised by a very low suspended-sediment load. The SPM particle-size distribution was modal close to the English coast (main mode 10-12 mu m); the remainder of the area was characterised by flat SPM distributions. Examination of the diatom communities in the SPM suggest:; that material resuspended in the intertidal zone and the estuarine environments was advected towards the offshore waters of the English Channel. Considerable variations in SPM concentrations occurred during a tidal cycle: maximum concentrations were sometimes up to 3 times higher than the minimum concentrations, Empirical orthogonal function (EOF) analysis of the SPM concentration time series indicates that, although the bottom waters were more turbid than the surficial waters, this was not likely to be the result of in situ sediment resuspension. Instead, the observed variations appear to be controlled mainly by advective mechanisms. The limited resuspension was probably caused by: (1) the limited availability of fine-grained material within the bottom sediments, and (2) 'bed-armouring' processes which protect the finer-grained fractions of the seabed material from erosion and entrainment within the overlying flow during the less energetic stages of the tide.

River mouth bar formation, riverbed aggradation and channel migration in the modern Huanghe (Yellow) River delta, China

This paper addresses the recent (1970s-1990s) processes of river mouth bar formation, riverbed aggradation and distributary migration in the Huanghe River mouth area, in the light of station-based monitoring, field measurements and remote sensing interpretation. The results show that the morphological changes of the river mouth bar have been closely associated with the largely reduced fluvial discharge and sediment load. Landforrn development such as bar progradation occurred in two phases, i.e. before and after 1989, which correspond to faster and lower bar growth rates, respectively. Fast riverbed aggradation in the mouth channel was strongly related to river mouth bar progradation. During 1976-1996, about 2.8% of the total sediment loads were deposited in the river channel on the upper to middle delta. Therefore, the river water level rose by a few meters from 1984 to 1996. The frequent distributary channel migration, which switched the radial channel pattern into the SE-directed pattern in the mid-1980s, was linked with mouth bar formation. Marine conditions also constrain seaward bar progradation. Furthermore, the history of river mouth bar formation reflects human impacts, such as dredging and dyking in order to stabilize the coastal area. (c) 2005 Elsevier B.V. All rights reserved.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.