Direct Inversion Method of Brittleness Parameters Based on Reweighted Lp-Norm

: Brittleness is an important factor that indicates shale properties, as well as fracturability, and it can be well-represented using the elasticity parameter E ρ . Seismic inversion allows direct access to the E ρ ofbrittleness parameters. Seismic inversion is a typical ill-posed problem that has an enormous multiplicity of solutions. In order to invert reservoir brittleness parameters more stably and reliably, a direct inversion method for determining brittleness parameters based on a reweighted Lp-norm is proposed, and the reweighted Lp method is introduced to brittleness parameter inversion for the ﬁrst time. The alternating direction method of multipliers (ADMM) is used to establish the inversion structure and to optimize the objective function in blocks, which effectively improves the convergence speed. We ﬁrst introduce a reweighted Lp method and establish a pre-stack inversion objective function based on the reweighted Lp method. Then, theoretical simulation data are applied to compare the inversion outcomes of the new method with those of the traditional method, and the effect of the method in this manuscript is veriﬁed. Finally, the feasibility of this method is further veriﬁed using actual data for experimental analysis. Through an analysis of the experimental results, we ﬁnd that this method can be well-applied to seismic pre-stack inversion calculation and provides a new direct inversion method for the determination of brittleness parameters for exploration geophysics.


Introduction
The brittleness index is a critical factor for indicating the properties of shale and the effect of fracturing.There are two main methods to evaluate brittleness: one is to evaluate brittleness using mineral content or experimental stress results, and the other is to express brittleness with elastic parameters.In terms of elastic parameters, many scholars have conducted studies to evaluate brittleness and have derived brittleness representation formulas [1][2][3].Seismic pre-stack inversion is an effective method of obtaining elastic parameters.The Zoeppritz equation is the basis of pre-stack inversion, but its equations are too complicated and inconvenient to use, so the majority of scholars have approximated it [4,5].At the same time, the Zoeppritz approximation equation provides an opportunity to directly invert brittleness parameters.Zong et al. [6] derived a YPD equation based on the Aki-Richards approximate formula; established a linear relationship between the reflection coefficient and Young's modulus, Poisson's ratio, and density; and provided a direct inversion method for brittleness parameters.Zhang et al. [7] derived approximate equations for compression and transition wave reflection coefficients based on Young's modulus multiplied by density (Eρ), as well as Poisson's ratio and density, providing a new way to obtain direct inversion.After quantitative studies on the lithology and brittleness of shale gas, it was concluded that the Eρ was a better indicator of shale properties [8].Wang [9] found that the Eρ could be used as a more sensitive elastic parameter to calculate the brittleness index and proposed a new formula to calculate brittleness using the Eρ instead of Young's modulus.Zhang et al. [10] proposed a binomial equation based on the relationship between Young's modulus, Poisson's ratio, and density.Ge et al. [11] derived a new three-property parametric reflection coefficient approximation equation for the prediction of the brittleness index.
Seismic inversion is a typical ill-posed problem.In 1943, Tikhonov explained that this problem could be solved, and in 1963, the Tikhonov regularization was proposed [12].The Tikhonov regularization is an additional term.The regularization term is a priori information that transforms an ill-posed problem into a conditionally posed problem to be solved, and the result of its solution is biased toward the prior information, thus mitigating the awkward posing of inverse problems.Rudin et al. [13] proposed total variation (TV) regularization and applied it to image denoising.Wang [14] proposed a P-wave impedance inversion method based on L1 regularization.Zhang et al. [15] proposed a sparse P-wave impedance inversion method with TV regularization.Liu et al. [16] proposed a P-wave impedance inversion method with L1 regularization.Li et al. [17] introduced anisotropic total variational regularization based on the Lp-norm (ATpV) to achieve impedance inversion, which led to smaller inversion errors and improved inversion results.Chen et al. [18] used the Lp-norm to express reflection coefficients and found that the Lp-norm improved the accuracy more than the L1-norm.Wu et al. [19] introduced mixed second-order fractional anisotropic total p-variation regularization (MS_ATpV) to P-wave impedance inversion, which reduced the number of multiple solutions and improved the inversion accuracy.Liu et al. [20] proposed hybrid total variational (HTV) regularization, which achieved better results for inversion.Currently, regularization methods are often used in seismic inversions to mitigate ill-posed problems [21][22][23][24][25][26][27].They can effectively improve the accuracy and stability of seismic inversion.
To better recover the sparsity of seismic inversions, Candes et al. [28] proposed a reweighted L1 method, which could better penalize non-zero coefficients.Cristian et al. [29] described a new method of sparse filtering based on a reweighted L1 algorithm.Ma et al. [30] applied a reweighted L1 method to synthetic aperture radar imaging.Feng et al. [31] proposed a reweighting method to improve the precision of inverse synthetic aperture radar imaging.Song et al. [32] proposed a new edge-preserving image-smoothing algorithm based on a reweighted L1 algorithm.He et al. [33] proposed an impedance inversion method based on a reweighted L1-norm constraint, which could obtain more accurate impedance boundaries and attenuate the pseudo-layer phenomenon.A reweighting algorithm constrains a regularization term by weight, which eliminates the effect of the amplitude of the regularized parametric term on the sparsity, yielding results with better sparsity than results with no weighting.
In this manuscript, we provide a direct inversion method for brittleness parameters based on a reweighted Lp-norm that introduces a reweighted Lp-norm to brittleness parameter inversion for the first time and recovers better result sparsity.The alternating direction method of multipliers (ADMM) is also used for establishing the inversion structure to optimize the objective function in blocks, which accelerates the convergence speed.In the following sections, the theory of the algorithm and the derivation process of the formula are introduced, and the computational effectiveness of the method is tested with theoretical and field data.

