Pore pressure prediction based on rock stress applied to Marmousi seismic data

This paper focuses on the modeling of a sedimentary basin where the exploration of oil and gas takes place. In the modeling, we calculate the vertical variation of pore and rock pressure, which serve as natural pumps for the accumulation of fluids, where we use post-migration surface seismic data and information as necessary. We compare two methods of pore pressure calculation. In the first case, we calculate stress due to vertical loading created by gravitational geological overburden. In the second case, we propose a more complex method to calculate the stress distribution based on the mechanics of solids under gravity loading, using the concept of the first tensor invariant and the linear behavior of Hooke’s law. We prove the use of the P and S velocities and density information to calculate rock, pore, and effective pressure distribution, useful for characterizing potential zones for oil and gas accumulation. The proposed method allows formulation of rock pressure from different calculations instead of as a simple overburden pressure value. The joint analysis of the calculated sections can be used to identify patterns and correlations, outline and characterize the target zones, and make practical geological conclusions. sandstone, and limestone. Several efforts are made for predicting pressure using data,


INTRODUCTION
The main purpose of this work is to calculate pore pressure in a sedimentary basin, where the exploration of oil and gas takes places. Besides, we also aim at calculating low-and high-pressure zones that serve as natural pumps that accumulate fluids in the subsurface. To meet this goal, we calculate changes in stress (and, therefore, pressure types) in the subsurface based on surface seismic data.
We compare two methods of pore pressure calculation. In the first case, we calculate stress due to the vertical loading created by the gravitational geological overburden. In the second case, we propose a more complex method to calculate the stress distribution based on the mechanics of solids under gravity loading, using the concept of the first tensor invariant and the linear behavior of Hooke's law.
We admit that the mechanics of solids is the proper method to describe the stress state in sedimentary basins and predict low-and high-pressure zones, where a reservoir is positioned in a low-pressure zone and capped by a high-pressure zone. Also, it is proper to consider a normal compaction trend (NCT) in the modeling and that the NCT curve is specific to a given lithology.
This work is more than a numerical experiment dedicated to pore pressure prediction, where the aim is to map pore pressure variations in the subsurface; this work is also to intended to show low-and high-pressure zones and correlate them to the mapped pore pressure profile. The method was applied to the Marmousi data described by Verrier and Branco (1972), Martin (2004), and Martin et al. (2006), among other authors. Their data supplies the information to the construction of stress prediction profiles using V p and V s velocities and density ρ distributions as the necessary data, as described by . Therefore, it is natural to classify the method as post-migration processing.
The method starts with the lithological identification for the construction of a porosity model, followed by least-squares regression to obtain the NCT. Then, we calculate the hydrostatic pressure, the alternative overburden weight, and pore and effective pressures.
The applied specific method is described from exploration company problems and the research literature classify pore pressure prediction as a relevant subject in oil and gas exploration. Besides, it consists of porosity and gravity loading phenomena, as described by Zhang (2013). Carcione and Helle (2002) and Carcione et al. (2003) address the subject of pore pressure based on a physical model that correlates porosity and elastic parameters that behave in the form of a depth-dependent exponential function, where the systematic use of the Hashin-Strickman model is made to describe upper and lower bounds for the rock parameters, as summarized by Mavko et al. (1999).
A method for calculating pore pressure is based on the NCT exponential fitting as described by Athy (1930) and Zhang (2011), where they call attention to the lithologies: clay, sandstone, and limestone. Several efforts are made for predicting pore pressure using porosity data, and Heppard et al. (1998) used an empirical porosity equation that is like Eaton's method using shale porosity results to predict pore pressure. We can also add the importance of this subject to engineering aspects, as given by Holbrook et al. (2005), where they describe the use of petrophysical data and stress-strain relationships for pore pressure prediction in real-time drilling. Schneider et al. (2009) also describe porosity-stress relationships to predict over-pressure in claystone.
The prediction of low-and high-pressure zones and pore pressure variation in sedimentary basins for gas and oil exploration has many theoretical aspects related to engineering, geology, and geophysics, aiming to characterize potential reservoirs in the subsurface. Furthermore, the prediction and monitoring of abnormal pore pressure is a concern to avoid serious drilling hazards. Abnormally high pore pressures can cause blowouts, besides inducing geologic disasters, such as mud volcano eruptions. Sibiryakov et al. (2020) describe applications of the theory of solid mechanics for pressure prediction in sedimentary basins, where the purpose is to show a mechanism for driving fluids from high-to low-pressure zones. Besides, it is necessary to analyze the participation of the so-called "effective pressure" in this fluid mechanism. As a first step in this direction, we start by defining the term "pressure", which is a phenomenon present in fluids (oil and gas) and solid porous rock (matrix, skeleton), being a concept directly related to stress and gravity loading. Thus, pressure stands for a scalar quantity, and stress for a vectorial field.

