Two-way Travel Time Research Articles

Along-strike variation in volcanic addition controlling post-breakup sedimentary infill: Pelotas margin, austral South Atlantic

Abstract. We investigate, using observations from seismic reflection data, the lateral variability in breakup extrusive magmatic addition along the strike of the Pelotas segment of the austral South Atlantic rifted margin and its control on post-rift accommodation space and sediment deposition. Our analysis of regional seismic reflection profiles shows that magmatic addition on the Pelotas margin varies substantially along strike from extremely magma-rich to magma-normal within a distance of ∼300 km. Using 2D flexural back-stripping, we determine the post-rift accommodation space above top volcanics. In the north, where volcanic seaward-dipping reflectors (SDRs) are thickest, the Torres High shows SDRs up to ∼20 km thick, and post-breakup water-loaded accommodation space is much less than in the south, where magmatic addition is normal and SDRs are thinner. We show that post-breakup accommodation space correlates inversely with SDR thickness, being less for magma-rich margins and more for magma-normal/intermediate margins. The Rio Grande Cone, with large sediment thickness, is underlain by small SDR thicknesses allowing large post-breakup accommodation space. A relationship is observed between the amount of volcanic material and the two-way travel time (TWTT) of first volcanics: first volcanics are observed between 1.2 and 2.2 s TWTT for the highly magmatic Torres High profile, while, in contrast, for the normally magmatic profiles in the south, first volcanics are observed between 4.2 and 6.5 s TWTT. The observed inverse relationship between post-breakup accommodation space and SDR thickness is consistent with predictions by a simple isostatic model of continental lithosphere thinning and magmatic addition melting during breakup. The methodology that we use in this paper provides a new approach for investigating the complex magmatic and sedimentary evolution of rifted continental margins.

Research on lunar regolith of the Chang'E-4 landing site: An automated analysis method based on deep learning framework

On January 3, 2019, the Chang'E-4 lander successfully landed within the Von Kármán crater, located in the South P ole-Aitken Basin (SPA) on the farside of the Moon (45.5°S, 177.6°E), marking the first soft landing on the lunar farside. The lander, equipped with the Lunar Penetrating Radar (LPR) system, aimed to provide insights into the structure and evolution of the Moon. Previous research often relied on manually identifying hyperbolic features to analyze the lunar shallow subsurface properties. This inefficient approach may lead to subjective biases, resulting in unstable outcomes. This research constructed an automatic analysis framework by integrating the Swin Transformer with a 3D velocity spectrum, which is then applied to analyze the properties of the Chang'E-4 LPR data. The experimental results indicate that the framework achieved a precision of 98.9 % and a recall of 96.7 % in hyperbolic feature identification, with an F1 of 0.9782 and AP of 94.8 %. Additionally, it has been experimentally validated that the framework can accurately invert hyperbolic features' two-way travel time and velocity. Finally, the framework is applied to analyze the lunar shallow subsurface structure and properties within the landing area of the Chang'E-4 mission.

Insights into the crustal and the magmatic feeding structure at the Payunia Volcanic Province highlighted by geophysical methods, in the retroarc of the Southern Central Andes

The Payunia Volcanic Province is a Quaternary volcanic plateau emplaced in the retroarc area in the northern Neuquén Mesozoic Basin, associated with a hydrocarbon system. At deeper levels, this basin is linked to different intrusive systems that developed in the retroarc region at different times, during the Jurassic, Cretaceous, Eocene, Miocene, and even the Quaternary which influenced the hydrocarbon system maturity. We analyzed the Moho structure through this retroarc region, as well as the crustal structure affected by different stages of regional extension. From continuous seismic noise data, we calculated the autocorrelograms to obtain the reflection response below each seismological station. This allowed imaging the surface of primary crustal reflectors and in a few stations the top of an asthenospheric anomaly (SWAP) found by magnetotelluric survey and in concordance with satellite magnetic data. The crustal reflectors were identified in all stations at a mean two-way travel time of about ∼8.5 s and ∼12.5 s using frequency bands of about 1.0–2.4 Hz. Therefore, this is the first geophysical research that estimates the depth of the magmatic system, hosted at the top of the lower crust and the Moho discontinuity. The deepest reflector, only recognized in 4 stations, was observed with a two-way travel time of 17.2 s to 19.6 s. We used a mean one-dimensional Vp model to obtain the corresponding reflector depths which constrain the two-dimensional forward gravity model that fits with the observed regional anomaly for the region. We finally established a relationship between the shallowest sublithospheric electrical conductivity anomalies determined in previous researches and the strong deep reflections observed in some of the seismological stations. This information may help to constrain geochemical and petrological models and re-evaluate the hydrocarbon system maturity of the northern Neuquén basin.

