Method and program product for determining water flux into serpentinized mantle peridotite

Abstract

The application discloses a method and program product for determining serpentinization water influx of a serpentized mantle peridotite. The method comprises the following steps: obtaining geophysical data of a target area, constructing a geological profile model representing the geological structure of the target area based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model; using a class profile method to inversely reconstruct geometric shape information of a fault in the geological profile model to obtain fault spatial distribution data; constructing a geodynamics model according to the geological profile model, determining seepage data of free water seeping downward from the upper part of the fault based on the geodynamics model; determining serpentinization water influx of a serpentinization reaction zone and a target water influx in the target area according to temperature and pressure field simulation results, mineral stability data of the target area and the seepage data. Therefore, the accuracy and reliability of serpentinization water influx calculation are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological exploration, and in particular to a method for determining water flux of serpentinized mantle peridotite and a program product. BACKGROUND

[0002] Subduction zones are important channels for water circulation in the Earth's interior. Water subduction transport significantly changes the physical properties of the mantle through processes such as fault seepage, mineral dehydration, and serpentinization, and has a profound impact on seismic activity, magma generation, and Earth habitability.

[0003] In related technologies, low wave velocity and high conductivity anomaly regions caused by serpentinization are mainly identified based on seismic tomography, marine magnetotelluric detection and other geophysical methods, and the water flux of the above typical regions can be preliminarily estimated. However, the current technical means mainly perform forward simulation through existing observation data, and the accuracy of the estimated water flux in the subduction zone is poor, and it is also difficult to effectively integrate existing observation constraints. SUMMARY

[0004] The present application provides a method for determining water flux of serpentinized mantle peridotite and a computer program product to solve the technical problem of insufficient identification accuracy of serpentinization water flux in the subduction zone in related technologies.

[0005] According to an aspect of the present application, a method for determining water flux of serpentinized mantle peridotite is provided, which comprises:

[0006] Obtaining geophysical data of a target region, constructing a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure includes at least one of a trench, an outer arc, an ocean basin, an oceanic crust and an upper mantle;

[0007] Reconstructing the geometric information of the fault in the geological profile model using a similar profile method to obtain fault spatial distribution data of the target region; wherein the geometric information includes at least one of fault throw data, dip angle data, spread data and activity direction data;

[0008] Constructing a geodynamic model according to the geological profile model, determining seepage data of free water seeping downward from the upper part of the fault based on the geodynamic model with the fault spatial distribution data as the boundary condition; wherein the seepage data includes seepage path, seepage speed and fault flux distribution;

[0009] The serpentinization water influx of the serpentinization reaction zone is determined according to temperature and pressure field simulation results, mineral stability data of the target region, and the permeation data, and the target water influx in the target region is determined according to the serpentinization water influx of the serpentinization reaction zone, wherein the target water influx includes a transient water influx and / or a cumulative water influx of the serpentinization water influx of the serpentinization reaction zone in the target region.

[0010] According to another aspect of the present application, a device for determining serpentinization mantle peridotite water influx is provided, which comprises:

[0011] A model construction module is configured to acquire geophysical data of a target region, construct a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embed initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure includes at least one of a trench, an outer arc, an ocean basin, an oceanic crust, and an upper mantle.

[0012] An inversion reconstruction module is configured to perform inversion reconstruction on geometric information of a fault in the geological profile model by using a plane section method to obtain fault spatial distribution data of the target region, wherein the geometric information includes at least one of fault throw data, dip angle data, spread data, and activity direction data.

[0013] A permeation data determination module is configured to construct a geodynamic model according to the geological profile model, determine permeation data of free water seeping downward from an upper part of the fault based on the geodynamic model with the fault spatial distribution data as a boundary condition, and wherein the permeation data includes a permeation path, a permeation speed, and a fault flux distribution.

[0014] A serpentinization water influx calculation module is configured to determine a serpentinization water influx of a serpentinization reaction zone according to temperature and pressure field simulation results, mineral stability data of the target region, and the permeation data, and determine a target water influx in the target region according to the serpentinization water influx of the serpentinization reaction zone, wherein the target water influx includes a transient water influx and / or a cumulative water influx of the serpentinization water influx of the serpentinization reaction zone in the target region.

[0015] According to another aspect of the present application, an electronic device is provided, which comprises:

[0016] At least one processor; and a memory connected with the at least one processor in communication; wherein the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor to enable the at least one processor to perform the method for determining water flux into serpentinized mantle peridotite according to any one of the embodiments of the application.

[0017] According to another aspect of the application, there is provided a computer-readable storage medium storing computer instructions for causing a processor to implement the method for determining water flux into serpentinized mantle peridotite according to any one of the embodiments of the application when executed by the processor.

[0018] According to another aspect of the application, the embodiments of the disclosure also provide a computer program product, comprising a computer program which, when executed by a processor, implements the method for determining water flux into serpentinized mantle peridotite according to any one of the embodiments of the disclosure.

[0019] The technical scheme of the embodiment of the present application comprises the following steps: obtaining geophysical data of a target region, constructing a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure comprises at least one of a trench, an outer arc, an ocean basin, an ocean crust and an upper mantle; thereby, by combining different geological structures and the geophysical data of the target region, the target region of different geological structures is simulated truly and comprehensively, so that a geological profile model more in line with the actual situation is obtained. Further, the geometric shape information of a fault in the geological profile model is reconstructed by inversion using a plane section method, so as to obtain fault spatial distribution data of the target region; wherein the geometric shape information comprises at least one of fault throw data, dip angle data, spread data and activity direction data. By reconstructing the geological profile model by inversion using the plane section method, the fault spatial distribution data with higher accuracy can be obtained based on the fine geological shape geometric information in the geological profile model, thereby improving the accuracy of the fault spatial distribution data. Then, an earth dynamics model is constructed according to the geological profile model, the fault spatial distribution data is taken as a boundary condition, and the penetration data of free water penetrating downward from the upper part of the fault is determined based on the earth dynamics model; wherein the penetration data comprises a penetration path, a penetration speed and a fault flux distribution; by limiting and constraining the penetration range of the free water, the penetration of the free water can be better simulated based on the fault spatial distribution data, so as to obtain the penetration data result of the free water penetrating downward from the fault in the target region, thereby enhancing the significance and effectiveness of the penetration data in the target region. Finally, the serpentinization reaction zone water influx of the serpentinization reaction zone is determined according to the temperature and pressure field simulation result, the mineral stability data of the target region and the penetration data, the target water influx in the target region is determined according to the serpentinization reaction zone water influx of the serpentinization reaction zone, wherein the target water influx comprises the instantaneous water influx and / or cumulative water influx of the serpentinization reaction zone water influx of the serpentinization reaction zone in the target region, so as to obtain the serpentinization reaction zone water influx of the serpentinization reaction zone with high precision by combining the geometric constraints and water penetration physical mechanisms of the target region, and further calculate the instantaneous water influx and / or cumulative water influx of the region, so as to more comprehensively evaluate the water influx change of the target region, thereby improving the expansibility and comprehensiveness of the technical scheme. Thus, by the above technical steps, a more accurate and detailed water influx change graph and analysis and evaluation result of the target region are obtained, thereby effectively improving the utilization efficiency of the geophysical data and the accuracy of the water influx calculation of the target region.

[0020] It is to be understood that the description of the background art is not an acknowledgement or any type of suggestion that this art formed part of the common general knowledge of persons having ordinary skill in the art. BRIEF DESCRIPTION OF DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed in the embodiments description. Obviously, the drawings described in the following embodiments are only some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained from these drawings without creative labor.