Solid mechanics and seismic modeling
The size of the hydrostatic pressure,P h , is controlled by the fluid column (h = z -z 0 ) from the surface (z 0 ) to the depth of interest (z), the fluid density [ρ f (z) ≈ ρ f ] and the gravity acceleration [g(z) ≈ g], expressed by the overload formula given by Equation 1.
The dimensions (and units) of pressure and stress are always [FL 2 ], but the physical concepts can be different.
The overburden weight plays the role of vertical normal stress (P z = σ zz ) and it is defined by the total weight of the rock formation and fluids down to the reference point, usually taken around the reservoir. This model does not integrate the lateral variation in physical parameters, like density. Therefore, considering that the bulk density varies only with the depth, the overburden stress can be calculated by integrating the density function to the depth z, i. e. (Eq. 2): May 25, 2022 quação 2 da maneira correta.
equações 6 e 7 da maneira correta. (2) where we included a liquid layer (over the rock basin formations) with water density ρ w , height of water column h w , seabed depth z w (for onshore drilling z w = 0 and h w = 0), and depth of interest z. The gravity effect g = g(z) is considered constant (very small variation) in the interval (z w , z).
It is relevant to note that for a horizontal plane stratified earth, the horizontal components of stress are null, and the loading action is due only to σ zz . This means that the structure is constrained laterally and that the normal are symmetrical components. Therefore, as an approximation, only gentle dips should be admitted in the model.
The formation rock bulk density function, ρ b (z), is admitted to vary only with depth (as a first approximation) and given by Equation 3: Where: ϕ = the porosity of the rock; ρ m = the matrix (or mineral) density; ρ f = the fluid (water, oil, or gas) density.
If the fluid is water, ρ f is the density of formation water, which usually varies from 1.0 to 1.08 g/cm 3 in sedimentary basins. The bulk density, ρ b , can also be obtained from the density log, but these measures are not always available, or only for a borehole representing a limited area of the basin; therefore, we can make use of empirical prediction models.
There is confusion about the use of the word "pressure" in the geological and geophysical literature; therefore, it is convenient to add a second word to complement the concept behind the word "pressure" to make the physical description clearer. This terminology issue may derive from the fact that stress and pressure have the same physical dimension (force per unit area). The term "pressure" is clear for the hydrostatic (a scalar quantity) case. However, for rocks, we have nine elastic stresses (a second-order tensor) in the mechanics of solids.
Some authors use the term rock pressure, P r , (instead of pressure), and we did the same in this paper, defined as the simple average of the normal stresses, and this corresponds to one-third of the first stress tensor invariant (I 1 = σ xx + σ yy + For the hydrostatic case, σ xx = σ yy = σ zz , then, Pr has a clear physical meaning as a scalar. For the 3D geological medium, admitting the case of a flat-stratified earth model, the stress is non-hydrostatic and submitted to the symmetry condition σ xx = σ yy . In general, we can consider that σ xx , σ yy , and σ zz can all have different values. We call the relevant vertical normal component stress σ zz = P z (Persen 1965).

Horizontal stress
Another problem with the concept of gravitational pressure on rock formations is that rock pressure (4), P r , is simply taken as the loading σ zz ; that is, P r ≡ σ zz , without taking into account horizontal stress components σ xx and σ yy . Since we include this effect in our definition of rock pressure, we can show that the horizontal normal stress components (σ xx and σ yy ) are sufficiently lower than σ zz and given by Equation 6:
where we consider the lateral symmetry σ xx = σ yy , and the velocity ratio = (Sibiryakov et al. 2020).
Thus, using an algebra method with Equations 4 and 6 and using Hooke's law (for the 3D case), we obtain an expression to represent the rock pressure as given by quações 6 e 7 da maneira correta.
where the seismic information enters through the γ ratio. Equation 7 shows that P rv is significantly reduced concerning the simple overlying load P z = σ zz , since γ < 1. For example, for the Poisson solid, μ = λ, γ = 0.57735; ≈ 1 2 , i. e., P rv corresponds approximately to half the value of P r , which is a big difference. This means the rock pressure, P rv , depends strongly on the velocity γ ratio . It should be clear that pressure on solids (rocks) does not include the effect of fluid (liquid and gas) in the different definitions.
From the rationale around Equations 4 and 7, the concept of rock pressure (or simply pressure in solids under gravity effect) does not have a clear physical meaning; nevertheless, it can be used to represent a loading force for driving fluids. A similar line of thought was honored by Cibin et al. (2011). A relevant observation is the limitations of the models since the layers are considered horizontal (or with a very smooth dip); also, fracturing and chemical transformations are not admitted in the modeling.
As a partial conclusion, loads P z (x, z), P r (x, z), and P rv (x, z) should be compared as models for the driving fluid mechanism in the effective pore pressure concept, where P z (x, z) is dependent on the density distribution, and P rv (x, z) needs the addition of information about V p (x, z) and V s (x, z) seismic velocities. But, these loads should be calculated without considering the presence and effect of fluids.
Effective and pore pressures Xu et al. (2006) and Li et al. (2009) showed in laboratory how transversal, V s , and longitudinal, V p , wave velocities in micro-inhomogeneous medium containing fluids depend on the difference between the rock matrix pressure [the three possibilities are P z , P r , or P rv , all denominated by P sca (and "sca" from scalar)] and the fluid pressure (or pore pressure, P p ).
The effective (or differential) pressure, P eff , is the net pressure that acts on the solid rock matrix (skeleton) formed by the grains. According to Terzaghi (1943)'s principle, P eff is defined as the difference between a general scalar overburden pressure, P sca , and the pore pressure P p , namely (Eq. 8): We must assert that P p is a natural scalar and that P sca is also a scalar under the definition of the mean of the first stress invariant. The experimental results from Xu et al. (2006) and Li et al. (2009) show that sometimes it is necessary to introduce a dimensionless variable factor, κ, of unknown physical nature, based only on empirical evidence to reduce the participation of P p for better data fitting, so that Equation 8 is rewritten as Equation 9: It should also be clear that the model does not consider fluid flow.
Our proposed alternative is to change Equations 9 and 8 and use the rock pressures P rv (and P r ) in Equations 8 and 9. We can obviously rewrite them as P eff = P rv -P p and P effκ = P rv -κP p (we do not estimate κ in this work). The meaning of effective pressure is related to the part of the stress tensor which produces measurable effects such as compaction of sedimentary rocks, or an increase in shearing resistance. If P eff is reduced, the compaction rate and pressure generation are also reduced; we can argue that P eff can be considered a driving fluid mechanism.
As a basic principle, there are bounds for the variation of pore pressure in a sedimentary basin. As a reference, pore pressure can vary from the hydrostatic pressure to more than two times the hydrostatic pressure in a rock formation. When the pore pressure is equal or approximately equal to the hydrostatic pressure, the sediments are considered normally pressured, and the pore pressure is called normal pore pressure.
Following this principle, a usual qualitative description of the geological model for the sedimentation phenomenon and pressure variation with depth is as follows, and depicted in Figure 1. Then, the question is how to transform this qualitative model into a quantitative description. The initial marine sediments are formed by unconsolidated material and have high permeability and porosity. This means that the water is not confined in the pore spaces, and has a pressure communication with the seafloor surface. While compaction occurs, the sediments are buried at greater depths, and the increase of overburden load acting in the grain contacts results in lower porosity, lower permeability, and fluid ejection. If the sedimentation rate is higher than the rate at which fluids are expelled from the pores, or seals avoid dewatering of rock formations, the fluids in the pores become over pressured and bear part of the load of the overburden sediments. This compaction disequilibrium generates abnormally high pore pressure zones, with higher porosity than the normally compacted formations. Figure 1 illustrates this simple description; but, here we look at details that can contribute to more realistic pressure and stress distributions in a sedimentary basin based on seismic information (Carcione and Gangi 2000). The adopted model from solid mechanics for the description of pressure variation with depth is related to the parameter variations that directly define jumps in the = ratio profile (V s and V p are the velocities of S and P waves, respectively). Therefore, this method will take us to a more consistent physical model to explain geological phenomena and, consequently, the engineering potential applications as described by Sibiryakov B.P. et al. (2015), Sibiryakov E.P. et al. (2015), and Vieira et al. (2017).
As further details, Sibiryakov (2014) and Sibiryakov and Sibiryakov (2018) addressed the problem of pore pressure and conditions for the disjoining pressure as the solution of a boundary value problem of the mixed type, with new computer technology for solving three-dimensional problems for the elastodynamic equation under stationary oscillations. The use of this new approach allowed the analysis of a few problems with the driving pressure in the granular medium, and they showed that the variation of properties of the static stress state under the influence of pore pressure depends essentially on the contact geometry and only a little on the type of boundary condition at the grain contact.

Estimation of porosity from density
The classical bulk density (or rock density) is defined as the total mass contained in a given unit volume of the rock. In a porous rock, this measure depends on the grain and fluid densities forming the rock, as expressed by the Equation 3. The bulk density is a relevant parameter for obtaining the overburden (or lithostatic) pressure.
In the present case, we are interested in estimating porosity from the density profile since it is a more direct measurement. Therefore, isolating for the porosity in the Equation 3, we obtain the Equation 10: The porosity estimate from the Equation 10 depends on the knowledge of the fluid density in the formations, ρ f , and the density of the matrix (or main mineral), ρ m , for a given rock type. Table 1 is an example of a list of matrix materials formed by common minerals. For the most solid granitic type of rocks, the matrix average density is around 2.65 g/cm 3 , also approximately for quartz. Clay minerals predominate in shale composition.
NCT and pore pressure from porosity and loading To predict abnormal pore pressure generated by the compaction disequilibrium model, we need to obtain the porosity under normal compaction conditions. Athy (1930) proposed a relation for the decrease of porosity with depth under normal pressure conditions (P z = σ zz ) with an exponential behavior given by Equation 11.
Where: ϕ n = the porosity in normal compacted formation; ϕ 0 = the porosity at the mudline (seafloor); c = a constant that represents the compaction factor; z = the true vertical depth below the sedimentary depositional mudline (x fixed).
The law expressed by Equation 11 is widely applied for shales because their normal porosity profiles generally show a concave downward curvature, as shown in Figure 1. Therefore, the main goal when analyzing compaction trends is to identify

4/12
and evaluate the parameters that control the mathematical expression for normally compacted sediments, as in the case of shales. In the data used in this work, we have a variety of limestones, sandstones, and shales.
There are several methods for pore pressure prediction from formation porosity, and here we followed the method proposed by Zhang (2011) with some modification, which derived from the NCT of porosity proposed by Athy (1930) using Equation 11. Following this method, the pore pressure, the overburden pressure, and the porosity relate to one another by the optional expressions given by Equations 12 and 13: Where: P p (rv) and P p (z) = the pore pressure models; P rv = our present proposed overburden pressure model; P h = the hydrostatic pressure (or the normal pore pressure); ϕ 0 = the initial porosity (at seabed level); ϕ = the porosity at the depth of interest z and c is compaction factor.
The main point in the interpretation for the Equations 12 and 13 is that the calculated pore pressure P p (z) has porosity as depth-dependent. If the porosity ϕ(z) (at a target depth z) is greater than the normal porosity ϕ n for the same depth, this zone corresponds to an over-pressure; otherwise, to an under-pressure.
From the concept of effective pressure, in this work, we compare the Equations 12 and 13 calculated with P sca = P rv and with P sca = P z and present as results. The vast literature uses only P z for the overburden loading in the geological description for the mechanics of sedimentation. However, from the concepts of mechanics of solids, this description is incomplete due to the horizontally developed stresses.

Flowchart for pressure distribution calculation
The Marmousi2 data (an update of Marmousi data) is a collection of seismic distribution of the P and S wave velocities [V p (x, z) and V s (x, z), respectively], and the distribution of density [ρ(x, z)] built from a region of the Cuanza basin (Angola), as described in the sequel. The basin is characterized by intense tectonics, a clear sequence of lithological units with well-defined layer geometries. In such a way, with the distribution of physical parameters and geometry, forward modeling allows different studies in a complex and important geological environment.
The following flowchart and block diagram in Figure 2 describe the main steps for the modeling of pore and effective pressures applied to the Marmousi2 seismic data, and it represents a summary of the methodology. We use a 3D grid (Δx, Δy, Δz) for all calculations, but the transversal Y direction is formed by only one cell (Δy), while along with the (x, z) coordinates we have (M x N) grid cells (brick wall type of profile).

Flowchart
• Field and seismic parameters of the migrated seismic section; • The input of V p (x, z), V s (x, z), ρ(x, z) distributions in an equal mathematical matrix format; • Separation of the main lithologies in the seismic section; • Calculation of the porosity distribution; • Least-square fitting of the NCT based on the chosen geological formation: definition of the parameters ϕ 0 (porosity at the mudline level) and c (compaction factor for the geological formation); • Calculation of different pressure models: hydrostatic P h (x, z), overburden P z (x, z), and rock P rv (x, z) as a function of the velocity γ(x, z) ratio; • Calculation in loop of the pore pressure P p

RESULTS
We present our findings in a sequence of color figures that also show the modeling steps taken in mapping the lowand high-pressure zones. In general, the target is to point to new potential oil and gas accumulations or extend production reservoirs.

Geological and seismic information
The Marmousi model was created at the Institut Français du Pétrole (IFP) in 1988, based on a profile through the North Figure 2. Block diagram showing the main steps for pore pressure prediction and the relationships between the input data, calculated parameters, and pressure sections.

5/12
Quenguela Trough in the Cuanza basin (Angola), as described in the work by Verrier and Branco (1972). The Marmousi model was developed to meet several research demands formulated at the time, to which both the model and the data would need to adhere. First, the model should be geologically plausible, preferably based on a real situation. Second, the model should be complex, containing several reflectors, steep dips, and strong velocity gradients in both vertical and lateral directions. Thus, the generated data should be so complex that it would not satisfy the assumptions of conventional processing, requiring depth migration methods to get good images.
The model used in this work is an upgrade of the original Marmousi model, now called Marmousi2, as described in detail by Martin (2004) and Martin et al. (2006). This update consisted of the reconstruction and spatial expansion of the original Marmousi model, preserving its structural complexity. Besides, they added new channels and hydrocarbon traps and kept the original model reservoirs.
The lithological composition of the Marmousi model is predominantly of shale units, with occasional sand layers. The central area has large faults, with an anticline composed of marl (carbonate-rich shale). In addition, this model has an unconformity overlapped by large salt layers separating the marls from the deeper anticlinal units, which are also mostly shales with some interweaving of sand layers (Fig. 3).
The Marmousi model has a series of hydrocarbon reservoirs contained in its structural model, displayed as gas (red) and oil (green) in Figure 3. In this model, all layers containing hydrocarbons are sand. The reservoirs are distributed within a complex fault zone at varying depths, also in the layers with simple structural settings and more structurally complex locations, such as folds.
All figures relate to the complex Marmousi2 used to apply the developed methodology for pore pressure prediction based on porosity. The general vertical exaggeration of the following figures is of the order 4:1, and the interpretation principle is that all computed sections should resemble the original data in Figure 3.
The available data has seismic velocities [V p (x, z) and V s (x, z)] and density [ρ(x, z)] distributions, as shown in Figures 4,  5 and 6, respectively, which are the necessary starting information for calculating stress distribution. Also, another relevant information is the γ(x, z) ratio distribution shown in Figure  7, which participates in the modeling and represents the contribution of the seismic information for the prediction of lowand high-pressure zones and pore pressure.
In the absence of borehole data, we adopted a strategy to create a model with the lithological distinction between the main lithologies as described by Martin et al. (2006) and depicted in Figure 3. Using the information for the velocities V p , V s , and Figure 3. Structural elements, formations, and lithologies of Marmousi2 according to Martin et al. (2006). In general, the vertical exaggeration is of the order of 4:1. 6/12 density ρ, we constructed a model giving a constant value for each type of lithology as described in the caption of Figure 8.
The porosity model used Equation 10 and the matrix of the model with the lithological separation as shown in Figure 8. Using the density model as an input, as displayed in Figure 6, the porosity estimate was applied at all points of density distribution, applying the matrix density values given in Table  1 for each type of lithology. Figure 9 shows a wide range of porosity values, which highlights the large contrast between the sandstone reservoirs and the shales in the region between 500 and 2,000 meters of depth. Figure 10 shows the hydrostatic pressure P h (x, z) corresponding to a linear trend according to Equation 1 where the constant value for water density used was ρ f = 1.01 g/ cm 3 . Figure 11 shows the variation of pressure P z according

Figure 8.
A lithological separation where each integer in the caption corresponds to a rock type in the following order: 1, water (dark blue); 2, shale (blue); 3, sandstone (light blue); 4, sandstone with hydrocarbons (yellow); 5, limestone (orange); 6, salt (red wine color).    Figure 11. Overburden pressure section, P z (x, z), according to Equation 2. We can see some very light variations with respect to Figure 10.
to Equation 2. The outline of map 11 is similar to 10, but it shows some details, as can be seen in the central region, where the complex fault zone is located. We should keep in mind the limitations due to the color scale of the figures. In the other hand, Figure 12 represents a partial conclusion of this work, where shows the rock pressure distribution P rv (x, z) according to our Equation 7. We notice in this figure a very different outline with respect to the previous described pressure sections 10 and 11.
The outline shown in Figure 12 is due to the contribution of lateral variations of horizontal stress, σ xx = σ yy , which depend on = ratio for calculating P rv . A more complete modeling would consider the condition of non-symmetry in σ xx , σ yy , and σ zz stresses.

Porosity profiles and NCT estimation
The next step was to select some ordered porosity profiles for specific analysis, taken from the calculated model displayed in Figure 9, and shown in Figure 13. These profiles are positioned starting at 10,500 m from the left border of the model, and they cross two reservoirs in different depths. The first reservoir is located in the complex central fault region, around 1,000 m depth, and the second is located at the top of the lower anticlinal and around 3,000 m depth. We selected for exemplification only the region with the presence of shale to find the best fitting exponential curve to the porosity values. That is, the purpose of this step is to find the NCT involving the shale formations. Figure 14 shows details of the selected profiles across the main reservoir, where the porosity was calculated using Equation 10. The green line represents the fit in the least-squares sense between Athy (1930)'s exponential model and the red dots (shale formations). The obtained values for the model were ϕ 0 = 0.43367 (fixed porosity for the mudline), and c (compaction factor) is variable (example, c = 0.0006773 m -1 for one line). Furthermore, we can see that Athy's exponential model (NCT) is suitable for a general representation of the variation of porosity with depth. Besides, based on this model, the line segments positioned to the right side of the green reference line indicate the condition of higher pore pressure and, to the left, of lower pore pressure. Also, the reservoir is located at a depth of around 3,000 meters, and we can see that the condition is systematic of higher pressure, but it does not stand out with high values as in other places of the section.

Pore pressure prediction sections and profiles
In the following figures, we are interested in looking for patterns, combinations, correlations, coherences, values, and caption scales, to compare, understand, and interpret the obtained sections and profiles. The combination of all information is the method for interpreting and characterizing the target zones (reservoirs, collectors). Figures 15 and 16 represent the sections of pore pressure calculated, respectively, by Equations 12 and 13 where we can see the good coherence between them. Also, that Figure 16 more clearly displays the details of the target reservoirs, which we consider an improvement and a partial conclusion of this work. Observe that the caption of Figure 15 indicates the presence of low pore pressure and even a strange negative value (see Equation 12) as a result of high contrast for the physical parameters of the salt formation. However, this effect does not appear for the model of the Equation 13.
The effective pressure is considered as a potential driving mechanism for fluid percolation, and a lower pressure zone, capped by a high-pressure zone, is the proper zone for fluid accumulation (Sibiryakov et al. 2020). Figures 17 and 18 represent sections of effective pressure calculated by Equation 8 with its variants 12 and 13, respectively. As in the previous case of Figures 15 and 16, Figure 18 more effectively reveals details with respect to the target reservoirs, and a better color gradient for the figure; besides, Figure 18 represents a partial conclusion in this work by more effectively highlighting the reservoirs. Figures 19 and 20 represent the seven profiles selected for detail analysis, as shown in Figure 13 across a target zone. The curves shown are calculated with their respective models where the hydrostatic, rock, pore, and effective pressures are included.
For the set of profiles in Figure 19, the reference bounds are linear and shown in blue [hydrostatic P h (z)] and red [overburden P z (z)] colors. Furthermore, the pore pressure P p (z) (z) and the effective pressure P eff (z) (z) vary around these reference bounds, but mainly around the hydrostatic trend and in the opposite sense (negative correlation: P eff = P sca -P p ). Also, for the target reservoir located at 3,000 meters depth, the pore pressure is more indicative.   Figure 9. The profiles serve for further analysis of the anticline region that includes the main reservoir (see Fig. 14). The orange lines are placed over the reservoir and the white lines outside and near the reservoir.      (rv) (z)]. Also, for the target reservoir located at 3,000 meters depth, the combination of the rock pressure P rv (z) with the pore pressure P p (rv) (z) is more indicative for the target.
The complementary Figures 21 and 22 depict the deviation between the pore pressure models, (12) and (13), and the hydrostatic model, (1). The first observation is that these findings are very coherent, as can be seen by the color distribution. The second is that the scales of the captions show the double magnitude interval, which is explained by the important presence of the γ ratio that represents the seismic data. Therefore, we can wind up that the modeling is consistent, with these figures showing only a relative value since the difference could be the opposite.
Disclaimer: This is an incomplete paper and, therefore, represents research in progress. It is not meant to represent professional results, but rather an academic opinion of the authors. The authors are not responsible for any practical business applications.