Anisotropy in shale and its impact on AVO modelling and prospect de-risk: a case study of a Central North Sea field

The inherent anisotropy in shale is an attribute associated with its clay minerals platelets alignment and layering and very common with transversely isotropic (TI) medium. This implies that velocity is fastest along the horizontal axis and slowest along the vertical axis, which could affect velocity signature from sonic tools measurements in deviated wells. This may lead to wrong amplitude variation with offset (AVO) signature, inaccurate depth conversions, seismic-to-well tie difficulty etc. Nineteen (19) wells in Forties field andthirteen (13) in Huntington field were investigated and interpreted for the presence of anisotropy effect within the caprock over the reservoir Forties sands of the Sele (shales) Formation. In order to analyse the Thomsen weak anisotropy parameters (ϵ, γ and δ) per field, wells with closer Sele depth range were considered. The well inclinations ( θ ), average compressional (α) and shear (β) velocities per well were utilised for forward and inverse modelling of Thomsen anisotropy equations. The estimated anisotropy (VTI) parameters and the P-wave ( α o ) and S-wave velocities (βo ) at zero inclinations in Sele Formation for the studied wells are: α o = 2448 ± 15.9 m/s, β o = 1003 ± 34.8 m/s, ϵ = 0.15 ± 0.02, γ = 0.33 ± 0.10, δ = 0.030 ± 0.015 in the Forties field Sele, while α o = 2775 ± 23.6 m/s, β o = 1306 ± 47.9 m/s, ϵ = 0.22 ± 0.02, γ = 0.43 ± 0.14, δ = 0.022 ± 0.018 in the Huntington field. AVO responses were compared, assuming both isotropic and anisotropic P–P reflectivity models. Result shows an increase in AVO gradient for all anisotropic case in both fields. AVO signatures were either completely or averagely masked in both fields using the isotropic assumptions. There is also an increase in AVO gradient for all anisotropic cases in both fields. Velocity decreases of about 1.2% and 3.5% were observed in the Sele Formation for the Forties and Huntington fields respectively for anisotropy. The corrected two-way traveltime (TWT) should in principle allow for a better seismic-to-well tie estimation. Similarly, using the faster deviated velocity (i.e. uncorrected velocity) in any depth conversion would result in under estimation of the actual depth across the field.

Detecting Reinforced Concrete Rebars Using Ground Penetrating Radars

A new algorithm is developed to automatically detect rebar locations and diameters of reinforced concrete structures using the ground penetrating radar technique. The study uses two-way travel time and biquadratic equations to formulate electromagnetic wave speed in reinforced concrete structures where hyperbolic signatures are approximated. Leveraging an established algorithm, a computer code has been developed to offer automated analysis of ground-penetrating radar data obtained from survey grids. Four reinforced concrete slabs were designed, fabricated, and tested to validate the developed evaluation approach. The proposed methodology demonstrates outstanding signal processing proficiency and reliably and effectively identifies rebar information.

Wykorzystanie metody radaru penetrującego grunt w badaniu gleby zanieczyszczonej węglowodorami w obszarze Navodari - Rumunia