[0022] Figure 1 is a flow chart of a method for determining water influx of serpentinized mantle peridotite according to the first embodiment of the present application;

[0023] Figure 2 is a flow chart of a method for determining water influx of serpentinized mantle peridotite according to the second embodiment of the present application;

[0024] Figure 3A is a flow chart of a method for determining water influx of serpentinized mantle peridotite according to the third embodiment of the present application;

[0025] Figure 3B is a schematic diagram of a multi-beam topographic map of the Mariana subduction zone which can be used in the third embodiment of the present application according to the third embodiment of the present application;

[0026] Figure 3C is a schematic diagram of a water influx plan view of the Mariana subduction zone which can be used in the third embodiment of the present application according to the third embodiment of the present application;

[0027] Figure 4 is a structural block diagram of a device for determining water influx of serpentinized mantle peridotite according to the fourth embodiment of the present application;

[0028] Figure 5 is a structural schematic diagram of an electronic device for implementing the fifth embodiment of the present application. DETAILED DESCRIPTION

[0029] In order to make the technical personnel in the art better understand the present application scheme, the following will combine the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only some embodiments of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor should be within the scope of protection of the present application.

[0030] It is to be understood that the terminology "target", "instantaneous", "cumulative" and the like used in the detailed description and the appended claims of the present application and the above abstract are used in the sense of the applicants to differentiate between similar objects and are not necessarily intended to describe a particular sequential or chronological order. It is to be understood that the use of such terms as data herein is interchangeable under appropriate circumstances, so that the embodiments of the application described herein can be embodied or carried out in other than the order described or shown herein. Furthermore, the term "comprising" and any of its derivatives is intended to cover the inclusion of the recited elements or steps, but not the exclusion of other elements or steps, e.g., a process, method, system, product or apparatus that comprises a list of steps or elements as limitations is not necessarily limited to only those steps or elements specifically recited, but can include additional steps or elements not expressly listed or inherent to such process, method, product or apparatus.

[0031] It should be noted that the modification of "one", "multiple" mentioned in the present disclosure is illustrative but not restrictive, and those skilled in the art should understand that "one or more" should be understood unless otherwise explicitly indicated in the context.

[0032] The names of the messages or information exchanged between the plurality of devices in the embodiments of the present disclosure are only for illustrative purposes, and are not intended to limit the scope of the messages or information.

[0033] It can be understood that the data involved in the technical solution (including but not limited to data itself, acquisition or use of data) should comply with the requirements of relevant laws, regulations and provisions.

[0034] Embodiment one

[0035] Figure 1 A flowchart of a method for determining water flux of serpentinized mantle peridotite is provided for the first embodiment of the present application. The present embodiment can be applicable to the scene of determining water flux of serpentinization in the subduction zone. The method can be executed by a device for determining water flux of serpentinized mantle peridotite, which can be realized in the form of hardware and / or software. Optionally, the device is realized by an electronic device, which can be a mobile terminal, a PC terminal or a server, etc. As shown in the figure, the method can specifically include: Figure 1

[0036] S110, acquiring geophysical data of a target area, constructing a geological profile model for presenting a geological structure of the target area based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure includes at least one of a trench, an outer arc, an ocean basin, an oceanic crust and an upper mantle.

[0037] ​The target area mainly refers to a subduction zone, which is usually a boundary area where the lithospheric plates of the earth converge. The geophysical data specifically includes at least one of multi-beam sounding data, multi-channel seismic reflection data, wide-angle seismic refraction data, seismic tomography data, and marine magnetotelluric measurement data. The geological profile model refers to a geological profile simulation model used to simulate and display the geological structure of the target area, which is mainly used to generate a two-dimensional or three-dimensional image of the real geological structure of the target area. The two-dimensional grid mode specifically refers to a matrix structure divided into rows and columns. The initial water content refers to the fault zone rock pore water volume data at the beginning of simulation. The rock porosity data refers to the proportion of multiple rock void spaces in the geological structure of the target area, which can determine the size of the rock fluid storage capacity.

[0038] Specifically, by obtaining multiple geophysical data in the target area, a geological profile model of the corresponding geological structure in the target area is constructed according to the multiple geophysical data. The initial water content and rock void data in the target area are input into the constructed geological profile model in a two-dimensional grid mode, thereby obtaining a geological profile model containing geological structure and rock data information. In this application, by combining geological structure and rock data information, a more realistic geological structure data can be simulated to achieve more accurate geological profile modeling of the target area, improving the real accuracy of the geological profile model.

[0039] S120, the geometry information of the fault in the geological profile model is reconstructed by using the class profile method to obtain the spatial distribution data of the fault in the target area; wherein the geometry information includes at least one of fault throw data, dip angle data, spread data and activity direction data.

[0040] The class profile method specifically refers to a method of compressing a three-dimensional geological structure into a two-dimensional plane expression. Specifically, the planar spread (horizontal direction) and profile shape (vertical direction) of the geological body are displayed on a single plane at the same time. The spatial distribution data of the fault specifically can include at least one of nucleation parameters, expansion parameters, penetration depth and permeation channel shape data.

[0041] Specifically, the geometry information of the fault in the geological profile model is reconstructed by using the class profile method to obtain the spatial distribution data of the fault after reconstruction. This technical solution can realize rapid visualization of complex geological structures by using the class profile method and reconstruction, thereby enhancing the distribution of the geometry of the fault in the fault space and effectively improving the accuracy of the display of the spatial distribution data of the fault.

[0042] S130, constructing an earth dynamics model according to the geological profile model, taking the fault spatial distribution data as a boundary condition, and determining, based on the earth dynamics model, penetration data of free water penetrating from the upper part of the fault downward; wherein the penetration data comprises a penetration path, a penetration speed, and a fault flux distribution.

[0043] The earth dynamics model can comprise at least one of an elastic earth dynamics model, an elastic-plastic earth dynamics model, and a viscoelastic earth dynamics model, and specifically, the earth dynamics model is mainly used to simulate and output the penetration data of free water in the fault space under a set constraint condition. The free water refers to seawater entering through the fault.

[0044] Specifically, the target type of earth dynamics model can be constructed according to the geological profile model by using a Lagrangian finite difference method or a finite element method, so as to simulate and calculate the penetration data of free water in the process of penetrating from the upper part of the fault to the lower part of the fault. Thus, the stress-strain field evolution process of free water in the fault penetration due to gravity deflection and geological structure is simulated, and accurate free water penetration data results are obtained.

[0045] S140, determining, according to the temperature and pressure field simulation results, the mineral stability data of the target region, and the penetration data, a serpentinization water influx of a serpentinization reaction zone, and determining a target water influx in the target region according to the serpentinization water influx of the serpentinization reaction zone, wherein the target water influx comprises a transient water influx and / or a cumulative water influx of the serpentinization water influx of the serpentinization reaction zone in the target region.

[0046] The temperature and pressure field simulation results can be simulated by a user in advance, and specifically, the temperature and pressure field simulation results can be obtained by fitting the actually obtained temperature and pressure data in the geological structure. The mineral stability data usually comprise at least one of a heat capacity coefficient, a temperature-pressure phase diagram, a reaction rate constant, an acid-base degree, and an oxygen fugacity. The serpentinization reaction zone refers to a region where free water reacts with rocks mainly composed of serpentine minerals. The serpentinization water influx specifically refers to water stored by rocks mainly composed of serpentine minerals after the reaction of free water with the rocks. The serpentinization water influx can refer to a vertical water flow per unit area through the reaction zone per unit time. The transient water influx of the serpentinization water influx specifically refers to the amount of water stored by the rocks at a certain moment, or a vertical water flow per unit area through the reaction zone at a certain moment. The cumulative water influx of the serpentinization water influx specifically refers to the amount of water continuously stored by the rocks, or a vertical water flow per unit area through the reaction zone in a certain period of time.