GENERAL CONCLUSIONS
We outlined throughout the text some partial conclusions for the different pressure sections and profiles. In this section, we provide general comments only.
We demonstrated in this work the possibility of using the proposed method that uses seismic information, V p (x, z) and V s (x, z) velocities and density ρ(x, z), to calculate rock, pore, and effective pressure distribution. These are useful for characterizing potential zones for oil and gas accumulation. However, it is necessary to have a lithological description of the basin, and possibly use empirical and petrophysical models. The proposed macro model establishes the use of the rock pressure P rv (x, z) in the models for different pressure calculations, instead of the simple overburden pressure P z (x, z). We can state that the present method has adaptative characteristics in the sense that it can be reviewed as new information is incorporated in the processing.
It is relevant to consider a combination of different modeling for the analysis of the target zones. Therefore, as a macro model, we should consider driving rock pressure P rv (x, z), pore pressure P p (rv) (x, z), and effective pressure P eff (rv) (x, z) to understand the fluid percolation mechanism between highand low-pressure zones in the basin, together with the fundamentals of petroleum geology.
For the application of the method, it is necessary to know the porosity distribution values. Therefore, in the absence of borehole logs, an estimation of this property is made using empirical relationships to produce an optimum porosity model for the basin.
The calculation of the NCT is used for the prediction of pore pressure. Therefore, lithological separation is necessary to identify the geological formations and adjust the NCT for the same type of rock. Namely, an ill-fit curve can create a false over-pressure or under-pressure zone.
A usual geological descriptive mechanism for pore pressure generation is the compaction unbalance, as shown in Figure  1, but this description is physically limited. Furthermore, our physical modeling for the present state can be guided by the principles of petroleum geology concerning the generation and accumulation of oil and gas. In addition, the pore pressure-porosity model employed in this work is closely related to detecting changes in compaction rates (ϕ 0 and c parameters), which are reflected in porosity values.
Our main conclusions in this work are directly related to Figures 12, 16, 18, and the companion profile figures, that are based on the P rv (x, z) pressure model. We can summarize our conclusions based on the calculated images by laying out the following sections: • The geological information highlighting the potential reservoirs and specifying the targets; • The rock pressure to show the low-and high-pressure zones; • The pore pressure P p  Pore pressure deviation section calculated by the difference between the pore pressure from overburden pressure and the hydrostatic model:

