後藤 忠徳

Blurred resistivity boundaries resulting from smoothness-regularized inversions of electrical resistivity tomography (ERT) data can lead to inaccurate interpretations of sharp boundary structures. To address this issue, various ERT inversion algorithms have introduced localized adjustments (localized discontinuities) in the regularization operator at positions where sharp boundaries are anticipated. Current approaches rely on prior information about sharp boundary locations, obtained from complementary geophysical, geological, and drilling data, to determine the positions and weights for these regularization adjustments. However, such prior information is frequently insufficient, limiting the application of localized regularization adjustments. Accordingly, we developed a sharp boundary inversion (SBI) algorithm using the Akaike Bayesian Information Criterion (ABIC) that determines the optimal positions and weights for localized regularization adjustments by testing various configurations and selecting the one that minimizes ABIC. A synthetic modeling study demonstrated that the SBI algorithm correctly delineated the sharp boundaries of a conductor. Its application to field data demonstrated that it delineated the sharp boundaries of a utility tunnel, and the size and horizontal position of the recovered tunnel were consistent with the estimated dimensions from the blueprint. As it does not rely heavily on prior information, the SBI algorithm can be applied to a wide range of geophysical survey data, even when prior knowledge of sharp boundary locations is limited.

Abstract Near-seafloor concentrated gas hydrates (GHs) containing large amounts of methane have been identified at various gas chimney sites. Although understanding the spatial distribution of GHs is fundamental for assessing their dissociation impact on aggravating global warming and resource potential, the spatial distribution of GHs within gas chimneys remains unclear. Here, we estimate the subseafloor distribution of GHs at a gas chimney site in the Japan Sea using marine electrical resistivity tomography data. The resulting two-dimensional subseafloor resistivity structure shows high anomalies (10–100 Ωm) within seismically inferred gas chimneys. As the resistivity anomalies are aligned with high amplitude seismic reflections and core positions recovering GHs, we interpret the resistivity anomalies are near-seafloor concentrated GH deposits. We also detect various distribution patterns of the high resistivity anomalies including 100-m wide and 40-m thick anomaly near the seafloor and 500-m wide anomaly buried 50 m below the seafloor, suggesting that GHs are heterogeneously distributed. Therefore, considering such heterogeneous GH distribution within gas chimneys is critical for in-depth assessments of GH environmental impacts and energy resources.

Abstract The physical properties of seafloor massive sulfides are crucial for interpreting sub-seafloor images from geophysical surveys, shedding light on the evolution of seafloor mineral deposits. While some studies have explored the relationship between electrical properties and the volume of conductive minerals in rocks from seafloor massive sulfide deposits, they primarily focused on artificial samples, leaving the characteristics of natural samples less understood. Moreover, there has been no comprehensive study detailing the general characteristics of electrical properties, particularly chargeability and relaxation time, in relation to the volumetric fraction of sulfides in rocks from massive sulfide mounds in typical hydrothermal areas. In this study, we employed complex conductivity measurements, elemental concentration analysis, and mineral content identification on to rock samples from the active hydrothermal zones of the Okinawa Trough in Japan. The complex conductivity observed was remarkably high, with a pronounced imaginary component and a broad frequency range. This is attributed to induced polarization extending beyond our measurement range. The rock samples were rich in conductive sulfide minerals such as pyrite, chalcopyrite, and galena. Using the Cole–Cole rock physics model, we established a correlation between rock chargeability and relaxation time coefficient with the volume fraction of conductive sulfide minerals, which deviated from previous findings. The intensity of induced polarization was notably higher than anticipated in earlier studies using artificial samples. Furthermore, we observed a distinct positive correlation between the coefficient of relaxation time and the increase in sulfide volume, likely due to the geometric characteristics of the sulfide minerals. Our findings suggest that rocks in massive sulfide mounds may generally construct sulfide clusters that lengthen the conductive path of the electrical carrier. Graphical Abstract

Deep-sea massive sulfide deposits formed by hydrothermal fluid circulation are potential metal resources. They can exist not only as mound manifestations on the seafloor (seafloor massive sulfides) but also as embedded anomalies buried beneath the seafloor (embedded massive sulfides). The distribution of embedded massive sulfides is largely unknown, despite their expected high economic value. Recent drilling surveys have revealed a complex model suggesting embedded massive sulfides coexist beneath seafloor massive sulfides. In the coexisting case, geophysical methods are required to distinguish and map both seafloor and embedded massive sulfides for accurate resource estimation. Marine controlled-source electromagnetic (CSEM) methods are useful for mapping massive sulfides as they exhibit higher electrical conductivity compared to the surrounding host rock. However, CSEM applications capable of distinguishing and mapping both massive sulfides are lacking. We employ a towed electric dipole transmitter with two types of receivers: stationary ocean bottom electric (OBE) and short-offset towed receivers. This combination utilizes differences in sensitivity: the towed receiver data are sensitive to seafloor massive sulfides and the stationary OBE receiver data are sensitive to embedded massive sulfides. Our synthetic data example demonstrates that the combined inversion of towed and OBE data can recover resistivities and positions of both massive sulfides more accurately than the existing inversion methods using individual applications. We perform the combined inversion of measured CSEM data obtained from the middle Okinawa Trough. The inversion models demonstrate that a combined inversion can map the location and shape of embedded massive sulfides identified during drilling more accurately than the inversion of individual datasets.

地下比抵抗構造は地下流体の分布や移動を議論する際に有益な情報であるが，流体の粘性の違いが地層の比抵抗に与える影響については，十分には議論されていない。本研究では実験室において，高粘性の泥水を堆積物試料に浸透させ，泥水飽和度と試料比抵抗の関係を調査した。試料として，砂礫混じりの川砂および細粒の砂を用いた。室内実験の結果，泥水使用時の比抵抗は，水道水を用いた時と比べて，低飽和域で高い値を示し，飽和係数も大きくなった。その要因について定性的な議論を行ったところ，泥水の非ニュートン流体（ビンガム流体）的性質により，泥水の毛管現象がほとんど起きなかった事が主な要因であると推測された。高粘性流体が探査対象の場合は，室内実験を通じて，流体の特性や飽和度と比抵抗の関係を調査しておく必要性が示唆された。