Ground penetrating radar (GPR) is a very useful geophysical method for use in hydrogeologic and near-surface mapping studies. It can be used to study contaminants in groundwater, subsurface faulting, and underground cavities (natural or man-made), all of which pose potentially dangerous geological hazards. The GPR technique is similar in principle to seismic reflection and sonar techniques. The propagation of the radar signal depends on the frequency-dependent electrical properties of the ground.Electrical conductivity of the soil or rock materials along the propagation paths introduces significant absorptive losses which limit the depth of penetration into the earth formations and is primarily dependent upon the moisture content and mineralization present.Reflected signals are amplified, and transformed to the audio-frequency range, recorded, processed, and displayed. From the recorded display, subsurface features such as soil/ soil, soil/rock, and unsaturated/saturated interfaces can be identified. In addition, the presence of floating hydrocarbons on the water table, the geometry of contaminant plumes, and the location of buried cables, pipes, drums, and tanks can be detected. The GPR data are presented as a two-dimensional depth profile along a scanned traverse line in which the vertical axis is two-way travel time measured in nanoseconds. The location of hydrocarbon contamination in the ground using the GPR method is based mainly on information taken from reflected signals. In the cases investigated in Romania contaminated sites (Navodari area), such signals were very rarely recorded. A long time after spillage, contamination takes the form of plumes with different size and distribution, which depends on the geological and hydraulic properties of the ground. The survey discussed in this paper was carried out using the GPR system-Noggin with two antennas (250 and 500mHz) Data collected were processed using software(EKKO_Project™ GPR Data Analysis) to produce 2D radargram in time scale. The presence of contaminant plumes as well as the water table are observed in the GPR sections at depths approximately of 0.5 to 1.5 m. In the GPR section, the oil contaminated layer exhibits discontinuous, subparallel, and chaotic high amplitude reflection patterns. Promising results were also obtained in the GPR survey where three obvious reflection patterns representing the top sand-silt layer, oil-contaminated zone and, the underlying thick soft clay were detected in all 2D radargrams of the GPR traverse lines.

Seismic data extrapolation based on multi-scale dynamic time warping

Seismic data reconstruction can provide high-density sampling and regular input data for inversion and imaging, playing a crucial role in seismic data processing. In seismic data reconstruction, a common scenario involves a significant distance between the source and the first receiver, which makes it unattainable to acquire near-offset data. A new workflow for seismic data extrapolation is proposed to address this issue, which is based on a multi-scale dynamic time warping (MS-DTW) algorithm. MS-DTW can accurately calculate the time-shift between two time series and is a robust method for predicting time-offset (t−x) domain data. Using the time-shift calculated by the MS-DTW as the basic input, predict the two-way traveltime (TWT) of other traces based on the TWT of the reference trace. Perform autoregressive polynomial fitting on TWT and extrapolate TWT based on the fitted polynomial coefficients. Extract amplitude information from the TWT curve, fit the amplitude curve, and extrapolate the amplitude using polynomial coefficients. The proposed workflow does not necessitate data conversion to other domains and does not require prior knowledge of underground geological information. It applies to both isotropic and anisotropic media. The effectiveness of the workflow was verified through synthetic data and field data. The results show that compared with the method of predictive painting based on local slope, this approach can accurately predict missing near-offset seismic signals and demonstrates good robustness to noise.

Acoustic impedance prediction based on extended seismic attributes using multilayer perceptron, random forest, and extra tree regressor algorithms

Acoustic impedance is the product of the density of a material and the speed at which an acoustic wave travels through it. Understanding this relationship is essential because low acoustic impedance values are closely associated with high porosity, facilitating the accumulation of more hydrocarbons. In this study, we estimate the acoustic impedance based on nine different inputs of seismic attributes in addition to depth and two-way travel time using three supervised machine learning models, namely extra tree regression (ETR), random forest regression, and a multilayer perceptron regression algorithm using the scikit-learn library. Our results show that the R2 of multilayer perceptron regression is 0.85, which is close to what has been reported in recent studies. However, the ETR method outperformed those reported in the literature in terms of the mean absolute error, mean squared error, and root-mean-squared error. The novelty of this study lies in achieving more accurate predictions of acoustic impedance for exploration.

An adaptive MAC-agnostic ranging scheme for underwater acoustic mobile networks

In the last ten years several simulation studies on Autonomous Underwater Vehicle (AUV) swarm fleet formation have been performed, and some preliminary sea demonstrations of proof-of-concept prototypes were carried out. However, their actual realization is hindered by the challenges imposed by the underwater acoustic channel and the difficulties of keeping track of the vehicles’ positions due to the long latency required by traditional Two-Way Travel-Time (TWTT) ranging measurements, that require a specific signaling, hence limiting the throughput of the underwater network. Although One-Way Travel-Time (OWTT) halves the latency, it requires a high precision oscillator, such as an atomic clock or an oven controlled crystal oscillator (OCXO), to be installed in each modem processing unit: while atomic clocks are still very expensive, OCXO are very power demanding, making their application to underwater acoustic networks not always possible, especially in the case of low cost vehicle swarms.In this paper we present a network protocol stack able to perform ranging and localization within the communication task in underwater acoustic networks, limiting the network overhead. Specifically, new layers have been added to the preexisting DESERT Underwater protocol stack to perform the ranging tasks without compromising the correct operation of the communication network. A ranging layer is placed on top of the Medium Access Control (MAC) scheme, allowing the latter to be changed according to the network topology and requirements. This MAC-agnostic ranging layer is further optimized by adapting the amount of data transmitted according to the channel state, and the ranging entries inserted in the data packets according to the least recent information transmitted, hence minimizing the Age of Information. Simulation results obtained with the DESERT Underwater Framework show how this layer allows all AUVs in the swarm to know their distance from every other node in the network, thus limiting the probability of vehicle collisions and allowing better mission coordination.