[0047] Specifically, according to the obtained temperature and pressure field simulation results, mineral stability data of the target area and permeation data, the serpentinization water inflow of the serpentinization reaction zone is determined. Then, by combining the rock density, the unit volume of the serpentinization reaction zone and the water content change per unit volume, the serpentinization water inflow of each calculation unit in the serpentinization reaction zone is further calculated according to the calculation requirement, so as to output the instantaneous water inflow and / or cumulative water inflow of the serpentinization water inflow in the serpentinization reaction zone. The technical scheme determines the high-precision serpentinization water inflow, instantaneous water inflow and / or cumulative water inflow by combining the geometric and geological structure constraints of the target area and the water permeation physical mechanism, so as to facilitate the user to more conveniently and simply evaluate the water inflow change of the target area.

[0048] In one embodiment, the serpentinization water inflow of the serpentinization reaction zone is determined according to the temperature and pressure field simulation results, the mineral stability data of the target area and the permeation data, which comprises: determining the serpentine stable zone of the target area according to the temperature and pressure field simulation results and the mineral stability data of the target area, determining the serpentinization reaction zone according to the preset temperature and pressure window, the serpentine stable zone and the permeation data, and determining the serpentinization water inflow of the serpentinization reaction zone according to the water and peridotite reaction kinetics model.

[0049] The serpentine stable zone refers to the physical and chemical condition space in which the serpentine mineral can exist for a long time without decomposition or conversion into other minerals. The preset temperature and pressure window refers to the temperature and pressure range interval in which a specific geological process can effectively occur. The peridotite reaction kinetics model refers to the rate equation for describing the serpentinization reaction of water and olivine.

[0050] Specifically, the serpentine stable zone in the target area is determined by the temperature and pressure field simulation results and the mineral stability data of the target area. Then, the serpentinization reaction zone is determined according to the preset temperature and pressure window, the serpentine stable zone and the permeation parameter. Finally, the serpentinization water inflow of the serpentinization reaction zone is calculated according to the rate equation of the serpentinization reaction corresponding to the water and peridotite reaction kinetics model. Thus, the water amount stored in the serpentinization reaction zone or the water amount passing through per unit area per unit time is quickly obtained by using the peridotite reaction kinetics model, thereby improving the calculation efficiency of the serpentinization water inflow.

[0051] In one embodiment, the determining the serpentinization water influx of the serpentinization reaction zone according to the water and peridotite reaction kinetics model comprises: constructing a serpentinization rate model according to the water and peridotite reaction kinetics model, and determining the serpentinization water influx of the serpentinization reaction zone according to the rock density of the serpentinization reaction zone, the reaction volume, the water content change per unit volume, and the serpentinization rate model.

[0052] The serpentinization rate model refers to a serpentinization reaction progress spatiotemporal distribution function constructed based on the kinetics model. The rock density specifically refers to the mass per unit volume of unaltered protolith in the reaction zone. The reaction volume specifically refers to the rock volume change rate caused by serpentinization. The water content change per unit volume specifically refers to the water consumption per unit volume of rock in the serpentinization reaction.

[0053] Specifically, the serpentinization rate model is constructed by the water and peridotite reaction kinetics model. The water amount stored in the serpentinization reaction zone or the vertical water flux per unit area per unit time is determined by inputting the rock density, the reaction volume, and the water content change per unit volume of the serpentinization reaction zone into the constructed serpentinization rate model through multi-parameter coupling calculation. The theoretical parameters are combined with the various theoretical parameters of the serpentinization reaction zone, thereby realizing the organic connection between the theoretical parameters and the actual application. The serpentinization water influx calculated by the serpentinization rate model is more consistent with the actual situation, the calculation efficiency and accuracy of the serpentinization water influx are improved, and the deviation of the model calculation result is reduced.

[0054] In one embodiment, after the determining the serpentinization water influx of the serpentinization reaction zone according to the temperature and pressure field simulation results, the mineral stability data of the target region, and the permeation data, the method further comprises: for an observation region containing seismic tomography in the target region, processing the seismic tomography data according to the Voigt-Reuss-Hill average method and the corresponding relationship between P-wave velocity and water content established by experiments, converting the measured P-wave velocity distribution into a water content spatial distribution map, and updating the serpentinization water influx of the observation region according to the water content spatial distribution map.

[0055] The seismic tomography refers to a geophysical imaging technology for inverting the velocity structure or other physical property parameter distribution of underground medium by using seismic wave data. The Voigt-Reuss-Hill average method is a classical method for estimating the equivalent elastic modulus of a multiphase composite material. The corresponding relationship between P-wave velocity and water content specifically refers to the empirical results or theoretical formula between seismic wave velocity and water content established by experiments. The water content spatial distribution map refers to the spatial distribution results of water content in rock in the form of two-dimensional or three-dimensional data volume.

[0056] Specifically, by analyzing the observation area containing the seismic tomography, the Voigt-Reuss-Hill average method is used, and the corresponding relationship between the P-wave velocity and the water content established by the experiment is used to analyze and process the seismic tomography data in the observation area. According to the observed seismic P-wave velocity, the corresponding water content spatial distribution map of the observation area is determined. Further, according to the determined water content spatial distribution map, the vertical water flow per unit area per unit time through the serpentinization reaction zone of the observation area is updated. Thus, the seismic characteristics of different target areas can be better considered, and more accurate water content calculation can be performed on the observation area containing the seismic tomography data, so as to further correct the serpentinization water inflow flux result of the target observation area in combination with the specific geological characteristics, so as to obtain a more accurate serpentinization water inflow flux result.

[0057] The technical scheme of the embodiment of the present application obtains geophysical data of a target region, constructs a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embeds initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure includes at least one of a trench, an outer arc, an ocean basin, an ocean crust and an upper mantle; thus, by combining different geological structures and the geophysical data of the target region, the target region of different geological structures is simulated truly and comprehensively, so that a geological profile model more in line with actual conditions is obtained. Further, the geometry information of a fault in the geological profile model is reconstructed by inversion using a plane section method, so as to obtain fault spatial distribution data of the target region; wherein the geometry information includes at least one of fault throw data, dip angle data, spread data and activity direction data. By reconstructing the geological profile model by inversion using the plane section method, the fault spatial distribution data with higher accuracy can be obtained based on the fine geological morphology geometry information in the geological profile model, so that the accuracy of the fault spatial distribution data is improved. Then, an earth dynamics model is constructed according to the geological profile model, the fault spatial distribution data is taken as a boundary condition, and the penetration data of free water penetrating from the upper part of the fault downward is determined based on the earth dynamics model; wherein the penetration data includes a penetration path, a penetration speed and a fault flux distribution; by limiting and constraining the penetration range of the free water, the penetration of the free water can be better simulated based on the fault spatial distribution data, so that the penetration data result of the free water penetrating from the upper part of the fault downward in the target region is obtained, and the significance and effectiveness of the penetration data in the target region are enhanced. Finally, the serpentinization reaction zone water influx is determined according to the temperature and pressure field simulation result, the mineral stability data of the target region and the penetration data, and the target water influx in the target region is determined according to the serpentinization reaction zone water influx, wherein the target water influx includes the instantaneous water influx and / or cumulative water influx of the serpentinization reaction zone water influx in the target region, so as to obtain the serpentinization reaction zone water influx with high precision by combining the geometric constraints and water penetration physical mechanism of the target region, and further calculate the instantaneous water influx and / or cumulative water influx of the target region, so as to more comprehensively evaluate the water influx change of the target region, and improve the expansibility and comprehensiveness of the technical scheme. Thus, through the above technical steps, a more accurate and fine water influx change graph and analysis and evaluation result of the target region are obtained, and the utilization efficiency of the geophysical data and the accuracy of the water influx calculation of the target region are effectively improved.

