A new efficient 3D non-hydrostatic free-surface flow model for simulating water wave motions

A new three-dimensional, non-hydrostatic free surface flow model is presented. For simulating water wave motions over uneven bottoms, the model employs an explicit project method on a Cartesian the staggered gird system to solve the complete three-dimensional Navier–Stokes equations. A bi-conjugated gradient method with a pre-conditioning procedure is used to solve the resulting matrix system. The model is capable of resolving non-hydrostatic pressure by incorporating the integral method of the top-layer pressure treatment, and predicting wave propagation and interaction over irregular bottom by including a partial bottom-cell treatment. Four examples of surface wave propagation are used to demonstrate the capability of the model. Using a small of vertical layers (e.g. 2–3 layers), it is shown that the model could effectively and accurately resolve wave shoaling, non-linearity, dispersion, fission, refraction, and diffraction phenomena.     

4D seismic time-lapse monitoring of an active cold vent,northern Cascadia margin

Two single-channel seismic (SCS) data sets collected in 2000 and 2005 were used for a four-dimensional (4D) time-lapse analysis of an active cold vent (Bullseye Vent). The data set acquired in 2000 serves as a reference in the applied processing sequence. The 4D processing sequence utilizes time- and phase-matching, gain adjustments and shaping filters to transform the 2005 data set so that it is most comparable to the conditions under which the 2000 data were acquired. The cold vent is characterized by seismic blanking, which is a result of the presence of gas hydrate in the subsurface either within coarser-grained turbidite sands or in fractures, as well as free gas trapped in these fracture systems. The area of blanking was defined using the seismic attributes instantaneous amplitude and similarity. Several areas were identified where blanking was reduced in 2005 relative to 2000. But most of the centre of Bullseye Vent and the area around it were seen to be characterized by intensified blanking in 2005. Tracing these areas of intensified blanking through the three-dimensional (3D) seismic volume defined several apparent new flow pathways that were not seen in the 2000 data, which are interpreted as newly generated fractures/faults for upward fluid migration. Intensified blanking is interpreted as a result of new formation of gas hydrate in the subsurface along new fracture pathways. Areas with reduced blanking may be zones where formerly plugged fractures that had trapped some free gas may have been opened and free gas was liberated.     

Application of a two-dimensional depth-averaged hydrodynamic tidal model

A two-dimensional hydrodynamic model application to the San Francisco Bay was performed using the Boundary-Fitted HYDROdynamic model (BFHYDRO). The model forcing functions consist of tidal elevations along the open boundary and fresh water flows from the Delta Outflow. The model-predicted surface elevations compare well with the observed surface elevations at five stations in San Francisco Bay. Mean error in the model predicted surface elevations and currents are less than 7 and 9%, respectively. Correlation coefficients for surface elevations and currents are higher than 0.94 and 0.95, respectively. The amplitudes and phases of the principal tidal constituents at 24 tidal stations in San Francisco Bay, obtained from a harmonic analysis of a 90-day simulation compare well with the observed data. The predicted amplitude and phase of the M₂ tidal constituent at these stations are respectively within 8 cm and 8° of the observed data. Maximum errors in the K₁ harmonic amplitudes and phases are less than 3 cm and 7° respectively. The asymmetric diurnal and semi-diurnal tidal ranges and spring and neap tidal cycles of the surface elevations and currents are well reproduced in the model at all stations.     

Fast static correction methods for high-frequency multichannel marine seismic reflection data: A high-resolution seismic study of channel-levee systems on the Bengal Fan

