Geophysical Prospecting最新文献

{"title":"An Improved Equivalent Circuit for Electric Field Sensors in Geophysical Exploration","authors":"Thomas Kalscheuer, Uwe Zimmermann, Henrik Sparr, Anton Palm Ekspong, Alexander Lindblad","doi":"10.1111/1365-2478.70047","DOIUrl":"https://doi.org/10.1111/1365-2478.70047","url":null,"abstract":"<p>In electromagnetic measurements, electric field sensors consist of two halves, with remote electrodes of negative and positive polarity coupled through wires and low pass filters to the differential inputs of an analogue-to-digital converter; the electrical ground of the analogue-to-digital converter is connected to the ground through a reference electrode. We present, analyse and evaluate improved equivalent circuits for such electric field sensors. This serves to identify the maximum contact resistances of the electrodes for which the recorded voltages are unaffected by system response effects over a given frequency range. In the first step, we verify a new equivalent circuit for one half of an electric field sensor by comparison to a previously published equivalent circuit. In contrast to the latter, our equivalent circuit accounts for the spatial variability of the electric field along an extended sensor cable, the finite impedance of the receiver input stage, the non-zero contact resistance of the reference electrode and residual cable on a winch. Furthermore, the cable is characterised by its resistance, self-inductance and capacitance to the ground and the ionosphere or the borehole fluid. Compared to the absolute value of the voltage, our results show that the system response affects the phase of the voltage at lower frequencies. In the next step, we develop an equivalent circuit for a complete electric field sensor connecting two sensor halves to an analogue-to-digital converter. We study both symmetric and asymmetric set-ups with identical and differing cable lengths, respectively, of the sensor halves. Over the whole frequency range, the amplitude gets the lower, the higher the sum of contact resistances of the remote electrodes is. In contrast, the phase is distorted only at higher frequencies. Generally, the contact resistance of the central reference electrode has little effect. For symmetric sensors, of the combinations of contact resistances of the remote electrodes that have the same sum, it is the combination of identical contact resistances that shows the lowest distortion. The distortion owing to different contact resistances of the remote electrodes is only slight and mostly in the amplitude at high frequencies. For asymmetric sensors, the benefits of using a differential analogue-to-digital converter input are no longer exploited. For instance, flipping the contact resistances of the remote electrodes leads to different responses at high frequencies. In borehole applications, it is of particular importance to account for the spatial variability of the electric field due to the skin effect, field propagation and the curvature of the borehole track. We consider an extended electric field sensor that is placed in a borehole at an inclination of <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>45</mn>\u0000 <msup>\u0000 <mspace></mspace>\u0000 ","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 8","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70047","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145224080","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Low-Frequency Extrapolation by Deep-Learning for Cross-Well Full-Waveform Inversion – Case Study From the Aquistore \u0000 \u0000 \u0000 CO\u0000 2\u0000 \u0000 $rm CO_{2}$\u0000 Storage Site","authors":"Amir Mardan,&nbsp;Don White","doi":"10.1111/1365-2478.70081","DOIUrl":"https://doi.org/10.1111/1365-2478.70081","url":null,"abstract":"<p>Full-waveform inversion (FWI) of seismic data is a powerful method for estimating high-resolution models of the subsurface. An accurate initial model and low-frequency data are necessary to avoid cycle skipping and perform a successful FWI. In the absence of this information, FWI is likely to fail due to convergence in local misfit minima. With the recent advancements in artificial intelligence, studies have shown that absent low-frequency data can be extrapolated using deep learning (DL). These studies have been mostly focused on surface seismic data whose frequency content is different from cross-well data. In this study, we assess the use of DL for low-frequency extrapolation for a cross-well survey that was done at the Aquistore <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>CO</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <annotation>$rm CO_{2}$</annotation>\u0000 </semantics></math> storage site in Saskatchewan. This assessment includes both numerical and field data examples. We extrapolate the low frequencies to increase the bandwidth of the acquired data at the Aquistore site and perform FWI. We evaluate the efficiency of this method by comparing the results with obtained velocity models from the conventional multiscale FWI. Our results for the Aquistore data show that the proposed strategy leads to an accuracy improvement of 39% and 20% in the model and data domains, respectively.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70081","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145146338","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Improved One-Way Reflection Waveform Inversion and Strategies for Optimal Offset Selection","authors":"Siamak Abolhassani,&nbsp;Dirk Jacob Verschuur","doi":"10.1111/1365-2478.70062","DOIUrl":"https://doi.org/10.1111/1365-2478.70062","url":null,"abstract":"<p>Conventional reflection waveform inversion solves a two-parameter seismic inverse problem alternately for subsurface reflectivity and acoustic background velocity as the model parameters. It seeks to reconstruct a low-wavenumber velocity model of the subsurface from pure reflection data cyclically, through alternating migration and tomography loops, such that the remodelled data fits the observed data. Low-resolution seismic images with unpreserved amplitudes, full-wave inconsistency in the short-offset data and cycle skipping in the long-offset are perceived as the main reasons for suboptimal tomographic updates and slow convergence in conventional reflection waveform inversion. In the context of one-way reflection waveform inversion, this paper addresses the listed limitations through four main components. First, it augments one-way reflection waveform inversion with a computationally affordable preconditioned least-squares wave equation migration algorithm to ensure high-resolution reflectors with preserved amplitudes. Second, the paper verifies how well the full-wave consistency condition in the short-offset data is satisfied in one-way reflection waveform inversion and suggests muting inconsistent short-offset residual waveforms in the tomography loop to attenuate their adverse imprint. Third, the paper suggests extending the migration offset beyond short offsets to improve both the illumination and the signal-to-noise ratio of the reflectors. Fourth, the paper presents a data-selection algorithm to exclude the damaging effect of the cycle-skipped long-offset data in the tomography loop. The effectiveness of the proposed one-way reflection waveform inversion algorithm is finally validated through three numerical examples, demonstrating its capability to recover high-fidelity tomograms.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70062","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145101307","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Forward Modelling of Electrical Resistivity and Induced Polarization Using the Spectral-Infinite-Element Method","authors":"Kiana Damavandi,&nbsp;Hom Nath Gharti","doi":"10.1111/1365-2478.70080","DOIUrl":"https://doi.org/10.1111/1365-2478.70080","url":null,"abstract":"<p>Accurate and efficient modelling of subsurface electrical properties is critical for a wide range of applications, including mineral exploration, environmental studies and hydrogeological investigations. Traditional numerical approaches often use low-order discretization and impose artificial boundary conditions to approximate the unbounded spatial domain. These approximations can lead to inaccuracies and computational inefficiency, particularly in geologically complex environments. In this study, we present a spectral-infinite-element method (SIEM) for forward modelling of electrical resistivity and induced polarization. The approach couples high-order spectral elements within the finite domain with a single outer layer of mapped infinite elements, enabling precise representation of far-field boundary conditions. To achieve optimal numerical performance, we employ two distinct quadrature schemes: Gauss–Legendre–Lobatto quadrature for the spectral elements and Gauss–Radau quadrature for the infinite elements. We first verify the accuracy of our method by comparing the computed electric potential from a buried charged block with direct numerical integration. We conducted a convergence study by refining the mesh and increasing the order of the interpolation polynomials. To further evaluate the robustness of SIEM, we benchmark its results for a layered earth model against an analytical solution and an open-source Python-based geophysical modelling library, SimPEG. The comparisons demonstrate the accuracy, convergence and efficiency of SIEM. Finally, we apply SIEM to a complex heterogeneous conductivity model incorporating topography, generating apparent resistivity and chargeability pseudo-sections to illustrate its practical applicability under realistic survey conditions.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70080","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145101213","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Constraining Contact-Depth Solutions Through the Enhanced Horizontal Gradient Amplitude","authors":"Luan Thanh Pham,&nbsp;Saulo Pomponet Oliveira,&nbsp;Erdinc Oksum,&nbsp;Hanbing Ai,&nbsp;Fabiana de Fátima Giacomini","doi":"10.1111/1365-2478.70078","DOIUrl":"https://doi.org/10.1111/1365-2478.70078","url":null,"abstract":"<div>\u0000 \u0000 <p>Accurate depth estimation is crucial for the quantitative interpretation of magnetic anomalies, which plays a significant role in geological mapping, mineral exploration and subsurface investigations. Traditional depth estimation techniques, such as the contact-depth (CD) and tilt-depth (TD) methods, often suffer from the generation of spurious solutions, especially when applied to complex geological environments. To address this, we propose an enhanced depth estimation technique, namely, the located contact-depth (LCD) method that integrates the CD technique with the enhanced horizontal gradient amplitude (EHGA). By utilizing points near the peaks of EHGA, a mask is generated to constrain the solutions from the CD method, effectively eliminating false solutions. Furthermore, a stable finite-difference technique for calculating vertical derivatives is used to improve the robustness and stability of the outputs. The proposed technique is tested on synthetic data, both with and without noise, as well as on real aeromagnetic data from the Galinge Fe-polymetallic deposit (China). The results demonstrate that our method provides depth estimates with improved reliability and accuracy compared to traditional methods, reducing the number of spurious solutions and enhancing precision around source boundaries. The result from the real example is in good agreement with known structures, highlighting the potential for deep mineral exploration in the Galinge Fe-polymetallic deposit.</p>\u0000 </div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145022016","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Simple Method for Estimating Shale Brittleness From Seismic Data: An Example From Offshore Norway","authors":"Md. Jamilur Rahman,&nbsp;Owen Rohwer Huff,&nbsp;Sajjad Ahmadigoltapeh,&nbsp;Jan Inge Faleide,&nbsp;Manzar Fawad,&nbsp;Nazmul Haque Mondol","doi":"10.1111/1365-2478.70076","DOIUrl":"https://doi.org/10.1111/1365-2478.70076","url":null,"abstract":"<div>\u0000 \u0000 <p>The brittleness of a caprock layer is an important property to measure, as it influences sealing integrity for subsurface fluid injection projects. However, brittleness is a complex function of various factors, such as mineral composition, diagenesis and effective stress, and consequently can vary spatially. Because well log-based methods for brittleness estimation are often spatially limited, we develop a simple workflow to estimate the brittleness of a shale formation directly from seismic data. The shale in question is the Draupne Formation of the Upper Jurassic age, which acts as a primary seal of Viking Group sandstones in the Horda Platform area of the northern North Sea for hydrocarbon extraction (e.g., the Troll field) and geological CO₂ storage (e.g., Smeaheia). First, well log data from 26 wells in the Horda Platform area is aggregated, focusing on the compressional sonic, bulk density and resistivity values of the Draupne Formation. This data are used to establish a linear model relating the acoustic impedance and elastic properties-estimated brittleness index (BI) of the Draupne Formation; these two quantities display a correlation of 0.86. By combining this acoustic impedance information with a wavelet extracted from field seismic data and using average acoustic properties for the relatively homogeneous underlying sandstone reservoir, synthetic seismograms corresponding to different BI values of the Draupne Formation are generated. The amplitudes extracted from the synthetic seismograms are then used to establish a quadratic model relating seismic amplitudes at the base Draupne reflection with the BI. Applying this quadratic model on a 2D seismic line from the Stord Basin (south of the Horda Platform) results in BI values that are close to elastic-properties-based values at wells which intersect the seismic line and an expected trend of increasing brittleness with respect to depth. This integrated method can be used as part of a workflow to characterize top seal effectiveness, which may be useful in fluid storage prospect evaluation.</p>\u0000 </div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145012255","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Image-Domain Least-Squares Migration Through Preconditioned Hessian","authors":"Wei Zhang,&nbsp;Mauricio D. Sacchi","doi":"10.1111/1365-2478.70071","DOIUrl":"https://doi.org/10.1111/1365-2478.70071","url":null,"abstract":"<div>\u0000 \u0000 <p>Image-domain least-squares migration (IDLSM), which typically employs a diagonally dominant Hessian with narrow bandwidth for the inverse problem, provides an efficient deconvolution strategy for subsurface reflectivity imaging. Conventional methods often rely on the adjoint of the Born/Kirchhoff modelling operator to compute the Hessian matrix. However, the adjoint-derived Hessian is highly ill-conditioned, leading to slow convergence during linear inversion and resulting in images with undesired resolution and amplitude fidelity. To overcome these limitations, this study introduces a novel IDLSM approach that integrates the state-of-the-art migration operator. We derive and compute the preconditioned Hessian matrix through a Kirchhoff migration engine with source-side and receiver-side illumination. The preconditioned Hessian matrix exhibits identical values along its main diagonal. This illumination compensation will explicitly reduce the condition number of the Hessian matrix and significantly improve the quality of migrated images in terms of amplitude fidelity. In addition, we remove redundant source wavelets from the migrated image and the Hessian matrix. As a result, these improvements will greatly accelerate the convergence of linear inversion solvers while enhancing the resolution and amplitude fidelity of the resulting images. Experiments on synthetic and field datasets demonstrate that the proposed IDLSM method retrieves high-fidelity reflectivity images with superior resolution and amplitude fidelity compared to conventional IDLSM techniques.</p></div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145012344","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Case Study of Dual-Signal Processing of DAS-VSP Vibrator Data From a 3D Survey in a Geothermal Reservoir","authors":"Flavio Poletto,&nbsp;Cinzia Bellezza,&nbsp;Gualtiero Böhm,&nbsp;Fabio Meneghini,&nbsp;Athena Chalari,&nbsp;Anna Stork,&nbsp;Mahmut Parlaktuna,&nbsp;Erdinç Şentürk,&nbsp;Deyan Draganov,&nbsp;Gijs van Otten,&nbsp;Sevket Durucan","doi":"10.1111/1365-2478.70073","DOIUrl":"https://doi.org/10.1111/1365-2478.70073","url":null,"abstract":"<p>The three-dimensional (3D) distributed acoustic sensing (DAS) vertical seismic profile (VSP) technique is an effective tool to characterize subsurface reservoirs, enabling the use of large and densely sampled borehole receiver arrays with many surface vibrator source points for onshore time-lapse monitoring. However, the processing of the DAS VSP signals for imaging purposes is based on a reliable wavefield separation, which may depend on the recognition and quality of the direct arrivals. To overcome this limitation for common-source gathers with poor signal-to- noise ratio or with interferences, we apply the dual-signal processing method, which allows us to estimate and separate the DAS wavefields by signals' combination without arrival picking. We present a case study of a 3D VSP DAS dataset recorded at a geothermal reservoir in Turkey, showing that the method, similar to a geophone and hydrophone combination, is robust and effective and can be advantageously integrated with the conventional processing. Supported by signal benchmarking, modelling and signal-to-noise ratio analysis, we treat common-source and common-receiver data. Our analysis shows the advantages and limitations of the proposed approach, valuable in the time-lapse perspective.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70073","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144998923","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A 3D Joint Tomographic Inversion of First-Arrival and Reflection Waves Based On the Adjoint State Method","authors":"Junjie Sun,&nbsp;Huachen Yang,&nbsp;Fei Ma,&nbsp;Jianzhong Zhang","doi":"10.1111/1365-2478.70077","DOIUrl":"https://doi.org/10.1111/1365-2478.70077","url":null,"abstract":"<div>\u0000 \u0000 <p>Seismic reflection traveltime tomography (RTT) is an effective technique for inverting subsurface low-frequency velocity models for prestack depth migration and full-waveform inversion of seismic data. However, the velocity model established using RTT demonstrates limited resolution for extremely shallow and deeply complex strata. Small-offset first-arrival wave effectively characterize velocity variations in shallow strata, whereas large-offset first-arrival wave can reflect the velocity distribution in deeper strata. Therefore, we propose a three-dimensional joint tomographic inversion of first-arrival and reflection waves based on the adjoint state method in this article. The method integrates first-arrival traveltime data, reflection traveltime data and slope data for inversion, enhancing the accuracy of the inversion model from shallow to deep. The adjoint state method is employed to calculate the gradient of the misfit function with respect to velocity and the spatial coordinates of reflection points, thereby reducing the computational memory requirements and improving the efficiency of velocity modelling. The results of synthetic data tests based on a theoretical model verify the accuracy and effectiveness of the proposed joint inversion method. The proposed method is applied to ocean bottom node and ocean bottom cable wide-line seismic data collected in an extremely shallow sea in eastern China, yielding a more accurate velocity model.</p>\u0000 </div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144935129","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Dynamic Streaming Potential Coupling Coefficient in Partially Saturated Porous Media","authors":"Luong Duy Thanh,&nbsp;Santiago G. Solazzi,&nbsp;Nguyen Manh Hung,&nbsp;Nguyen Van Nghia,&nbsp;Phan Van Do,&nbsp;Damien Jougnot","doi":"10.1111/1365-2478.70072","DOIUrl":"https://doi.org/10.1111/1365-2478.70072","url":null,"abstract":"<div>\u0000 \u0000 <p>The seismoelectric effect is an electrokinetic phenomenon that arises when seismic waves propagate in water-containing geological formations. Given that seismoelectric signals are sensitive to the hydraulic properties of the probed porous medium, they have the capability to provide important information during subsurface characterization efforts. In this work, we present a physics-based model for the dynamic streaming potential coupling coefficient (SPCC) in partially saturated porous media. For this, we conceptualize the porous medium as a partially saturated bundle of capillary tubes. We take into account the variation of pore size to relate the capillary pressure to the water saturation in the porous medium of interest. We then up-scale the streaming current and conduction current within the saturated capillaries under oscillatory flow conditions from pore to sample scale. The results show that the dynamic SPCC is not only a function of water saturation and the probing frequencies but also of the properties of water, mineral–water interfaces and other microstructural parameters of the porous medium. We analyse and explain the characteristics of the dynamic SPCC for two different pore size distributions (PSD): fractal and lognormal. Results show that the PSD characteristics have a strong effect on the dynamic SPCC responses. The proposed model has a remarkable ability to replicate experimental data available in the literature. In addition, it is observed that the lognormal distribution can provide a better agreement with experimental data for sand samples, which display a relatively narrow PSD. The findings of this study provide a valuable basis for interpreting seismoelectric signals under partially saturated conditions. Our proposed technique can be applied to any PSD, regardless of the complexity, providing a flexibility that is not present in alternative models found in the literature.</p></div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 7","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144923624","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}