[0058] Embodiment two

[0059] Figure 2A flowchart of a method for determining water flux of serpentinized mantle peridotite provided in Embodiment Two of the present application, the scheme in the present embodiment is based on the above-mentioned embodiments, and is refined to the technical scheme of determining the permeation data of free water penetrating downward from the upper part of the fault based on the geodynamic model. Optionally, the permeation data of free water penetrating downward from the upper part of the fault based on the geodynamic model with the fault spatial distribution data as the boundary condition comprises: performing yield stress field inversion based on the geodynamic model with the fault spatial distribution data as the boundary condition to obtain the yield stress distribution information of the fault; and determining the permeation data of free water penetrating downward from the upper part of the fault based on the fault spatial distribution data and the yield stress distribution information of the fault. The specific implementation can be seen from the description of the present embodiment. Among them, the same or similar technical features as the foregoing embodiments will not be described here. As shown in Figure 2 The method can specifically include:

[0060] S210, acquiring geophysical data of a target region, constructing a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure comprises at least one of a trench, an outer arc, an ocean basin, an oceanic crust and an upper mantle.

[0061] S220, reconstructing geometric information of a fault in the geological profile model by using a class profile method to obtain fault spatial distribution data of the target region, wherein the geometric information comprises at least one of fault throw data, dip angle data, distribution data and activity direction data.

[0062] S230, constructing a geodynamic model according to the geological profile model, and performing yield stress field inversion based on the geodynamic model with the fault spatial distribution data as the boundary condition to obtain yield stress distribution information of the fault.

[0063] Wherein, the yield stress field inversion refers to the process of backstepping the yield stress distribution of underground rock according to the geological data and the mechanical parameter condition. The yield stress refers to the stress value when the rock starts to deform plastically. The yield stress distribution information refers to the distribution information of the yield stress in the fault space.

[0064] Specifically, the geodynamic model is constructed by a geological profile, and the fault spatial distribution data is taken as the calculation range of the model, so that the geodynamic model is used to inversely calculate the yield stress distribution of the rock in the fault space based on preset mechanical parameters, to obtain the yield stress distribution information at different positions in the fault space. Thus, by calculating the yield stress distribution information at different positions in the fault space, the activity of the fault, and the stability and destructiveness of the fault can be evaluated.

[0065] In an embodiment, the yield stress distribution information at least includes a yield stress state. Accordingly, the yield stress field inversion based on the geodynamic model and taking the fault spatial distribution data as the boundary condition to obtain the yield stress distribution information of the fault includes: using the Mohr-Coulomb yield criterion, setting the shear strength and the internal friction angle, and using the stress field output by the geodynamic model to inversely calculate the yield stress state of the fault to obtain the yield stress distribution information of the fault.

[0066] The Mohr-Coulomb yield criterion is a mechanical criterion for describing the yield (failure) of materials (such as rock and soil) under shear stress. The shear strength is the ability of the material to resist shear failure, that is, the maximum stress value of the material under shear stress without failure. The internal friction angle is an angle describing the friction characteristics between the particles in the material, and is used to reflect the ability of the material to resist relative sliding under shear stress. The stress field is the basic data for inversely calculating the yield stress state of the fault.

[0067] Specifically, the Mohr-Coulomb yield criterion, the set shear strength, and the internal friction angle are input into the geodynamic model. Thus, the stress field of the rock material is output, and the yield stress state at different positions of the fault is calculated according to the output stress field data and the Mohr-Coulomb criterion. Further, the effective permeability distribution data of the fault along the depth direction can be determined according to the fault spatial distribution data and the yield stress state of the fault, and the permeation direction of the free water along the depth direction of the fault is controlled by introducing an anisotropic tensor. Through the above technical steps, the yield stress state of the target fault can be effectively determined, and the permeability of the free water at different positions of the fault is obtained by combining the fault space and the introduced anisotropic tensor, so that the flow path of the free water in the rock can be more accurately predicted.

[0068] In the embodiment of the present application, the yield stress distribution information of the fault can be determined in the yield state by combining one or more of the parameters such as viscosity softening, elastic modulus, Poisson's ratio, and rheological parameters.

[0069] S240, determining the permeation data of the free water penetrating from the upper part of the fault downward according to the fault spatial distribution data and the yield stress distribution information of the fault.

[0070] Specifically, the permeation path, the permeation speed and the fault flux distribution of the free water permeating from the upper part of the fault to the lower part of the fault are determined by the fault spatial distribution data and the yield stress distribution information of the fault.

[0071] In one embodiment, the permeation data of the free water permeating from the upper part of the fault to the lower part of the fault is determined according to the fault spatial distribution data and the yield stress distribution information of the fault, comprising: determining the effective permeability distribution data of the fault along the depth direction according to the fault spatial distribution data and the yield stress distribution information of the fault, and introducing an anisotropic tensor to control the permeation direction; constructing a multi-dimensional seepage water path model evolving with time based on the Darcy's law and the mass conservation control equation of the fault, and constructing a path constraint condition according to the fault spatial distribution data of the fault; wherein the fault spatial distribution data comprises nucleation parameters, expansion parameters, penetration depth and permeation channel morphology; the permeation data of the free water permeating from the upper part of the fault to the lower part of the fault is simulated according to the initial water content, the maximum permeation depth, the permeation direction, the permeability and stress coupling data, the fluid driving pressure difference and lithology control data, the path constraint condition and the multi-dimensional seepage water path model.

[0072] The effective permeability distribution data refers to the direction-dependent permeation ability of the free water combined with stress-fracture data coupling. The anisotropic tensor specifically refers to a matrix used to describe the difference in the permeation rate direction of the free water. The multi-dimensional seepage water path model specifically refers to a partial differential equation set coupling the Darcy's law and the mass conservation, which is used to solve the spatio-temporal evolution results of pressure and saturation. The path constraint condition refers to the forced permeation boundary of the free water determined according to the fault spatial distribution data. The maximum permeation depth refers to the limit depth of the free water permeation under the condition of the set temperature and pressure. The permeability and stress coupling data are used to represent the empirical function of the effective stress changing the fracture opening degree. The fluid driving pressure difference specifically refers to the pressure difference generated by the fluid flow in the upper and lower ends of the fault. The lithology control data usually refers to the data directly related to the type of rock and its physical, chemical and mechanical properties.

[0073] Specifically, first, according to the fault spatial distribution data and the yield stress distribution information of the fault, the effective permeability distribution data of the fault along the depth direction is determined, and the anisotropic tensor is introduced to control the permeation direction of the free water along the depth direction of the fault. Second, a multi-dimensional seepage water path model evolving with time is constructed based on the Darcy's law and the mass conservation control equation of the fault, and the path constraint condition in the free water permeation process is determined according to the fault spatial distribution data, wherein the fault spatial distribution data includes nucleation parameters, expansion parameters, penetration depth and permeation channel morphology. Finally, the initial water content in the fault rock, the maximum penetration depth, the penetration direction, the stress coupling data of the permeability, the fluid driving pressure difference and the lithology control data, and the path constraint condition are input into the multi-dimensional seepage water path model, so as to obtain the permeation data of the free water from the upper part to the lower part of the fault by combining the influence factors in the free water permeation process and using the multi-dimensional seepage water path model. Through the above technical steps, the coupling of various data can be effectively combined, and the special properties of the fault structure are considered, so that more accurate and real free water permeation simulation results are obtained by using the multi-dimensional seepage water path model, and the prediction accuracy of the free water permeation mode is improved.

