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.

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 and system for forming a seismic image of a subterranean region are disclosed. The method includes determining an initial plan for a seismic survey with a value for each member of a set of acquisition parameters and acquiring a first seismic dataset from a first portion of the seismic survey based on the initial plan. The method further includes transmitting the first seismic dataset to a seismic processor, determining a first seismic image from the first seismic dataset by performing expedited seismic processing and determining a first updated plan for the seismic survey based on the first seismic image and acquiring a second seismic dataset from a second portion of the seismic survey based on the first updated plan. The method still further includes transmitting the second seismic dataset to the seismic processor and determining the seismic based on the first seismic dataset and the second seismic dataset.

A method includes receiving well log data for a plurality of wells. A flag is generated based at least partially on the well log data. The wells are sorted into groups based at least partially on the well log data, the flag, or both. A model is built for each of the wells based at least partially on the well log data, the flag, and the groups.

Methods and systems, including computer programs encoded on a computer storage medium can be used for an integrated methodology that can be used by a computing system to automate processes for generating, and updating (e.g., in real-time), subsurface reservoir models. The methodology and automated approaches employ technologies relating to machine learning and artificial intelligence (AI) to process seismic data and information relating to seismic facies.

Systems and methods to estimate a pore pressure of source rock include a pore pressure estimation processor, an executable, or both, and are operable to (i) calculate an estimate pore pressure based on overburden gradient data, a compaction velocity profile, hydrocarbon maturity, and an unloading velocity profile, (ii) determine a total organic content (TOC) estimate of the source rock based on a bulk density at a vertical depth measured using the density logging tool, (iii) determine a correction factor based on (a) the TOC estimate and (b) vitrinite ratio R.sub.o data, and (iv) update the estimated pore pressure in real-time based on the correction factor.

A computer implemented method for analyzing a reservoir grid modeling a reservoir geological formation is provided in which the reservoir grid corresponds to a 3D grid of cells associated to respective values of at least one geological property. The method includes obtaining a 4D seismic image of the reservoir geological formation. A skeleton of the 4D seismic image is calculated, and the skeleton extends between at least one origin and a plurality of extremities. Each point of the skeleton is associated to a value of the at least one geological property of the reservoir grid. Flow time values are calculated for a fluid flowing from the origin to the extremities along the skeleton, based on the at least one geological property values associated to the points of the skeleton. The reservoir grid is calculated based on the flow time values.

Implementations for monitoring seismic data recorded in a marine survey of a subterranean formation for coverage gaps are described herein. Implementations include computing Fresnel sum operators for Fresnel zones of the subterranean formation based on a Kirchhoff migration impulse response at horizons of a representative plane layer model of a survey area of the subterranean formation. Implementations also include computing an acceptability map of the survey area based on the Fresnel sum operators. The acceptability map reveals coverage gaps in the survey area. Geoscientist may use the acceptability map to infill seismic data in areas of the survey area that correspond to the coverage gaps.

A method includes acquiring seismic data of a region that utilizes multiple seismic energy sources and seismic energy receivers where the seismic data include blended seismic data for a number of emissions from a corresponding number of the multiple seismic energy sources; determining spatially distributed coherent noise properties for the region using the blended seismic data; via the spatially distributed coherent noise properties, modeling coherent noise as at least two coherent noise models for at least two of the emissions from a corresponding at least two of the multiple seismic energy sources; via the coherent noise models, attenuating coherent noise in a portion of the blended seismic data to generate coherent noise attenuated blended seismic data; deblending the coherent noise attenuated blended seismic data to generate deblended seismic data; and rendering an image of at least a portion of the region to a display using the deblended seismic data.

A method for establishing a geostress field distribution of slopes in canyon areas includes: obtaining a persistence ratio of a fracture surface based on a structural plane trace length and a rock bridge length of the fracture surface, and then obtaining a fracture stage of a crack according to progressive failure characteristics of rock mass, combining a character of the fracture surface to obtain magnitude and direction of a maximum principal stress, and establishing the geostress field distribution. The method is simple to operate, does not need to carry out geostress testing, does not need a large amount of manpower and material resources, does not need redundant fund investment, and can simply and effectively obtain geostress field data. Moreover, combining with the geostress field inversion technology, a large-scale geostress field distribution condition can be obtained, which can provide a basis for engineering site selection and engineering rock mass stability determination.