Figure 22.
Pore pressure deviation section calculated by the difference between the pore pressure from rock pressure and the hydrostatic model: • Then, display the profiles over the target zones, and focus on the low-pressure zones with relative high pore pressure [P p (rv) (x, z)], when these conditions correspond to low effective pressure [P eff (rv) (x, z)]; We analyze these trends taking the hydrostatic linear trend [P h (x, z)] as reference; • To extend the reservoir target zone, it can be done laterally and vertically.
The combined analysis of these three figures will identify patterns and correlations, outline and characterize the target zones, and make conclusions.
The modeling of 3D data requires more computation power. However, a complete application of this methodology needs 3D data where we can calculate the stress distributions in full, and no symmetry is adopted.

APPENDIX
We used least squares regression to calculate the NCT, shown in Figure Where: ϕ n (z) = the NCT curve; z = the depth below the mudline level; ϕ 0 = the porosity at the mudline; c = the compaction factor, both estimated from the curve fitting.
Then, applying the natural logarithm operator to the observed and model (Equation 14) where y (m) = lnϕ n , a = lnϕ 0 , b = -c. Figure 14 shows y (obs) by red dots and y (m) by a green line.
The misfit between the observed and calculated data is defined by Equation 17.  Depending on the description and classification of the given data, the optimizing functional can have a different expression. In this case, the estimated parameters can be different, but not that different from the least-square fit.