[0074] In one embodiment, optionally, the geological structure includes an outer arc region; after the permeation data of the free water from the upper part to the lower part of the fault is simulated according to the initial water content, the maximum penetration depth, the penetration direction, the stress coupling data of the permeability, the fluid driving pressure difference and the lithology control data, the path constraint condition and the multi-dimensional seepage water path model, it further includes: in the outer arc region, determining the permeation data of the fracture based on the mass conservation control equation of the fracture of Darcy's law, and updating the permeation data of the free water from the upper part to the lower part of the fault according to the permeation data of the fracture.

[0075] Wherein, the outer arc region is part of the island arc structure, and the island arc can be generally divided into two units of outer arc and inner arc, and the outer arc is located on the side of the inner arc close to the trench. The mass conservation control equation of the fracture of Darcy's law is a mathematical expression that satisfies Darcy's law and mass conservation law at the same time when fluid flows in the fracture medium, which is used to represent the relationship between pressure distribution and fluid flow in the fracture.

[0076] Specifically, in the outer arc area of the geological structure, by substituting Darcy's law into the mass conservation control equation, a partial differential equation describing the pressure distribution in the fracture and the fluid flow is determined, and the permeation data of the fracture is calculated. Then, the permeation data of the fracture is input into the multi-dimensional seepage water path model, so as to update the permeation data output by the multi-dimensional seepage water path model, and the output result of the model is corrected. The output result of the multi-dimensional seepage water path model is corrected and optimized by the actual rock fracture, so that the motion track of the free water and other fluids can be more accurately predicted, and the prediction accuracy of the fluid permeation data is effectively improved.

[0077] S250, determine the serpentinization water influx of the serpentinization reaction zone according to the temperature and pressure field simulation result, the mineral stability data of the target area and the permeation data, and determine the target water influx in the target area according to the serpentinization water influx of the serpentinization reaction zone, wherein the target water influx includes the instantaneous water influx and / or cumulative water influx of the serpentinization water influx of the serpentinization reaction zone in the target area.

[0078] The technical scheme of the embodiment of the present application can obtain the yield stress distribution information of the fault by taking the fault spatial distribution data as the boundary condition and performing yield stress field inversion based on the geodynamic model, so that the stress distribution in the target fault can be analyzed according to the fault spatial distribution data to determine the stability and destructiveness of the fault. Furthermore, the permeation data of the free water penetrating downward from the upper part of the fault can be determined according to the fault spatial distribution data and the yield stress distribution information of the fault, so that the influence of the free water on the fault during the penetration process of the free water in the fault can be determined, and the flow and penetration information of the free water in the fault can be determined according to the fault spatial distribution data and the yield stress distribution information of the fault, so that more real and effective free water permeation data can be obtained by combining the penetration influence of the free water on the fault. In summary, the flow and penetration information of the free water in the fault can be more comprehensively determined by combining the free water in the limited fault space and the yield stress distribution information of the fault in the technical scheme, so that the prediction accuracy of the free water permeation data is improved.

[0079] Embodiment three

[0080] The embodiment three of the present application provides a flow chart of a method for determining serpentinization mantle peridotite water influx, in order to better introduce the technical scheme provided by the embodiment of the present application, as shown in Figure 3A 、 Figure 3B and Figure 3C , wherein the same or similar technical features as the foregoing embodiments are not described again. The specific implementation includes:

[0081] 1. Using the existing high-resolution geophysical data (e.g. multi-beam bathymetry, multi-channel seismic reflection, wide-angle seismic refraction, seismic tomography and marine magnetotelluric measurement) in the study area, a two-dimensional or quasi-three-dimensional geologically realistic profile model is established, which presents the trench-outer arc-ocean basin-oceanic crust-upper mantle system.

[0082] 2. Using the method of quasi-planar section, the geometry of normal faults (fault throw, dip angle, distribution, and active direction) is restored and quantified, and a fault distribution database (i.e. spatial distribution data of faults) is established, which provides structural boundary conditions for subsequent dynamic calculations.

[0083] 3. The initial water content and rock porosity distribution are labeled synchronously, providing initial values for water migration and reaction calculations.

[0084] 4. Based on the finite difference or finite element method, an elastic, elastic-plastic or viscoelastic geodynamic model is constructed according to the geological profile model, simulating the stress-strain field evolution process of faults in the subduction zone caused by gravity deflection and tectonic loading.

