A method includes obtaining total stresses and pore pressures of each porous medium of a formation, determining a first and second set of effective stresses for the formation, determining an individual collapse and fracturing mud weight for each porous medium of the formation using a first set of associated failure criteria, wherein the first set of associated failure criteria are based on the first set of effective stresses, determining an overall collapse and fracturing mud weight for the formation using a second set of associated failure criteria, wherein the second set of associated failure criteria is based on the second set of effective stresses, determining a mud weight window for the formation using the individual collapse mud weight, the individual fracturing mud weight, the overall collapse mud weight, and the overall fracturing mud weight, and transmitting a command to a drilling system based on the mud weight window.

Methods and systems for identifying near-surface heterogeneities in a subterranean formation using surface seismic arrays can include: recording raw seismic data using sensors at ground surface; applying a band bass filter to the raw seismic data using a central frequency; picking a phase arrival time for the filtered data; generating an initial starting phase velocity model for tomographic inversion from the raw seismic data; applying tomographic inversion to the filtered data to generate a dispersion map associated at the central frequency; repeating the applying a band bass filter, picking a phase arrival time, generating an initial starting velocity model, and applying tomographic inversion steps for each of a set of central frequencies; and generating a three-dimensional dispersion volume representing near-surface conditions in the subterranean formation by combining the dispersion maps.

Methods and systems for identifying near-surface heterogeneities in a subterranean formation using surface seismic arrays can include: recording raw seismic data using sensors at ground surface; applying a band bass filter to the raw seismic data using a central frequency; picking a phase arrival time for the filtered data; generating an initial starting phase velocity model for tomographic inversion from the raw seismic data; applying tomographic inversion to the filtered data to generate a dispersion map associated at the central frequency; repeating the applying a band bass filter, picking a phase arrival time, generating an initial starting velocity model, and applying tomographic inversion steps for each of a set of central frequencies; and generating a three-dimensional dispersion volume representing near-surface conditions in the subterranean formation by combining the dispersion maps.

A seabed object detection system is provided. The system can include a receiver array including streamers. The system can include a plurality of receivers coupled with the streamers. The system can include a receiver array cross-cable to couple with the first streamer and to couple with the second streamer. The receiver array cross-cable can be disposed at a first depth of a body of water. The system can include a first diverter and a second diverter coupled with the receiver array cross-cable. The system can include a source array including a first source and a second source. The source array can be coplanar to the receiver array. The system can include a source array cross-cable to couple with the first source and to couple with the second source, the source array cross-cable disposed at a second depth of the body of water.

A method or system is configured for determining a seismic attenuation quality factor Q for intervals of subsurface formations by performing actions including receiving vertical seismic profile traces. The actions include filtering the vertical seismic profile traces with an inverse impulse response of a downhole receiver. The actions include transforming the vertical seismic profile data from the particle motion measured by the downhole receiver to the far-field particle motions represented by the source wavelet. The actions include determining a ratio of the spectral amplitudes of the direct arrival event of the transformed vertical seismic profile data and the source Klauder wavelet. A quality factor Q is generated representing an attenuation of the seismic signal between the source at ground level surface and the downhole receiver.

A method or system is configured for determining a seismic attenuation quality factor Q for intervals of subsurface formations by performing actions including receiving vertical seismic profile traces. The actions include filtering the vertical seismic profile traces with an inverse impulse response of a downhole receiver. The actions include transforming the vertical seismic profile data from the particle motion measured by the downhole receiver to the far-field particle motions represented by the source wavelet. The actions include determining a ratio of the spectral amplitudes of the direct arrival event of the transformed vertical seismic profile data and the source Klauder wavelet. A quality factor Q is generated representing an attenuation of the seismic signal between the source at ground level surface and the downhole receiver.

A system can determine a location for future wells using machine-learning techniques. The system can receive seismic data about a subterranean formation and may determine a set of seismic attributes from the seismic data. The system can block the set of seismic attributes into a set of blocked seismic attributes by distributing the set of seismic attributes onto a geo-cellular grid representative of the subterranean formation. The system can apply a trained machine-learning model to the set of blocked seismic attributes to generate a composite seismic parameter. The system can distribute the composite seismic parameter in the subterranean formation to characterize formation locations based on a predicted presence of hydrocarbons.

A media Parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface, comprising the following steps: importing an initial forward model, importing anisotropic parameters; and setting a space step and a time step according to the initial forward model parameters; and then starting a stepped discretization of a free surface of the initial forward model; and using a corrected constitutive relationship to correct a first level parameter of the initial forward model; and bringing the corrected constitutive relationship into a displacement stress equation, and the influence of the free surface can be introduced in the case of the anisotropic media after series of operation. The present disclosure can make an accurate numerical simulation of a wave field near the discontinuous surface, and the accurate numerical simulation will contribute to the extraction and analysis of information from the seismic data.

A media Parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface, comprising the following steps: importing an initial forward model, importing anisotropic parameters; and setting a space step and a time step according to the initial forward model parameters; and then starting a stepped discretization of a free surface of the initial forward model; and using a corrected constitutive relationship to correct a first level parameter of the initial forward model; and bringing the corrected constitutive relationship into a displacement stress equation, and the influence of the free surface can be introduced in the case of the anisotropic media after series of operation. The present disclosure can make an accurate numerical simulation of a wave field near the discontinuous surface, and the accurate numerical simulation will contribute to the extraction and analysis of information from the seismic data.

This disclosure describes a system and method for generating images and location data of a subsurface object using existing infrastructure as a source. Many infrastructure objects (e.g., pipes, cables, conduits, wells, foundation structures) are constructed of rigid materials and have a known shape and location. Additionally these infrastructure objects can have exposed portions that are above or near the surface and readily accessible. A signal generator can be affixed to the exposed portion of the infrastructure object, which induces acoustic energy, or vibrations in the object. The object with affixed signal generator can then be used as a source in performing a subsurface imaging of subsurface objects, which are not exposed.