Sediment thickness across Australia from passive seismic methods

SUMMARY Around the world the Earth’s crust is blanketed to various extents by sediment. For continental regions, knowledge of the distribution and thickness of sediments is crucial for a wide range of applications including seismic hazard, resource potential and our ability to constrain the deeper crustal geology. Excellent constraints on the sediment thickness can be obtained from borehole drilling or active seismic surveys. However, these approaches are expensive and impractical in remote continental interiors such as central Australia. Recently, a method for estimating the sediment thickness using passive seismic data, the collection of which is relatively simple and low cost, was developed and applied to seismic stations in South Australia. This method uses receiver functions, specifically the time delay of the P-to-S converted phase generated at the sediment-basement interface, relative to the direct-P arrival, to generate a first-order estimate of the thickness of sediments. In this work, we expand the analysis to the vast array of over 1500 seismic stations across Australia, covering an entire continent and numerous sedimentary basins that span the entire range from Precambrian to present day. We compare with an established yet separate method to estimate the sediment thickness, which utilizes the autocorrelation of the radial receiver functions to ascertain the two-way traveltime of shear waves reverberating in a sedimentary layer. Across the Australian continent the new results match the broad pattern of expected sedimentary features based on the various geological provinces. We are able to delineate the boundaries of many sedimentary basins, such as the Eucla and Murray Basins, which are Cenozoic, and the boundary between the Karumba Basin and the mineral rich Mount Isa Province. Contrasts in seismic delay time across these boundaries are upwards of 0.4 $\, \mathrm{s}$. The delay signal is found to diminish to <0.1 $\, \mathrm{s}$ for older Proterozoic basins, likely due to compaction and metamorphism of the sediments over time. As an application of the method, a comparison with measurements of sediment thickness from local boreholes allows for a straightforward predictive relationship between the delay time and the cover thickness to be defined. This offers future widespread potential, providing a simple and cheap way to characterize the sediment thickness in underexplored areas from passive seismic data.

The characteristics and distribution of the Moho reflections in the South China Sea Oceanic Basin

The distribution and characteristics of the Moho reflections (M-reflections) in oceanic basins record the seafloor spreading process and provide crucial insights into the formation and evolution of the oceanic basins. Approximately 20,000 km of multichannel seismic reflection sections are used to characterize the M-reflections and the crustal structures in the South China Sea oceanic subbasins. The identified M-reflections are mainly distributed in the northwest subbasin, the east subbasin, and the northern and southern flanks of the southwest subbasin (SWSB). They exhibit a banded distribution pattern in two-way traveltime along the north–south direction, with shallower depths observed in the central region and gradually deepening toward the northern and southern sides. In addition, the M-reflection amplitudes generally weaken, and the reflection events become less continuous as the igneous oceanic crust thins along with the seafloor spreading prior to the second ridge jump event. Ultimately, the M-reflection becomes indistinguishable in proximity to fossil spreading centers. These observations suggest a decreasing magma supply toward the cessation of the spreading. Despite the significance of Moho interface imaging at the crust-mantle boundary, the distinct M-reflections are absent in the SWSB and the vicinity of the Zhenbei-Huangyan Seamounts Chain, which are attributed to several factors, especially the vague nature of Moho interfaces.

Seismic data, cross-section balancing, and hydrocarbon prospectivity of the western Sulaiman fold belt, northwest Pakistan