[0085] 5. Using the fault structure reconstructed by the quasi-planar section method as the geometric boundary, the yield stress field is inverted by coupling the Mohr-Coulomb failure criterion, viscosity softening and rheological parameters (density, Poisson's ratio, viscosity) and other mechanical parameters.

[0086] 6. According to the spatial distribution of faults and the yield stress state output by dynamic simulation, the effective permeability distribution of faults along the depth direction is calculated , and an anisotropic tensor is introduced to control the directional change of permeability.

[0087] 7. The nucleation, expansion, penetration depth and permeation channel morphology of the deflection-induced fault are dynamically tracked to determine whether it can enter the serpentine stable zone, providing physical mechanism constraints for subsequent water entry path determination.

[0088] 8. Based on Darcy's law and the mass conservation control equation, a multi-dimensional seepage water entry path model is constructed that evolves with time: .

[0089] 9. Considering the initial water content, porosity evolution, permeability-stress coupling, fluid driving pressure difference and lithology control, the path, speed and flux distribution of free water permeating from the upper part of the fault are simulated.

[0090] 10. In the outer arc area of the subduction zone, the normal fault system caused by plate bending usually has a complex permeability structure, and the fault zone not only includes the main fracture surface that penetrates, but also considers the crack network, and the mathematical control equation is the mass conservation of the crack system:

[0091] ;

[0092] where: fracture porosity, fracture permeability, fracture pressure, intersection coupling coefficient, source-sink term of fracture.

[0093] 11. Combined with the simulation results of temperature-pressure field and mineral stability data, as well as the path and speed of water infiltration from the upper part of the fault, the spatial distribution of the serpentinization stable domain is determined, and the serpentinization rate model is constructed by coupling the water-olivine reaction kinetics. .

[0094] 12. Combined with the serpentinization rate and other parameters, the following serpentinization water flux estimation equation is used:

[0095] ;

[0096] wherein: the rock density is represented; the reaction unit volume is represented; the water content change per unit volume is represented; the serpentinization rate function is represented.

[0097] 13. In the observation constraint area, the Voigt-Reuss-Hill average method is called, and based on the relationship between the P-wave velocity change and the degree of serpentinization measured by experiments, the spatial distribution of water content is quantitatively inverted.

[0098] 14. The multi-dimensional serpentinization area water content map and the cumulative water flux profile are output.

[0099] It should be noted that the above steps are the main implementation steps of the technical solution, which can be referred to as the method flowchart shown in Figure 1 . On the other hand, the steps of the present technical solution are described in more detail by combining actual examples. The specific implementation process of the example is as follows:

[0100] 1. In the 18°-20°N area of the Mariana arc, the multi-source high-resolution geophysical data called include multi-beam sounding data (as shown in Figure 3B , which represents the multi-beam topography of the Mariana subduction zone), multi-channel seismic reflection profile, wide-angle seismic profile, seismic tomography P-wave velocity profile and marine magnetotelluric (MT) profile. Based on the above data, a two-dimensional geological profile is constructed along 18.5°N, and sedimentary, basalt layer, gabbro layer, lower crust layer and upper mantle olivine layer are divided.

[0101] 2. Identify 6-8 normal faults using MCS profile, and restore their geometric parameters by the method of quasi-equilibrium profile: fault displacement: 0.8-2.5 km; dip angle: 35-50°; extension depth: the longest fault penetrates the oceanic crust bottom boundary and enters the upper mantle about 3 km; the sliding direction is consistent with the stress direction (nearly perpendicular to the strike of the subduction zone). Establish a structural database for all fault data.

[0102] 3. Embed the porosity and water content field into the model initial value in a two-dimensional grid manner. The initial porosity of the sedimentary layer is 0.35, the oceanic crust is 0.1-0.15, the lower crust and peridotite is 0.03-0.05; the initial water content is 2.0 wt% for the oceanic crust and 0.2 wt% for the upper mantle.

[0103] 4. Use finite difference to build a two-dimensional elastic-plastic geodynamic model to calculate the bending and stress response process of the subducting plate.

[0104] 5. Use the Mohr-Coulomb yield criterion to determine whether the fault has yielded by setting the shear strength and internal friction angle using the stress field output by the model; determine the maximum penetration depth, and some faults penetrate the oceanic crust bottom boundary and contact the upper mantle.

[0105] 6. Superimpose the fault morphology and stress field obtained by inversion to screen out areas with high tensile stress, fault penetration depth , and permeability threshold as the main seepage path to determine the channel into the serpentine stable domain.

[0106] 7. According to the stress distribution obtained by simulation, set the permeability field according to the following relationship:

[0107] ;

[0108] where: ; ; is the fault activity (initially 1.0, and decreases to 0.2 after closure).

[0109] 8. Use the finite element method to solve: . The initial value of the P field is controlled by the hydrostatic pressure and the plate loading pressure.

[0110] 9. Pore evolution and lithology coupling control. The pore evolution depends on the local stress, and the porosity in the tensile stress area increases linearly with time ( ); the permeability is coupled with the porosity: ; the influence of lithology on permeability is adjusted by adjusting the coefficient.

[0111] 10. In the outer arc region of the subduction zone, the fracture mesh is considered, and the seepage exchange in the fracture domain is realized by solving the mathematical control equation of the mass conservation of the fracture system. Exemplarily, the following is taken for the Mariana Trench ; .

[0112] 11. On the basis of the simulation of the geodynamics and thermal structure of the Earth, the temperature-pressure (T-P) distribution field in the subducted plate is calculated, and the serpentinization stable zone (usually corresponding to a temperature of 300-600°C) is determined according to the mineral phase stability theory. In combination with the simulation results of the seepage path, the seepage volume unit entering the temperature-pressure window is identified as the serpentinization reaction active zone, and the reaction kinetics model of the olivine-water system is introduced. The reaction rate is constructed by using an Arrhenius-type expression, that is, , wherein is an empirical coefficient, is an activation energy, is a gas constant, and the function is used to quantify the serpentine generation rate under different conditions and is used to couple the water flux estimation.

[0113] 12. On the basis of the reaction zone determination in step 11, a unit water inflow flux estimation model based on the serpentinization rate function is established, and the following expression is used:

[0114] ; wherein, rock density; reaction unit volume; water content change per unit volume; is the temperature-pressure controlled reaction rate function defined in step 11. Through point-by-point flux calculation of each calculation unit in the reaction zone, the instantaneous and cumulative water inflow flux maps of the serpentinization zone are accumulated and output. (As shown in Figure 3C , which represents the plan view of the water inflow flux in the Mariana subduction zone).

[0115] 13. To verify the reliability of the serpentinization reaction and water inflow flux in the simulation, the Voigt-Reuss-Hill average method is used in the seismic tomography observation area, and the P-wave velocity-serpentine content relationship obtained by experiments is applied to the inversion of the water content spatial map from the measured P-wave velocity distribution. The typical mapping relationship is that when Vp is 7.5 km / s, the water content is <1 wt%; when Vp is 6.5 km / s, the water content is about 4 wt%; and when Vp is 6.1 km / s, the water content is more than 6 wt%. The inversion results are compared with the water content distribution calculated in step 12, and the spatial consistency and physical reasonableness of the water inflow prediction of the model are verified.

[0116] The technical scheme restores the formation process of the flexural normal fault and realizes the modeling of the temporal and spatial evolution of the permeability, and can truly describe the control ability of the fault on the water inlet path and intensity at different stages. The full physical field coupling modeling of mechanical stress, heat conduction, seepage migration and the like is realized, and the scientific credibility of the water flux calculation is significantly improved.

[0117] Embodiment four

[0118] Figure 4 A structural schematic diagram of a device for determining serpentinization mantle peridotite water inlet flux provided by the fourth embodiment of the present application. The device can be realized by software and / or hardware, and can be configured in an electronic device. As shown in the figure, the device comprises a model construction module 401, an inversion reconstruction module 402, a permeation data determination module 403 and a serpentinization water inlet flux calculation module 404. Figure 4

[0119] The model construction module 401 is configured to obtain geophysical data of a target region, construct a geological profile model for presenting a geological structure of the target region based on the geophysical data, and embed initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure comprises at least one of a trench, an outer arc, an ocean basin, an oceanic crust and an upper mantle; the inversion reconstruction module 402 is configured to perform inversion reconstruction on geometric information of a fault in the geological profile model by using a similar cross-section method to obtain fault spatial distribution data of the target region; the geometric information comprises at least one of fault throw data, dip angle data, distribution data and activity direction data; the permeation data determination module 403 is configured to construct a geodynamic model according to the geological profile model, determine permeation data of free water seeping downward from the upper part of the fault based on the geodynamic model with the fault spatial distribution data as a boundary condition; the permeation data comprises a permeation path, a permeation speed and a fault flux distribution; and the serpentinization water inlet flux calculation module 404 is configured to determine a serpentinization water inlet flux of a serpentinization reaction zone according to a temperature and pressure field simulation result, mineral stability data of the target region and the permeation data, and determine a target water inlet flux in the target region according to the serpentinization water inlet flux of the serpentinization reaction zone, wherein the target water inlet flux comprises a transient water inlet flux and / or a cumulative water inlet flux of the serpentinization water inlet flux of the serpentinization reaction zone in the target region.

[0120] ​The technical scheme of the embodiment of the present application, the model construction module 401 obtains geophysical data of a target area, constructs a geological profile model for presenting a geological structure of the target area based on the geophysical data, and embeds initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure includes at least one of a trench, an outer arc, an ocean basin, an ocean crust and an upper mantle; thereby, by combining different geological structures and the geophysical data of the target area, the target area of different geological structures is simulated truly and comprehensively, so that a geological profile model more in line with actual conditions is obtained. Further, the inversion reconstruction module 402 reconstructs the geometric shape information of the fault in the geological profile model by using the class plane section method to obtain the spatial distribution data of the fault in the target area; wherein the geometric shape information includes at least one of fault throw data, dip angle data, spread data and activity direction data. By using the class plane section method to reconstruct the geological profile model, the spatial distribution data of the fault with higher accuracy can be obtained based on the fine geological shape geometric information in the geological profile model, thereby improving the accuracy of the spatial distribution data of the fault. Then, the permeation data determination module 403 constructs an earth dynamics model according to the geological profile model, takes the spatial distribution data of the fault as a boundary condition, and determines the permeation data of free water from the upper part of the fault downward based on the earth dynamics model; wherein the permeation data includes a permeation path, a permeation speed and a fault flux distribution; by limiting and constraining the permeation range of free water, the permeation of free water can be better simulated based on the spatial distribution data of the fault to obtain the permeation data result of free water from the upper part of the fault downward in the target area, thereby enhancing the significance and effectiveness of the permeation data in the target area. Finally, the serpentinization water influx calculation module 404 determines the serpentinization water influx of the serpentinization reaction zone according to the temperature and pressure field simulation result, the mineral stability data of the target area and the permeation data, and determines the target water influx in the target area according to the serpentinization water influx of the serpentinization reaction zone, wherein the target water influx includes the instantaneous water influx and / or cumulative water influx of the serpentinization water influx of the serpentinization reaction zone in the target area. In this way, by combining the geometric constraints and water permeation physical mechanisms of the target area, the serpentinization water influx of the serpentinization reaction zone with high precision is obtained, and then the instantaneous water influx and / or cumulative water influx of the target area is calculated, so that the change of the water influx of the target area is more comprehensively evaluated, and the expansibility and comprehensiveness of the technical scheme are improved. Thus, through the above technical steps, the change graph of the water influx of the target area with higher accuracy and the analysis and evaluation result are obtained, thereby effectively improving the utilization efficiency of the geophysical data and the accuracy of the water influx calculation of the target area.

[0121] On the basis of each of the optional technical solutions above, the permeation data determination module 403 optionally further includes a yield stress distribution information determination unit and a permeation data determination unit. The yield stress distribution information determination unit is configured to take the fault spatial distribution data as a boundary condition, perform yield stress field inversion based on the geodynamic model, and obtain yield stress distribution information of the fault. The permeation data determination unit is configured to determine permeation data of free water from the upper part of the fault downward according to the fault spatial distribution data and the yield stress distribution information of the fault.

[0122] On the basis of each of the optional technical solutions above, the permeation data determination unit can further include an effective permeability distribution data determination unit, a path constraint condition determination unit, and a permeation data determination subunit. The effective permeability distribution data determination unit is configured to determine effective permeability distribution data of the fault along the depth direction according to the fault spatial distribution data and the yield stress distribution information of the fault, and introduce an anisotropic tensor to control the permeation direction. The path constraint condition determination unit is configured to construct a multi-dimensional seepage water inlet path model evolving over time based on Darcy's law and a mass conservation control equation of the fault, and construct a path constraint condition according to the fault spatial distribution data of the fault; wherein the fault spatial distribution data includes nucleation parameters, expansion parameters, penetration depth, and permeation channel morphology. The permeation data determination subunit is configured to simulate permeation data of free water from the upper part of the fault downward according to the initial water content, the maximum permeation depth, the permeation direction, stress coupling data of permeability, fluid driving pressure difference and lithology control data, the path constraint condition, and the multi-dimensional seepage water inlet path model.

[0123] On the basis of each of the optional technical solutions above, the yield stress distribution information determination unit can be specifically configured to use the Mohr-Coulomb yield criterion, set shear strength and internal friction angle, use a stress field output by the geodynamic model, and invert the yield stress state of the fault to obtain the yield stress distribution information of the fault.

[0124] On the basis of each of the optional technical solutions described above, optionally, the geological structure includes an outer arc region, and the device can further include an outer arc region permeation data determination module. The outer arc region permeation data determination module is configured to, after the permeation data of the free water that penetrates downward from the upper part of the fault is simulated according to the initial water content, the maximum permeation depth, the permeation direction, the stress coupling data of the permeation rate, the fluid driving pressure difference and lithology control data, the path constraint condition, and the multi-dimensional seepage water inflow path model, determine the permeation data of the cracks in the outer arc region based on the mass conservation control equation of the cracks according to Darcy's law, update the permeation data of the free water that penetrates downward from the upper part of the fault according to the permeation data of the cracks.

[0125] On the basis of each of the optional technical solutions described above, optionally, the serpentinization water inflow flux calculation module 404 can further include a serpentinization water inflow flux determination unit. The serpentinization water inflow flux determination unit is configured to determine a serpentine stability zone of the target region according to the temperature and pressure field simulation results and the mineral stability data of the target region, determine a serpentinization reaction zone according to a preset temperature and pressure window, the serpentine stability zone, and the permeation data, and determine a serpentinization water inflow flux of the serpentinization reaction zone according to a water and peridotite reaction kinetics model.

[0126] On the basis of each of the optional technical solutions described above, optionally, the serpentinization water inflow flux determination unit can further include a serpentinization water inflow flux determination subunit. The serpentinization water inflow flux determination subunit is configured to construct a serpentinization rate model according to the water and peridotite reaction kinetics model, and determine the serpentinization water inflow flux of the serpentinization reaction zone according to the rock density, reaction volume, unit volume water content change of the serpentinization reaction zone, and the serpentinization rate model.

[0127] On the basis of each of the optional technical solutions described above, optionally, the device can further include a serpentinization water inflow flux updating module. The serpentinization water inflow flux updating module is configured to, for an observation region containing seismic tomography in the target region, process seismic tomography data according to a Voigt-Reuss-Hill average method and a corresponding relationship between P-wave velocity and water content established through experiments, convert a measured P-wave velocity distribution into a water content spatial distribution map, and update the serpentinization water inflow flux of the observation region according to the water content spatial distribution map.

[0128] On the basis of each of the optional technical solutions described above, optionally, the geophysical data includes at least one of multi-beam sounding data, multi-channel seismic reflection data, wide-angle seismic refraction data, seismic tomography imaging data, and marine geoelectric measurement data.

[0129] The training device of the geological identification model provided in the embodiment of the present application can execute the method for determining water flux of serpentinized mantle peridotite provided in any embodiment of the present application. The training device has the corresponding function modules and beneficial effects of the execution method. Technical details not described in detail in the present embodiment can be referred to the method for determining water flux of serpentinized mantle peridotite in any embodiment of the present application.

[0130] Embodiment five

[0131] Figure 5 A structural schematic diagram of an electronic device 10 that can be used to implement embodiments of the present application is shown. The electronic device is intended to represent various forms of digital computers, such as laptops, desktops, tablets, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device can also represent various forms of mobile devices such as personal digital processors, cellular telephones, smart phones, wearable devices (e.g., headsets, glasses, watches, etc.), and other similar computing devices. The components shown here, their connections and relationships, and their functions, are meant to be examples only, and are not meant to limit implementations of the present application described and / or claimed in this document.

[0132] As shown in Figure 5 The electronic device 10 includes at least one processor 11, and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11, wherein the memory stores a computer program that can be executed by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0133] A plurality of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, a speaker, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network, such as the Internet, and / or various telecommunication networks.