Reweighted Lp Method
First, the Zoeppritz approximate equation [7] used in this manuscript is introduced, and its expression is as follows: where θ is the angle of incidence; Eρ is Young's modulus multiplied by density; σ is Poisson's ratio; ρ is density; k is the ratio of shear to compressional wave velocity (k = V S V P ); and R _Eρ , R _σ , and R _ρ are the reflectivity of the Eρ, Poisson's ratio, and density, respectively.To obtain the required parameters from the seismic data, it is necessary to establish the relationship between the seismic data and the parameters.Taking ρ as an example, the reflection coefficient of ρ can be approximated as follows: where R _ρ is represented by the following matrix: where D is the difference matrix, and L ρ is the natural logarithm of density.
According to Equation (1) with the use of a wavelet kernel matrix W(θ), the seismic record S(θ) can be shown as follows: . . .
where L Eρ is the natural logarithm of the Eρ, L σ is the natural logarithm of Poisson's ratio, and L ρ is the natural logarithm of density.The forward equation based on an approximate linear formula can be abbreviated as follows: where B is the matrix consisting of Equations ( 2)-( 4), and L is the matrix of the natural logarithms of the Eρ, Poisson's ratio, and density.As shown in the introduction, we first introduced a weighted Lp-norm into the inversion of brittleness parameters.This method can solve the shortcomings of L1 regularization in inverse objective function.The regularization term is weighted and constrained by the Lp-norm, which yields results with better sparsity and improves the inversion accuracy.The reweighted Lp method for inversion parameters can be defined as follows: where 1 is the L1 − norm, λ is the regularization factor, and w i is the weight value.
where p is the Lp − norm(0 ≤ p < 1), and ε is a constant greater than 0. ε can ensure the stability of the algorithm, that is, if L has a value of zero in a round of iterations, there is a non-zero value in the corresponding position in the next iteration.

Direct Inversion Method of Brittleness Parameters Based on Reweighted Lp
In constructing the objective function, a reweighted Lp method was introduced for the first time to constrain the inverse objective function.Meanwhile, the ADMM method was used to minimize the objective function.According to Equation (8), the traditional objective function can be expressed as follows: where 2 is the L2 − norm, S obs is the observed data, adding the initial model L 0 constraint and the reweighted Lp method constraint to the objective function.The regularization term is a priori information.The ill-posed problem is transformed into a conditionally, wellposed problem for solution.The solution results are biased toward the prior information to reduce the uncertainty of the inversion.The reweighted Lp method contains amplitude information, and the constrained objective function can better restore sparsity.Objective function can be expressed as follows: where µ is the initial model regularization parameter.The matrix R is introduced to replace the matrix DL.The objective function is further shown as follows: Then, the secondary penalty term and dual variable are introduced to change the above objective function into a constrained optimization problem: where η is the difference of the Lagrangian constraint parameters, and C is the dual term.We used the ADMM algorithm to minimize the objective function.Based on the ADMM algorithm, the objective function can be decomposed into the sub-objective functions of L, R, and C, as follows: 16) We used the reweighted Lp method to construct the inversion objective function.The reweighted Lp algorithm is shown in Equations ( 9) and (10), and the inversion objective function is constructed in Equation (14).The ADMM algorithm decomposes the objective function into three sub-objective functions, namely, Equations ( 15)- (17).The inversion results are obtained by solving the sub-objective function.Equations ( 7) and ( 8) are forward problems, and Equation ( 14) is an inverse problem. Figure 1 shows the inversion flow according to which the theoretical and actual data inversion experiments were carried out.
reweighted Lp algorithm is shown in Equations ( 9) and (10), and the inversion objective function is constructed in Equation (14).The ADMM algorithm decomposes the objective function into three sub-objective functions, namely, Equations ( 15), (16), and (17).The inversion results are obtained by solving the sub-objective function.Equations ( 7) and ( 8) are forward problems, and Equation ( 14) is an inverse problem. Figure 1 shows the inver sion flow according to which the theoretical and actual data inversion experiments were carried out.