Very high-frequency marine multichannel seismic reflection data generated by small-volume air- or waterguns allow detailed, high-resolution studies of sedimentary structures of the order of one to few metres wavelength. The high-frequency content, however, requires (1) a very exact knowledge of the source and receiver positions, and (2) the development of data processing methods which take this exact geometry into account. Static corrections are crucial for the quality of very high-frequency stacked data because static shifts caused by variations of the source and streamer depths are of the order of half to one dominant wavelength, so that they can lead to destructive interference during stacking of CDP sorted traces. As common surface-consistent residual static correction methods developed for land seismic data require fixed shot and receiver locations two simple and fast techniques have been developed for marine seismic data with moving sources and receivers to correct such static shifts. The first method – called CDP static correction method – is based on a simultaneous recording of Parasound sediment echosounder and multichannel seismic reflection data. It compares the depth information derived from the first arrivals of both data sets to calculate static correction time shifts for each seismic channel relative to the Parasound water depths. The second method – called average static correction method – utilises the fact that the streamer depth is mainly controlled by bird units, which keep the streamer in a predefined depth at certain increments but do not prevent the streamer from being slightly buoyant in-between. In case of calm weather conditions these streamer bendings mainly contribute to the overall static time shifts, whereas depth variations of the source are negligible. Hence, mean static correction time shifts are calculated for each channel by averaging the depth values determined at each geophone group position for several subsequent shots. Application of both methods to data of a high-resolution seismic survey of channel-levee systems on the Bengal Fan shows that the quality of the stacked section can be improved significantly compared to stacking results achieved without preceding static corrections. The optimised records show sedimentary features in great detail, that are not visible without static corrections. Limitations only result from the sea floor topography. The CDP static correction method generally provides more coherent reflections than the average static correction method but can only be applied in areas with rather flat sea floor, where no diffraction hyperbolae occur. In contrast, the average static correction method can also be used in regions with rough morphology, but the coherency of reflections is slightly reduced compared to the results of the CDP static correction method.     

悬浮沙激光散斑图像处理的研究

本文以统计模式识别为基础，阐述了图像处理在悬浮沙粒粒径流量中的实际应用，提出了用图像层析来确定放大系数的方法，提出了测量的准确度。     

A Sensitivity Analysis of the Waterline Method of Constructing a Digital Elevation Model for Intertidal Areas in ERS SAR scene of Eastern England

A sensitivity analysis of the waterline method of constructing a Digital Elevation Model (DEM) of an intertidal zone using remote sensing and hydrodynamic modelling is described. Variation in vertical height accuracy as a function of beach slope is investigated using a set of nine ERS Synthetic Aperture Radar (SAR) images of the Humber/Wash area on the English east coast acquired between 1992 and 1994. Waterlines from these images are heighted using a hydrodynamic tide-surge model and interpolated using block kriging. On 1:500 slope beaches, an average block height estimation standard deviation of 18–22 cm is achieved. This rises to 27 cm on 1:100 slope beaches, and 32 cm on 1:30 slope beaches. The average heighting error at different slopes is decomposed into components due to waterline heighting error, inadequate sensor resolution and interpolation inaccuracy. It is shown that, at 1:500 slope, waterline heighting error and interpolation inaccuracy are the main error sources, whilst at 1:30 slope, errors due to inadequate sensor resolution become dominant. The ability of the technique to generate intertidal DEMs for almost the entire coastal zone in a complete ERS SAR scene covering 100×100 km is demonstrated.     

走航测深数据的自动采集

本文介绍应用IBM微机实现走航实时水深及定位数据的自动采集系统。给出了联机系统方框图及采集、传送数据的程序流程图。并对IBM-PC/XT微机与测深仪之间的接口、数码转换和数据处理作了简要的叙述。     

潜艇疲劳载荷的概率模型

探讨了建立潜艇疲劳载荷概率模型的方法，引入了正态分布和两参数威布尔分布两种概率模型。两参数威布尔分布较适合于用来描述潜艇下潜深度分布的概率特征。可根据潜艇的设计参数、任务及航行区域等因素来选取最可能的分布形式获得潜深分布密度函数。潜艇疲劳热点部位的应力(应变)幅值分布，可由应力(应变)和潜深的关系通过相应的变换得到。