[0134] The processor 11 can be various general and / or special purpose processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 performs various methods and processes described above, such as the method of determining water flux into serpentinized mantle peridotite.

[0135] In some embodiments, the method of determining water flux into serpentinized mantle peridotite can be implemented as a computer program tangibly embodied in a computer readable storage medium, such as the storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded onto the RAM 13 and executed by the processor 11, one or more steps of the method of determining water flux into serpentinized mantle peridotite described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to perform the method of determining water flux into serpentinized mantle peridotite by any other suitable means, such as by means of firmware.

[0136] Various implementations of the systems and techniques described above can be realized in digital electronic circuitry, integrated circuitry, a field programmable gate array (FPGA), an application specific integrated circuit (ASIC), a system on a chip (SOC), a programmable logic device (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include implementation in one or more computer programs that are executable and / or interpretable on a programmable system including at least one programmable processor, which can be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.

[0137] Computer programs used to implement the methods of the application can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program, when executed by the processor of the machine, implements the functions / acts specified in the flowcharts and / or block diagrams. The computer program can be executed entirely on a machine, partially on a machine, partially on a machine as a stand-alone software package, and partially on a remote machine or a server.

[0138] In the context of the present application, a computer-readable storage medium can be a tangible medium that can contain or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. A computer-readable storage medium can include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium can be a machine-readable signal medium. More specific examples of a machine-readable storage medium will include one or more lines of a program of instructions in a transitory signal, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0139] To provide for interaction with a user, the systems and techniques described here can be implemented on an electronic device having a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the electronic device. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form, including acoustic, speech, or tactile input.