Experiments
In order to verify the feasibility of the reweighted Lp method, in this section, the reweighted Lp method and the traditional method were each applied to the inversion of the theoretical model and actual data, respectively, for a series of comparative analyses.

Theoretical Model Experiment
The model used in the theoretical experiments was part of the Marmousi2 model.To carry controlled tests, we generated a synthetic dataset that was used as observationa data in the present work.Such a dataset is generated by using the same mechanism to solve the forward problem. Figure 2 shows the angle stack data volume composed of re flection coefficient and wavelet convolution.The observed data made a synthesis accord ing to Equation (7). Figure 2(a), (b), and (c) show the angle superposition data volumes for the incident angles of 10°, 20°, and 30°, respectively.The theoretical models are shown in Figure 3, which includes three parameters: Figure 3

Experiments
In order to verify the feasibility of the reweighted Lp method, in this section, the reweighted Lp method and the traditional method were each applied to the inversion of the theoretical model and actual data, respectively, for a series of comparative analyses.

Theoretical Model Experiment
The model used in the theoretical experiments was part of the Marmousi2 model.To carry controlled tests, we generated a synthetic dataset that was used as observational data in the present work.Such a dataset is generated by using the same mechanism to solve the forward problem. Figure 2 shows the angle stack data volume composed of reflection coefficient and wavelet convolution.The observed data made a synthesis according to Equation (7). Figure 2a-c show the angle superposition data volumes for the incident angles of 10 • , 20 • , and 30 • , respectively.The theoretical models are shown in Figure 3, which includes three parameters: Figure 3a shows the Eρ theoretical model, Figure 3b shows the Poisson's ratio model, and Figure 3c shows the density theoretical model.The models have 800 traces in the cross-sections, and each trace includes 500 sampling points.
We used the traditional method and the weighted Lp method to invert the theoretical models.Figures 4a, 5a and 6a show the inversion results based on the traditional method for the Eρ, Poisson's ratio, and density.Figures 4b, 5b and 6b show the inversion results based on the reweighted Lp method for the Eρ, Poisson's ratio, and density.Compared with the theoretical models, it was demonstrated that the inversion results of the reweighted Lp method had good continuity in the horizontal orientation and clear stratigraphic boundaries in the vertical orientation.The new method could adequately reflect the geological characteristics of the models and was highly consistent with the theoretical models.Compared with the traditional method, the reweighted Lp method showed approximately the same characteristics.However, through careful observation, it was found that the reweighted Lp method had higher accuracy and better stability.In the black rectangular boxes, the reweighted Lp method had better continuity and was closer to the theoretical models, which further proved the significance of the method's improvement.We used the traditional method and the weighted Lp method to invert the theoretical models.We used the traditional method and the weighted Lp method to invert the theoretical models.reflect the geological characteristics of the models and was highly consistent with the theoretical models.Compared with the traditional method, the reweighted Lp method showed approximately the same characteristics.However, through careful observation, it was found that the reweighted Lp method had higher accuracy and better stability.In the black rectangular boxes, the reweighted Lp method had better continuity and was closer to the theoretical models, which further proved the significance of the method's improvement.
(a) (b) In order to reflect the difference between the two methods more clearly, we calculated the absolute error values of inversion for the two methods.Figures 7, 8  reflect the geological characteristics of the models and was highly consistent with the theoretical models.Compared with the traditional method, the reweighted Lp method showed approximately the same characteristics.However, through careful observation, it was found that the reweighted Lp method had higher accuracy and better stability.In the black rectangular boxes, the reweighted Lp method had better continuity and was closer to the theoretical models, which further proved the significance of the method's improvement.
(a) (b) In order to reflect the difference between the two methods more clearly, we calculated the absolute error values of inversion for the two methods.Figures 7, 8, and 9 show the absolute error values of inversion for the two methods.Figures 7(a  reflect the geological characteristics of the models and was highly consistent with the theoretical models.Compared with the traditional method, the reweighted Lp method showed approximately the same characteristics.However, through careful observation, it was found that the reweighted Lp method had higher accuracy and better stability.In the black rectangular boxes, the reweighted Lp method had better continuity and was closer to the theoretical models, which further proved the significance of the method's improvement.
(a) (b) In order to reflect the difference between the two methods more clearly, we calculated the absolute error values of inversion for the two methods.Figures 7, 8, and 9 show the absolute error values of inversion for the two methods.Figures 7(a In order to reflect the difference between the two methods more clearly, we calculated the absolute error values of inversion for the two methods.Figures 7-9 show the absolute error values of inversion for the two methods.Figures 7a, 8a and 9a show the absolute error values of inversion for the traditional approach for the Eρ, Poisson's ratio, and density, respectively.Figures 7b, 8b and 9b show the absolute error values of inversion for the reweighted Lp method for the Eρ, Poisson's ratio, and density, respectively.It can be seen from the three figures that the errors of the inversion results for the proposed method were significantly reduced, and the results were more stable, especially those in the black rectangular boxes.The above experiments show that this research method could improve inversion accuracy and had important significance for improvement.
especially those in the black rectangular boxes.The above experiments show that this research method could improve inversion accuracy and had important significance for iment.
(a) (b) Then, we analyzed the inversion results of the 150th trace and the 500th trace.Figure 10 shows the inversion results of the 150th trace, and Figure 12 shows the inversion results of the 500th trace, where the black, green, blue, and red curves represent the actual model, initial model, inversion outcome of the traditional approach, and inversion outcome of the reweighted Lp method, respectively.Figure 11 shows the absolute error values of inversion for the 150th trace, and Figure 13 shows the absolute error values of inversion for the 500th trace, where the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.From the individual trace inversion outcomes and absolute errors, it can be seen that the line of the inversion outcome of the reweighted Lp method was Then, we analyzed the inversion results of the 150th trace and the 500th trace.Figure 10 shows the inversion results of the 150th trace, and Figure 12 shows the inversion results of the 500th trace, where the black, green, blue, and red curves represent the actual model, initial model, inversion outcome of the traditional approach, and inversion outcome of the reweighted Lp method, respectively.Figure 11 shows the absolute error values of inversion for the 150th trace, and Figure 13 shows the absolute error values of inversion for the 500th trace, where the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.From the individual trace inversion outcomes and absolute errors, it can be seen that the line of the inversion outcome of the reweighted Lp method was Then, we analyzed the inversion results of the 150th trace and the 500th trace.Figure 10 shows the inversion results of the 150th trace, and Figure 12 shows the inversion results of the 500th trace, where the black, green, blue, and red curves represent the actual model, initial model, inversion outcome of the traditional approach, and inversion outcome of the reweighted Lp method, respectively.Figure 11 shows the absolute error values of inversion for the 150th trace, and Figure 13 shows the absolute error values of inversion for the 500th trace, where the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.From the individual trace inversion outcomes and absolute errors, it can be seen that the line of the inversion outcome of the reweighted Lp method was Then, we analyzed the inversion results of the 150th trace and the 500th trace.Figure 10 shows the inversion results of the 150th trace, and Figure 11 shows the inversion results of the 500th trace, where the black, green, blue, and red curves represent the actual model, initial model, inversion outcome of the traditional approach, and inversion outcome of the reweighted Lp method, respectively.Figure 12 shows the absolute error values of inversion for the 150th trace, and Figure 13 shows the absolute error values of inversion for the 500th trace, where the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.From the individual trace inversion outcomes and absolute errors, it can be seen that the line of the inversion outcome of the reweighted Lp method was consistent with the line of the actual model.The method used was close to the actual inversion results, with high accuracy and stability, which further proves the validity of the method.consistent with the line of the actual model.The method used was close to the version results, with high accuracy and stability, which further proves the vali method.In order to carry out the noise resistance experiment, Gaussian white noi added to seismic records, and seismic records with a signal-to-noise ratio of 5 we thesized.Figure 14 shows the angle stack data volume after adding noise, and 14(a), (b), and (c) express the angle stack data volume with incident angles of 10°, 2 30°, respectively.Figures 15, 16, and 17 show the inversion results of the weigh method and the traditional method after adding noise.Figure 15    In order to carry out the noise resistance experiment, Gaussian white noi added to seismic records, and seismic records with a signal-to-noise ratio of 5 we thesized.Figure 14 shows the angle stack data volume after adding noise, and 14(a), (b), and (c) express the angle stack data volume with incident angles of 10°, 2 30°, respectively.Figures 15, 16, and 17 show the inversion results of the weigh method and the traditional method after adding noise.Figure 15  In order to carry out the noise resistance experiment, Gaussian white noise was added to seismic records, and seismic records with a signal-to-noise ratio of 5 were synthesized.Figure 14 shows the angle stack data volume after adding noise, and Figure 14a-c express the angle stack data volume with incident angles of 10 • , 20 • , and 30 • , respectively.Figures 15-17 show the inversion results of the weighted Lp method and the traditional method after adding noise.Figure 15a,b express the Eρ inversion results of the two methods.Figure 16a,b show the Poisson's ratio inversion results of the two methods.Figure 17a,b show the density inversion results of the two methods.It can be seen from the figures that the inversion results of the traditional approach had more outliers, which was more obvious in the black rectangular boxes.The experiments show that, compared with the traditional method, the reweighted Lp method had stronger antinoise ability.We then calculated the absolute error values of inversion for the two methods separately, as shown in Figures 18, 19, and 20. Figure 18(a) and (b) show the  E absolute error values of inversion for the two methods.Figure 19(a) and (b) show the Poisson's ratio absolute error values of inversion for the two methods.Figure 20(a) and (b) show the density absolute error values of inversion for the two methods.From the absolute error plots, it can be seen that the error values of the new method were smaller, with relatively few anomalies near the fault area and higher accuracy.After that, in order to analyze the experimental results for antinoise ability from a microscopic perspective, data for the 150th trace and the 500th trace were extracted to compare the two methods.The inversion results of the two methods are expressed in Figure 21 and Figure 23, and the absolute error values of inversion are expressed in Figure 22 and Figure 24.In Figure 21 and Figure 23, the black, green, blue, and red curves represent the actual model, initial model, inversion outcomes of the traditional approach, and inversion outcomes of the reweighted Lp method, respectively.In Figure 21 and Figure 23, the inversion outcomes of the reweighted Lp method were in better agreement with the actual model.In Figure 22 and Figure 24, the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.The accuracy of the two methods is more clearly shown in Figure 22 and Figure 24, and the inversion results of the new method could maintain the minimum error value and were more stable at most sampling points.The antinoise ability experiments show that the reweighted Lp method had stronger noise resistance and better precision, which further illustrates the significance of the method's improvement.We then calculated the absolute error values of inversion for the two methods separately, as shown in Figures 18-20.Figure 18a,b show the Eρ absolute error values of inversion for the two methods.Figure 19a,b show the Poisson's ratio absolute error values of inversion for the two methods.Figure 20a,b show the density absolute error values of inversion for the two methods.From the absolute error plots, it can be seen that the error values of the new method were smaller, with relatively few anomalies near the fault area and higher accuracy.After that, in order to analyze the experimental results for antinoise ability from a microscopic perspective, data for the 150th trace and the 500th trace were extracted to compare the two methods.The inversion results of the two methods are expressed in Figures 21 and 22, and the absolute error values of inversion are expressed in Figures 23 and 24.In Figures 21 and 22, the black, green, blue, and red curves represent the actual model, initial model, inversion outcomes of the traditional approach, and inversion outcomes of the reweighted Lp method, respectively.In Figures 21 and 22, the inversion outcomes of the reweighted Lp method were in better agreement with the actual model.In Figures 23 and 24, the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.The accuracy of the two methods is more clearly shown in Figures 23 and 24, and the inversion results of the new method could maintain the minimum error value and were more stable at most sampling points.The antinoise ability experiments show that the reweighted Lp method had stronger noise resistance and better precision, which further illustrates the significance of the method's improvement.We then calculated the absolute error values of inversion for the two methods separately, as shown in Figures 18, 19, and 20. Figure 18(a) and (b) show the  E absolute error values of inversion for the two methods.Figure 19(a) and (b) show the Poisson's ratio absolute error values of inversion for the two methods.Figure 20(a) and (b) show the density absolute error values of inversion for the two methods.From the absolute error plots, it can be seen that the error values of the new method were smaller, with relatively few anomalies near the fault area and higher accuracy.After that, in order to analyze the experimental results for antinoise ability from a microscopic perspective, data for the 150th trace and the 500th trace were extracted to compare the two methods.The inversion results of the two methods are expressed in Figure 21 and Figure 23, and the absolute error values of inversion are expressed in Figure 22 and Figure 24.In Figure 21 and Figure 23, the black, green, blue, and red curves represent the actual model, initial model, inversion outcomes of the traditional approach, and inversion outcomes of the reweighted Lp method, respectively.In Figure 21 and Figure 23, the inversion outcomes of the reweighted Lp method were in better agreement with the actual model.In Figure 22 and Figure 24, the blue and red curves represent the absolute error of inversion for the traditional approach and the absolute error of inversion for the reweighted Lp method, respectively.The accuracy of the two methods is more clearly shown in Figure 22 and Figure 24, and the inversion results of the new method could maintain the minimum error value and were more stable at most sampling points.The antinoise ability experiments show that the reweighted Lp method had stronger noise resistance and better precision, which further illustrates the significance of the method's improvement.To consider the effect of the regularization factor λ on the experimental results, dif- ferent values were adjusted for the inversion experiments.The root mean square error (RMSE) of the inversion results was also calculated, and the value of λ was determined when the root mean square error was smallest.In Figure 25, a graph of the variation in   To consider the effect of the regularization factor λ on the experimental results, dif- ferent values were adjusted for the inversion experiments.The root mean square error (RMSE) of the inversion results was also calculated, and the value of λ was determined when the root mean square error was smallest.In Figure 25, a graph of the variation in To consider the effect of the regularization factor λ on the experimental results, different values were adjusted for the inversion experiments.The root square error (RMSE) of the inversion results was also calculated, and the value of λ was determined when the root mean square error was smallest.In Figure 25, a graph of the variation in the RMSE error with A is shown.Figure 26 shows the convergence curve of the objective function.To consider the effect of the regularization factor  on the experimental results, different values were adjusted for the inversion experiments.The root mean square error (RMSE) of the inversion results was also calculated, and the value of  was determined when the root mean square error was smallest.In Figure 25, a graph of the variation in the RMSE error with A is shown.Figure 26 shows the convergence curve of the objective function.

Field Data Experiment
The field data were collected from Sichuan, China.The angle stack data volume is shown in Figure 27, with incidence angles of 9°, 18°, and 27° in Figure 27  To consider the effect of the regularization factor  on the experimental results, different values were adjusted for the inversion experiments.The root mean square error (RMSE) of the inversion results was also calculated, and the value of  was determined when the root mean square error was smallest.In Figure 25, a graph of the variation in the RMSE error with A is shown.Figure 26 shows the convergence curve of the objective function.

Field Data Experiment
The field data were collected from Sichuan, China.The angle stack data volume is shown in Figure 27, with incidence angles of 9°, 18°, and 27° in Figure 27  .

Field Data Experiment
The field data were collected from Sichuan, China.The angle stack data volume is shown in Figure 27, with incidence angles of 9 • , 18 • , and 27 • in Figure 27a-c, respectively.The section had 800 traces with a vertical time length of 750 ms and a sampling interval of 2 ms.The two red lines represent the two layers of geological interpretation, and the black vertical line is the well trajectory.The inversion method was applied to the field data, and Figures 28-30 show the inversion results of the reweighted Lp and the traditional method, respectively.Figure 28a,b express the Eρ inversion results of the two methods.Figure 29a,b show the Poisson's ratio inversion results of the two methods.Figure 30a,b show the density inversion results of the two methods.For the sake of observation, the black rectangular boxes are enlarged in the figures.In these areas, compared with the traditional approach, the reweighted Lp method had better transverse continuity, while the vertical resolution was effectively improved.In addition, the overall effect in all the other regions was greatly improved and could effectively portrayed stratigraphic features.In these areas, compared with the traditional approach, the reweighted Lp method had better transverse continuity, while the vertical resolution was effectively improved.In addition, the overall effect in all the other regions was greatly improved and could effectively portrayed stratigraphic features.For the purpose of verifying the reliability of the inversion method, we compared the logged data with the inversion data of the well location trace.The inversion comparison chart is shown in Figure 31.In order to more obviously reflect the difference between the two methods, we calculated the absolute error values of inversion for the two methods, as shown in Figure 32.As can be seen in the two figures mentioned above, both methods were in good agreement with the logged data.However, when carefully observed, the inversion results of the reweighted Lp method were closer overall to the logged data and were more stable.The black rectangular boxes indicate the reservoir location.In this area, compared with the traditional method, the reweighted Lp method could invert the characteristics of the reservoir more accurately, with less error and higher accuracy.This further authenticated the validity of the method and the significance of the improvement.For the purpose of verifying the reliability of the inversion method, we compared the logged data with the inversion data of the well location trace.The inversion comparison chart is shown in Figure 31.In order to more obviously reflect the difference between the two methods, we calculated the absolute error values of inversion for the two methods, as shown in Figure 32.As can be seen in the two figures mentioned above, both methods were in good agreement with the logged data.However, when carefully observed, the inversion results of the reweighted Lp method were closer overall to the logged data and were more stable.The black rectangular boxes indicate the reservoir location.In this area, compared with the traditional method, the reweighted Lp method could invert the characteristics of the reservoir more accurately, with less error and higher accuracy.This further authenticated the validity of the method and the significance of the improvement.For the purpose of verifying the reliability of the inversion method, we compared the logged data with the inversion data of the well location trace.The inversion comparison chart is shown in Figure 31.In order to more obviously reflect the difference between the two methods, we calculated the absolute error values of inversion for the two methods, as shown in Figure 32.As can be seen in the two figures mentioned above, both methods were in good agreement with the logged data.However, when carefully observed, the inversion results of the reweighted Lp method were closer overall to the logged data and were more stable.The black rectangular boxes indicate the reservoir location.In this area, compared with the traditional method, the reweighted Lp method could invert the characteristics of the reservoir more accurately, with less error and higher accuracy.This further authenticated the validity of the method and the significance of the improvement.

Conclusions
In this research, we proposed a direct inversion method for brittleness parameters based on a reweighted Lp-norm.The regularization constraint in this method could reduce the number of multiple solutions to the inversion problem and could improve inversion stability.The reweighted Lp method effectively restored the sparsity of the inversion results and improved inversion accuracy by constraining the regularization term.In addition, the inversion framework established by the ADMM method improved the convergence speed.
The method was validated by theoretical models and field data.The inversion results expressed that the reweighted Lp method could reduce the inversion error, improve the vertical resolution, enhance the lateral continuity, and describe the reservoir information more accurately.The proposed methodology can be applied in different contexts, both in

Conclusions
In this research, we proposed a direct inversion method for brittleness parameters based on a reweighted Lp-norm.The regularization constraint in this method could reduce the number of multiple solutions to the inversion problem and could improve inversion stability.The reweighted Lp method effectively restored the sparsity of the inversion results and improved inversion accuracy by constraining the regularization term.In addition, the inversion framework established by the ADMM method improved the convergence speed.
The method was validated by theoretical models and field data.The inversion results expressed that the reweighted Lp method could reduce the inversion error, improve the vertical resolution, enhance the lateral continuity, and describe the reservoir information more accurately.The proposed methodology can be applied in different contexts, both in post-stacked data inversion processes [34] and also in inverting pre-stacked data [35], and it provides a new direct inversion method for brittleness parameters for exploration geophysics.
(a) shows the  E theoretical model Figure 3(b) shows the Poisson's ratio model, and Figure 3(c) shows the density theoretica model.The models have 800 traces in the cross-sections, and each trace includes 500 sam pling points.

Figure 4 (
a), 5(a), and 6(a) show the inversion results based on the traditional method for the  E , Poisson's ratio, and density.Figure 4(b), 5(b), and 6(b) show the inversion results based on the reweighted Lp method for the  E , Poisson's ratio, and density.Compared with the theoretical models, it was demonstrated that the inversion results of the reweighted Lp method had good continuity in the horizontal orientation and clear stratigraphic boundaries in the vertical orientation.The new method could adequately
Figure 4(a), 5(a), and 6(a) show the inversion results based on the traditional method for the  E , Poisson's ratio, and density.Figure 4(b), 5(b), and 6(b) show the inversion results based on the reweighted Lp method for the  E , Poisson's ratio, and density.Compared with the theoretical models, it was demonstrated that the inversion results of the reweighted Lp method had good continuity in the horizontal orientation and clear stratigraphic boundaries in the vertical orientation.The new method could adequately
, and 9 show the absolute error values of inversion for the two methods.Figures 7(a), 8(a), and 9 (a) show the absolute error values of inversion for the traditional approach for the  E , Poisson's ratio, and density, respectively.Figures 7(b), 8(b), and 9(b) show the absolute error values of inversion for the reweighted Lp method for the  E , Poisson's ratio, and density, respectively.It can be seen from the three figures that the errors of the inversion results for the proposed method were significantly reduced, and the results were more stable,

Figure 4 .
Figure 4. Eρ inversion results based on (a) traditional method and (b) reweighted Lp method.
), 8(a), and 9 (a) show the absolute error values of inversion for the traditional approach for the  E , Poisson's ratio, and density, respectively.Figures7(b), 8(b), and 9(b) show the absolute error values of inversion for the reweighted Lp method for the  E , Poisson's ratio, and density, respectively.It can be seen from the three figures that the errors of the inversion results for the proposed method were significantly reduced, and the results were more stable,

Figure 5 .
Figure 5. Poisson's ratio inversion results based on (a) traditional method and (b) reweighted Lp method.

Figure 6 .
Figure 6.Density inversion results based on (a) traditional method and (b) reweighted Lp method.

Figure 7 .Figure 8 .Figure 9 .
Figure 7.  E absolute error of inversion based on (a) traditional method and (b) reweighted Lp method.

Figure 9 .
Figure 9. Density absolute error of inversion based on (a) traditional method and (b) reweighted Lp method.

Figure 10 .
Figure 10.Comparison of inversion results for the 150th trace.

Figure 11 .
Figure 11.Absolute error comparison of the 150th inversion results.

Figure 11 .
Figure 11.Absolute error comparison of the 150th inversion results.

Figure 13 .
Figure 13.Absolute error comparison of the 500th inversion results.
(a) and (b) expr  E inversion results of the two methods.Figure 16(a) and (b) show the Poisson inversion results of the two methods.Figure 17(a) and (b) show the density invers

Figure 13 .
Figure 13.Absolute error comparison of the 500th inversion results.
(a) and (b) expr  E inversion results of the two methods.Figure 16(a) and (b) show the Poisson inversion results of the two methods.Figure 17(a) and (b) show the density invers sults of the two methods.It can be seen from the figures that the inversion results

Figure 13 .
Figure 13.Absolute error comparison of the 500th inversion results.

Figure 16 .
Figure 16.Poisson's ratio inversion results after adding noise based on (a) reweighted Lp method and (b) traditional method.Figure 16.Poisson's ratio inversion results after adding noise based on (a) reweighted Lp method and (b) traditional method.

Figure 17 .
Figure 17.Density inversion results after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 18
Figure 18. E inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 17 .
Figure 17.Density inversion results after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 17 .
Figure 17.Density inversion results after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 18
Figure 18. E inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 18 .
Figure 18.Eρ inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 19 .
Figure 19.Poisson's ratio inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 20 ..Figure 21 .
Figure 20.Density inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 19 .Figure 19 .
Figure 19.Poisson's ratio inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 20 ..Figure 21 .
Figure 20.Density inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 20 .
Figure 20.Density inversion error after adding noise based on (a) traditional method and (b) reweighted Lp method.

Figure 19 .
Figure 19.Poisson's ratio inversion error after adding noise based on (a) traditional method reweighted Lp method.

Figure 20 ..Figure 21 .
Figure 20.Density inversion error after adding noise based on (a) traditional method and weighted Lp method.

Figure 21 .
Figure 21.Comparison of inversion results for the 150th trace after adding noise.

Figure 22 .
Figure 22.Comparison of inversion results for the 500th trace after adding noise.

Figure 24 .
Figure 24.Absolute error comparison of the 500th trace inversion results after adding noise.

Figure 22 .
Figure 22.Comparison of inversion results for the 500th trace after adding noise.

Figure 24 .
Figure 24.Absolute error comparison of the 500th trace inversion results after adding noise.

Figure 24 .
Figure 24.Absolute error comparison of the 500th trace inversion results after adding noise.

20 .Figure 24 .
Figure 24.Absolute error comparison of the 500th trace inversion results after adding noise.
Appl.Sci.2022,12,  x FOR PEER REVIEW 16 of 20 respectively.The section had 800 traces with a vertical time length of 750 ms and a sampling interval of 2 ms.The two red lines represent the two layers of geological interpretation, and the black vertical line is the well trajectory.The inversion method was applied to the field data, and Figures28, 29, and 30 show the inversion results of the reweighted Lp and the traditional method, respectively.Figure28(a) and (b) express the  E inversion results of the two methods.
Figure 29(a) and (b) show the Poisson's ratio inversion results of the two methods.
Figure 30(a) and (b) show the density inversion results of the two methods.For the sake of observation, the black rectangular boxes are enlarged in the figures.

Figure 28 .
Figure 28.Eρ inversion results based on (a) traditional method and (b) reweighted Lp method.

Figure 29 .Figure 30 .
Figure 29.Poisson's ratio inversion results based on (a) traditional method and (b) reweighted Lp method.

Figure 30 .
Figure 30.Density inversion results based on (a) traditional method and (b) reweighted Lp method.

20 .Figure 31 .
Figure 31.Inversion results of the trace near the well and the logged well.Figure 31.Inversion results of the trace near the well and the logged well.

Figure 31 ..Figure 31 .Figure 32 .
Figure 31.Inversion results of the trace near the well and the logged well.Figure 31.Inversion results of the trace near the well and the logged well.

Figure 32 .
Figure 32.Absolute error comparison of the inversion results.