Over the past few decades, hydrocarbon exploration has considerably improved our understanding of foreland deformation. This study used the interpretation of seismic reflection profiles and cross-section balancing techniques to understand the passive-roof duplex geometry of foreland structures in the Western Sulaiman Fold Belt (SFB), Pakistan. Seismic reflection profiles were interpreted to show the presence of hinterland-vergent fault-propagation and fault-bend folds as prospects over a detachment located at a depth of ∼6 km in the pelitic Triassic succession. Cross-section balancing showed intraplate deformation in the hanging-wall flat of the frontal duplex. Clear seismic reflectors above 3 s two-way travel time (TWT), estimates of stratigraphic and tectonic thicknesses with field geology and isopach mapping were used to locate the décollement in the pelitic Permo-Triassic strata at a depth of ∼11 km and interpret stacked duplexes at the mountain front. The foreland deformation of the western SFB was compared with that of the central/eastern SFB. The comparison shows the presence of discovered gas fields, with intraplate deformation of the duplex in the western, and pop-ups in the roof-sequence in eastern SFB depicting the complexity of deformation for the assessment of hydrocarbon prospectivity in the fold belt. The 33 km deformed balanced cross-section was restored to ∼71 km, with 38 km (54%) shortening in the foreland. Approximately 12 km (17%) of shortening was passively accommodated by the hinterland-vergent out-of-sequence passive backthrust in the frontal horse. This description of foreland deformation using the example of the SFB adds to our understanding of the structures of the mountain front that are recognized as zones of prolific hydrocarbon prospects.

2D Seismic Structural Interpretation of Yamama Formation in Al-Fao Area, Basrah, Iraq

A structural study was conducted in Al-Fao area, southeastern Iraq using two-dimensional seismic data. The seismic lines were divided into three surveys (F, Kl, and Sa). Al-Fao is an exploratory area that does not contain any wells; the well information was based on Siba field, located north of the study area. The primary reservoir target in the lower Cretaceous was Yamama Formation. Two reflectors,top and base Yamama Formation are defined using synthetic seismogram in the time domain for Sb well. A top Yamama Formation was defined as Z-crossing, and a base Yamama Formation as S-crossing. The results showed that the strata tended to slope from the south and southwest to the north and northeast, and the area was affected by multiple structural orogenies. Siba field is an elongated fold whose axis extends northeast-southwest; it consists of a structural nose and northeastern dome. The Two-Way Travel Time map showed higher values in the southeast and east, indicating sloping reflectors in those directions. The depth maps of Yamama Formation showed the deepest part in the southeast and east, suggesting a plunging structure in those directions. The study also revealed the presence of numerous closures as structural noses, which could serve as potential structural traps.

Three-dimensional seismic evidence for depositional undulations nucleated around pre-existing relief, South China Sea

Sediment undulations are widely distributed along continental shelf and slope worldwide. The driving mechanisms of undulations are crucial for assessing the long-term stability of continental slope, which, however, remain poorly understood. We used high-resolution 3D seismic data to investigate the internal architecture and origin of a field of seafloor and subsurface undulations along the slope of South China Sea (SCS). The undulations have narrow troughs and broad ridges. The reflection geometries in the troughs between undulations tend to be continuous and not offset by faults. These troughs terminate at different subsurface depths, commonly ranging from 50 to 400 ms in two-way travel time (TWTT) below the seafloor, often at locations related to subsurface geological structures, such as buried faults and steps within mass transport deposits (MTDs). These subsurface terminations therefore suggest that the undulations are probably depositional features not linked to creep-related basal detachment. The troughs between undulations display an absence of growth features, while the undulations increase in amplitude from base to top (characteristics that are also incompatible with gravitational deformation). Erosional and depositional features are commonly found in the undulations; thus, we interpret the narrow troughs between undulations to represent the areas where sediment-laden bottom currents undergo passage from supercritical to subcritical conditions. These troughs nucleated mostly from fault scarps and seafloor steps induced by past slope failures. Previous alternative interpretations report that these features originate from subsurface creep; alternatively, the present interpretation as sediment transport structures has significant implications for seafloor management and engineering projects.

Prediction of acoustic impedance and density porosity using seismic inversion, and geostatistical methods on Krishna Godavari basin, India: A case study