[0140] The systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a user computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), blockchain network, and the Internet.

[0141] The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. A server can be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system, to solve the defects of large management difficulty and weak business scalability in traditional physical host and VPS service.

[0142] In particular, the processes described above with reference to the flow charts can be implemented as computer software programs in accordance with embodiments of the application. For example, embodiments of the application include a computer program product comprising a computer program carried on a non-transitory computer readable medium, the computer program comprising program code for executing the methods illustrated by the flow charts. In such embodiments, the computer program can be downloaded and installed from a network via the communication unit 19, or installed from the storage unit 18, or installed from the ROM 12. When the computer program is executed by the processor 11, the above-described functions defined in the methods of the embodiments of the application are performed.

[0143] It should be understood that the various forms of flow shown above can be re-ordered, added to, or deleted from without departing from the scope of the present application. For example, the steps recited in the present application can be performed in parallel, in series, or in a different order, without departing from the desired results of the technical solutions of the present application, and this is not limited herein.

[0144] The above detailed description does not constitute a limitation on the protection scope of the present application. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent replacements, and improvements made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method of determining the water flux of serpentinized mantle peridotite, characterized in that, The method comprises: acquiring geophysical data of a target area, constructing a geological profile model for presenting a geological structure of the target area based on the geophysical data, and embedding initial water content and rock porosity data into the geological profile model in a two-dimensional grid manner, wherein the geological structure comprises at least one of a trench, an outer arc, an ocean basin, an oceanic crust and an upper mantle; reconstructing geometric information of a fault in the geological profile model by using a class profile method to obtain spatial distribution data of the fault in the target area, wherein the geometric information comprises at least one of fault throw data, dip angle data, spread data and activity direction data; constructing an earth dynamics model according to the geological profile model, determining penetration data of free water penetrating downward from an upper part of the fault based on the earth dynamics model with the spatial distribution data of the fault as a boundary condition, wherein the penetration data comprises a penetration path, a penetration speed and a fault flux distribution; determining a serpentinization inflow flux of a serpentinization reaction zone according to a temperature and pressure field simulation result, mineral stability data of the target area and the penetration data, and determining a target inflow flux in the target area according to the serpentinization inflow flux of the serpentinization reaction zone, wherein the target inflow flux comprises a transient inflow flux and / or a cumulative inflow flux of the serpentinization inflow flux of the serpentinization reaction zone in the target area.

2. The method of claim 1, wherein, The method of determining the penetration data of the free water penetrating downward from the upper part of the fault based on the earth dynamics model with the spatial distribution data of the fault as the boundary condition comprises: performing a yield stress field inversion based on the earth dynamics model with the spatial distribution data of the fault as a boundary condition to obtain yield stress distribution information of the fault; determining the penetration data of the free water penetrating downward from the upper part of the fault according to the spatial distribution data of the fault and the yield stress distribution information of the fault.

3. The method of claim 2, wherein, The method of determining the penetration data of the free water penetrating downward from the upper part of the fault according to the spatial distribution data of the fault and the yield stress distribution information of the fault comprises: determining effective permeability distribution data of the fault along a depth direction according to the spatial distribution data of the fault and the yield stress distribution information of the fault, and introducing an anisotropic tensor to control a penetration direction; constructing a time-evolving multi-dimensional seepage inflow path model based on Darcy's law and a mass conservation control equation of the fault, and constructing a path constraint condition according to the spatial distribution data of the fault, wherein the spatial distribution data of the fault comprises nucleation parameters, expansion parameters, penetration depth and penetration channel morphology; determining the penetration data of the free water penetrating downward from the upper part of the fault according to the initial water content, the maximum penetration depth, the penetration direction, permeability and stress coupling data, fluid driving pressure difference and lithology control data, the path constraint condition and the multi-dimensional seepage inflow path model.

4. The method of claim 2, wherein, The yield stress distribution information at least includes a yield stress state; and the yield stress field inversion based on the geodynamic model and taking the fault spatial distribution data as a boundary condition to obtain the yield stress distribution information of the fault, includes: The yield stress distribution information of the fault is obtained by setting a shear strength and an internal friction angle, using the stress field output by the geodynamic model, and inverting the yield stress state of the fault according to the Mohr-Coulomb yield criterion.

5. The method of claim 3, wherein, The geological structure includes an outer arc region; after the permeation data of the free water permeating downward from the upper part of the fault is simulated according to the initial water content, the maximum permeation depth, the permeation direction, the stress coupling data, the fluid driving pressure difference and lithology control data, the path constraint condition and the multi-dimensional seepage water inflow path model, the geological structure further includes: In the outer arc region, the permeation data of the fracture is determined based on the mass conservation control equation of the fracture according to Darcy's law, and the permeation data of the free water permeating downward from the upper part of the fault is updated according to the permeation data of the fracture.

6. The method of claim 1, wherein, The serpentinization water inflow flux of the serpentinization reaction zone is determined according to the simulation results of the temperature and pressure field, the mineral stability data of the target region and the permeation data, including: The serpentine stable zone of the target region is determined according to the simulation results of the temperature and pressure field and the mineral stability data of the target region, and the serpentinization reaction zone is determined according to the preset temperature and pressure window, the serpentine stable zone and the permeation data, and the serpentinization water inflow flux of the serpentinization reaction zone is determined according to the water and peridotite reaction kinetics model.

7. The method of claim 6, wherein, The serpentinization water inflow flux of the serpentinization reaction zone is determined according to the water and peridotite reaction kinetics model, including: The serpentinization rate model is constructed according to the water and peridotite reaction kinetics model, and the serpentinization water inflow flux of the serpentinization reaction zone is determined according to the rock density, reaction volume, unit volume water content change of the serpentinization reaction zone and the serpentinization rate model.

8. The method of claim 1, wherein, After the serpentinization water inflow flux of the serpentinization reaction zone is determined according to the simulation results of the temperature and pressure field, the mineral stability data of the target region and the permeation data, the method further includes: For the observation region containing seismic tomography in the target region, the seismic tomography data is processed according to the Voigt-Reuss-Hill average method and the corresponding relationship between the P-wave velocity and the water content established by experiments, the measured P-wave velocity distribution is converted into a water content spatial distribution map, and the serpentinization water inflow flux of the observation region is updated according to the water content spatial distribution map.

9. The method of claim 1, wherein, The geophysical data includes at least one of multi-beam sounding data, multi-channel seismic reflection data, wide-angle seismic refraction data, seismic tomography data and marine geoelectric measurement data.

10. A computer program product comprising a computer program, characterized in that, The computer program, when executed by a processor, implements the method for determining the serpentinization water inflow flux of the serpentine mantle peridotite as claimed in any one of claims 1-9.

Images