Seismic inversion is a geophysical technique used to estimate subsurface acoustic impedance models from seismic reflection data. The basic concept behind seismic inversion is to establish mathematical assumptions that can relate seismic responses to geological formation properties. This study addresses two types of seismic inversion methods, namely Band-limited inversion (BLI) and Linear programming inversion (LPI) techniques for the delineation of gas reservoirs. The study presents the successful use of the inversion in a real example from a gas sand reservoir in KG-Basin, India. Inversion has led to the removal of ambiguity and revealing clear information about the target area. This work aims to use and compare the results of these two generic seismic post-stack inversion methods for characterizing the reservoir. The inversion reveals the presence of low acoustic impedance around 2320-2440 ms and is inferred as a sandstone layer with a potential hydrocarbon reservoir. One can improve the analysis of the elastic properties for reservoir characterization by using rock physics knowledge to understand better reservoir properties such as water saturation, porosity, and shale volume. This study also uses the Rock Physics Template (RPT) and plots it with elastic properties to aid the efficient interpretation. RPT helps to identify the lithology and fluid content in the area. Results from the seismic study suggest a comparatively better analysis done through the LPI technique over BLI. In the last, the Multi-attribute transform as a part of the Geostatistical method is used to predict porosity away from the borehole. The method shows a very high porosity section and clearly shows the distribution of porosity in and around the region. The interpretation shows that a very high porosity zone (15–25%) is present that is in cross-correlation with a low impedance zone (4000-8000 m/s*g/cc) and confirms the presence of a hydrocarbon zone in-between 2320 and 2440 ms two-way travel time. All seismic inversion and geostatistical methods can retrieve subsurface information with very high-resolution images but the analysis based on linear programming is found slightly better in each step and depicts good prediction results. The study finally gives a flowchart that combines seismic inversion and geostatistical methods to extract petrophysical parameters from seismic data and can help to interpret hydrocarbon-bearing formation in any virgin area.

Excavation and Construction Technology of Diversion Tunnel under Complex Geological Conditions

During the construction of a diversion tunnel, geological problems often include faults, fragile strata, hard rock formations, karst landforms, etc., which may have adverse effects on the excavation and construction of the diversion tunnel. Based on the analysis of the engineering overview, this study designed a new construction technology for the excavation of water diversion tunnels in hard rock layers and high-karst areas. Based on tunnel seismic prediction (TSP) technology to achieve advanced geological prediction, combined with actual geological conditions, construction difficulties are analyzed. Then, detection technology is used to collect two-way travel time, amplitude, and waveform data. By processing and analyzing the detection image, the spatial orientation and length of the main tunnel during the construction of the diversion tunnel are calculated. After completing the construction ventilation and wind, water, and electricity layout, the excavation construction procedure is designed. In the specific excavation design, the tunnel curtain excavation technology, tunnel body excavation and support technology, excavation grouting technology, important unfavorable geological tunnel section excavation technology, upper/lower flat section excavation technology, upper/lower curved section excavation technology, and vertical shaft section excavation technology were elaborated. Finally, a plan was made for the reuse process of slag material and a construction quality control system was established. During the testing process, it was found that the antidamage coefficient of the side wall was above 0.9 after using the technology described in this article. Therefore, it indicates that the excavation construction technology designed in this article can ensure the support capacity of the side wall of the diversion tunnel, which is suitable for the excavation of the main tunnel during the construction of the diversion tunnel.

Methane-flow system within the Nyegga pockmark field, offshore mid-Norway

We investigated fluid seepage within the Nyegga pockmark field (600–900 m water depths) off mid-Norway from Remotely Operated Vehicle dives at the so-called CNE sites (CNE01 to CNE17). The seafloor morphology of some of these features corresponds to pockmarks and adjacent ridges, with the latter being the focus of present seepage activity. These structures are underlain by chimneys above a gas-charged zone with, in some cases, a substantial body of hydrate-invaded sediment (down to 1.3 s in two-way travel time at CNE03). Present-day methane-rich fluid seepage through the seabed is indicated by chemosynthetic fauna, in particular Siboglinidae polychaetes (Oligobrachia haakonmobiensis webbi and Sclerolinum contortum), microbial mats and associated Rissoidae gastropod (Alvania sp.) grazers, and confirmed by measured in situ bottom-water methane anomalies, up to 2,130 nL/L. No free-gas bubble emissions were observed or acoustically identified. The presence of authigenic carbonates reveals past seepage with very low δ13C values (down to −58‰) indicating that the major source of carbon was methane carried by the venting fluids. The ages of major periods of methane venting are provided by vesicomyid bivalve shells (Isorropodon nyeggaensis) present in two sedimentary layers, 14,930 and 15,500 14C yr BP (ca. 17,238 and 17,952 cal yr BP), respectively, corresponding to the time of Melt Water Pulse IA. The seafloor morphology and pattern of seepage -chemosynthetic fauna and microbial mat distribution and dissolved methane concentration-are remarkably heterogeneous. Pore-water chemistry profiles in a gravity core taken only 40 m from major seepage sites indicate no seepage and anaerobic methane oxidation at a sub-bottom depth of about 2 m. Present-day seepage from the studied pockmark-chimney fluid-flow system charged with gas hydrate is dominated by the advection of methane solution in pore water. Some of this methane could result from the dissolution of hydrate in the chimney, most of which would have formed during an earlier period (post-LGM times) of history of the chimney, when it was venting free gas. However, the presence of free gas beneath this chimney is probably why the water entering the chimney is already saturated with methane and the process of hydrate formation in the chimney continues today.

Underwater Localization of AUVs in Motion Using Two-Way Travel Time Measurements With Unknown Sound Velocity

For an autonomous underwater vehicle (AUV) underwater localization system based on acoustic two-way travel time measurements, the motion of the AUV during the interrogation-reception time interval, unknown underwater sound velocity, and beacon calibration error affect the performance of the positioning solution. Ignoring these factors significantly deteriorates the performance of the localization system. This paper focuses on the localization problem of AUV in motion with an unknown sound velocity and beacon calibration error. The Cramér–Rao lower bound of the concerned localization problem is first derived as the performance bound of an unbiased estimator. Next, through multi-stage processing, we design a closed-form solution that can attain the corresponding performance bound under small measurement noise conditions. Furthermore, considering that the AUV motion velocity is typically much lower than the sound velocity in practical applications, a simpler and more efficient simplified closed-form solution is proposed, which can guarantee localization performance under the condition that the AUV motion speed is low, as well as small measurement noise. Finally, the theoretical performance and superiority of the proposed solutions were verified through numerical simulation.

Field study of the dielectric constant of concrete: A parameter less sensitive to environmental variations than electrical resistivity

In this work, a field study was conducted to evaluate the curing effectiveness of concrete slabs cured by air and compound. Two non-destructive testing methods including ground penetrating radar (GPR) and electrical resistivity measurement were used. The dielectric constant value developments of air and compound cured concrete slabs were determined by GPR measurement. The compressive strength development was tested by concrete cylinders. The hydration degrees of air and compound cured slabs were determined by Ca(OH)2 (CH) content through thermogravimetric analysis (TGA) by a modified calculation method. The microstructure development of air and compound cured concrete slabs were characterized by mercury intrusion porosimetry measurement (MIP). It has found that the dielectric constant value of compound cured concrete slab was always higher than air cured concrete slab. The second reflected pulse peak amplitude of GPR waveform was less influenced by weather conditions, making it a reliable indicator of the dielectric constant value. On the other hand, the dielectric constant value determined by the first reflected pulse peak amplitude was greatly influenced by the weather conditions. The dielectric constant value determined by the two-way travel time method exhibited strong correlations with compressive strength, CH content and porosity. However, the electrical resistivity development of air and compound cured concrete slabs were not consistent with compressive strength, CH content and microstructure development. Our findings suggested that the GPR measurement can be practically used to evaluate the curing effectiveness of concrete and can be a good indicator of “hydration potential” from a development perspective.

Automatic quantitative recognition method for vertical concealed cracks in asphalt pavement based on feature pixel points and 3D reconstructions

Non-destructive detection and characterisation of concealed cracks in asphalt pavements is a key aspect of ground-penetrating radar (GPR) road maintenance. This paper proposes a method based on feature pixel points to quantify and calculate the vertical height of cracks. Firstly, the team used gprMax to conduct numerical simulations to study the GPR image characteristics of concealed cracks in asphalt pavement with varying lengths and widths. Subsequently, the relationship between the pixel value of the crack area and the two-way travel time was established to obtain the relationship between the vertical height of the crack and the pixel. This study combined with the deep learning model (YOLOv5) allows for the calculation of the vertical height of a crack while simultaneously recognizing it, with the minimum error being only 1.3 %. Finally, concealed cracks were visualized in three-dimensional(3D) by slicing the results of vertical crack simulations of asphalt pavement and observing cloud images, and the principle of computed tomography(CT) was employed to reconstruct 3D models of cracked asphalt pavement and estimate the vertical height of the cracks. This method achieved a minimum error of only 2.9 %. The research presents a theoretical framework for recognizing and more intuitive characterizing concealed cracks in asphalt pavement accurately, enabling it to be used in practical engineering applications.