A three-dimensional rock facies modeling method constrained by deep-sea channel configuration

Abstract

The present invention discloses a three-dimensional modeling method of rock phase constrained by deep-sea channel configuration, including obtaining geological background, three-dimensional seismic, well logging curve, and drilling coring data in the target area; establishing a fine sequence stratigraphic framework of the target area based on the core-well logging-seismic comprehensive calibration template; combining the sequence framework interface and the theory of deep-sea channel configuration mode, carrying out the characterization and description of the configuration interface of deep-sea channels of different levels in the target area; encoding the configuration elements on a single well; comprehensively considering the structural characteristics of the target area, taking the top and bottom interfaces and fault distribution of the channel in the target area as constraints, establishing a structural geological model of the target area; taking the boundaries of different configuration units as constraints, establishing a three-dimensional geological model of the deep-sea channel configuration; taking the well filling lithofacies type and proportion as constraints, establishing a three-dimensional geological model of the rock phase of the deep-sea channel in the target area. The present invention is a three-dimensional model construction method that can accurately and finely reflect the sand body morphology, spatial distribution and superposition style of the deep-sea channel reservoir.

Description

Technical Field

[0001] The invention relates to a rock phase three-dimensional modeling method constrained by deep-sea channel configuration, and belongs to the field of oil and gas field development geology. Background Art

[0002] As an important sediment transport channel connecting the onshore and deep-sea environments, deep-sea water channels are also one of the main places for deep-sea turbidite sandstone deposition. They contain rich and valuable oil and gas resources and are the focus of research in the field of deep-sea sedimentation. However, although channel-type reservoirs have excellent porosity and permeability conditions, actual drilling of oil reservoirs shows that their plane morphology, internal structure, and filling lithology are complex and changeable. The sedimentary characteristics of channel reservoirs within a range of hundreds of meters horizontally may change dramatically. In addition, the spatial superposition relationship of sand bodies inside channel reservoirs is diverse, and the internal directional connectivity is complex, which brings severe challenges to the efficient development of this type of oil reservoir.

[0003] The original intention of establishing a three-dimensional model of an oil reservoir is to clarify the morphological characteristics, spatial distribution patterns and mutual superposition of underground reservoir sand bodies, and to provide geological guidance for the later development of oil reservoirs and the potential of remaining oil. Looking at the existing reservoir modeling methods, the more common ones include: (1) Reservoir interpolation modeling methods, which use interpolation methods to use existing data label areas to complete data in other areas of the reservoir. There are single data label extension interpolation and multi-data label interpolation, and then establish a three-dimensional geological model of the target reservoir. It is a deterministic modeling method. According to the difference in its interpolation algorithm, it is divided into mathematical statistics interpolation method and Kriging interpolation method; (2) Reservoir random modeling method, based on known information and random functions as theory, it applies random simulation methods (such as sequential indicator simulation) to generate a reservoir model that conforms to actual knowledge. According to the difference in research objects, it is divided into two categories: target-based method and pixel-based method. However, the deep-sea environment has the characteristics of small number of wells, sparse well network, and complex heterogeneity of reservoirs between wells. There are few available data labels, and the accuracy of the algorithm model of the relationship between well labels is low. As a result, the application of existing geological modeling methods in deep-sea channel-type reservoirs has the problems of low interpolation accuracy, large deviation, and inconsistency between interpolation results and geological understanding, which seriously restricts the efficient development of this type of oil reservoir. Summary of the invention

[0004] In order to overcome the defects in the prior art, the present invention aims to provide a three-dimensional modeling method of rock phases constrained by deep-sea channel configuration, which can accurately reflect the sand body morphology, spatial distribution and superposition style of deep-sea channels of different configuration levels underground in three-dimensional space.

[0005] The technical solution provided by the present invention to solve the above technical problems is: a three-dimensional rock phase modeling method constrained by deep-sea channel configuration, comprising the following steps:

[0006] Step S10, obtaining geological background, 3D seismic, well logging curves, and drilling coring data of the target area;

[0007] Step S20, establishing a fine sequence stratigraphic framework for the target area based on the core-logging-seismic integrated calibration template;

[0008] Step S30, combining the sequence framework interface and deep-sea channel configuration model theory, characterizing and describing the deep-sea channel configuration interfaces of different levels in the target area;

[0009] Step S40, encoding configuration elements on a single well based on the characterization result of the deep-sea channel configuration interface;

[0010] Step S50, based on the deep-sea channel configuration interface representation and configuration element coding results, the structural characteristics of the target area are integrated, and the top and bottom interfaces of the channel and the fault distribution in the target area are used as constraints to establish a structural geological model of the target area;

[0011] Step S60: On the basis of the structural geological model, a three-dimensional geological model of the deep-sea channel configuration is established with the boundaries of different configuration units as constraints;

[0012] Step S70: On the basis of the three-dimensional geological model of the deep-sea channel configuration, a three-dimensional geological model of the rock phase of the deep-sea channel in the target area is established with the type and proportion of the well filling rock phase as constraints.

[0013] A further technical solution is that the specific steps of step S20 are:

[0014] Step S21, based on the large-scale sandstone-mudstone boundary and cycle changes on the logging curve, establish the third-level sequence interface interpretation framework on the well;

[0015] Step S22: Based on the interpretation framework of the third-order sequence interface on the well, the position of the third-order sequence interface of the seismic profile is determined by using the well-seismic joint analysis, the sequence framework interpretation of the entire target area is carried out, and a fine sequence stratigraphic framework of the target area is established.

[0016] A further technical solution is that the specific steps of step S30 are:

[0017] Step S31, under the constraint of the third-order sequence interface, based on the guidance of the channel system sedimentary architecture model, a template for well-seismic identification of the channel system is established, and then the identification and characterization of the hierarchical architecture interface of the channel system is carried out in the whole area;

[0018] Step S32: Under the constraints of the channel system configuration interface, based on the guidance of the composite channel sedimentary configuration model and in combination with the logging and seismic response characteristics, the composite channel hierarchical configuration edge of the target area is identified by using the vertical staging and lateral delimitation methods;

[0019] Step S33: Under the constraint of the composite channel configuration interface, based on the guidance of the single channel sedimentary configuration model and combined with the single channel well-seismic calibration results, the single channel hierarchical configuration boundary identification and characterization are carried out in the study area.

[0020] A further technical solution is that the specific process of step S40 is: based on the characterization result of the deep-sea channel configuration interface, corresponding interpretation elements are projected on a single well and corresponding labels are established, and then the configuration labels on all wells in the target area are discretized.

[0021] A further technical solution is that the specific process of step S50 is: based on the deep-sea channel configuration interface characterization and configuration element coding results, a third-level sequence top and bottom interface isoT0 map is established, and a structural geological model of the target area is formed using the isoT0 map and faults as a skeleton.

[0022] A further technical solution is that the specific process of step S60 is:

[0023] Step S61: Based on the structural geological model of the target area, combined with the well seismic characterization and plane attribute extraction results, the boundary of the waterway system is characterized and the corresponding plane boundary polygon is established. Then, through the deterministic assignment method, the waterway system code CS is input under the constraint of the waterway system plane polygon. Then, under the constraint of the profile waterway system interface, the boundary morphology of the waterway system is adjusted manually to obtain the final three-dimensional geological model of the waterway system that conforms to the geological understanding and geological model.

[0024] Step S62: Based on the three-dimensional geological model of the water channel system, the plane position projection point of the lateral boundary of the composite water channel is determined with the well-seismic configuration characterization result as a hard constraint, and the plane map of seismic attributes extracted from the sublayer where the composite water channel is located is used as a soft constraint to characterize the plane boundary of the composite water channel and establish the corresponding plane boundary polygon. Through a deterministic assignment method, the composite water channel code CC is input under the constraint of the composite water channel plane polygon; then, under the constraint of the composite water channel profile interface, the composite water channel boundary morphology is adjusted by manual intervention, the cut-over surface is depicted in the cut-over area of ​​the composite water channel, and the composite water channel boundary code CC-Boundary is input, so as to obtain a final three-dimensional geological model of the composite water channel that conforms to the geological understanding and geological model;

[0025] Step S63, based on the three-dimensional geological model of the composite waterway, taking the well-seismic configuration characterization result as a hard constraint, determining the plane position projection point of the lateral boundary of the single waterway, and taking the plane map of seismic attributes extracted from the sublayer where the single waterway is located as a soft constraint, characterizing the plane boundary of the single waterway and establishing the corresponding plane boundary polygon, and inputting the single waterway code C under the constraint of the single waterway plane polygon through a deterministic assignment method; then, under the constraint of the single waterway profile boundary, manually adjusting the single waterway boundary morphology, characterizing the cut-over surface in the single waterway cut-over area, and inputting the single waterway boundary code C-Boundary, and finally obtaining the single waterway three-dimensional geological model that conforms to the geological understanding and geological model; under the constraint of the structural model of the target area, and based on the structural geological model of the target area, deterministic assignment is performed in layers, so as to construct a three-dimensional geological model of the waterway system;

[0026] A further technical solution is that the specific process of step S70 is:

[0027] Step S71: Based on the type and proportion of lithofacies filling inside the waterway analyzed in the target area, corresponding lithofacies labels are established and discretized;

[0028] Step S72: Under the constraints of a single water channel three-dimensional geological model, starting from the lithofacies types and proportions identified on the well, set the distribution ratio of the main body and edge of the water channel, and then use the truncated Gaussian simulation method to simulate the three-dimensional distribution of turbidite channel lithofacies, and establish a three-dimensional geological model of the deep-sea channel lithofacies in the target area.

[0029] The present invention has the following beneficial effects: The present invention is a three-dimensional model construction method that can accurately and finely reflect the sand body morphology, spatial distribution and superposition style of deep-sea channel-type reservoirs. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 It is a flowchart of the present invention;

[0031] Figure 2 This is the result of the joint characterization of well and seismic data of the third-level sequence interface in the target area;

[0032] Figure 3 This is the result diagram of the characterization of the interface configuration of the waterway system in the target area;

[0033] Figure 4 This is the result diagram of characterization of composite water channel interface configuration in the target area;

[0034] Figure 5 This is the result diagram of the characterization of the interface configuration of a single water channel in the target area;

[0035] Figure 6 This is the single-well interpretation result diagram of the sedimentary architecture elements of the GEN-P03 well;

[0036] Figure 7 It is the geological model map of the deep-sea channel reservoir structure in the target area;

[0037] Figure 8 It is a three-dimensional geological model of the waterway system level in the target area (I / J / K directions);

[0038] Fig. 9 It is a three-dimensional geological model map of the composite waterway level in the target area;

[0039] Fig.10 A three-dimensional geological model diagram of a single waterway configuration in the target area;

[0040] Fig.11 This is the structural elements of the turbidite channel on the well and the results of the lithofacies interpretation;

[0041] Fig.12 This is the result of deep-sea channel lithofacies simulation. DETAILED DESCRIPTION

[0042] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0043] like Figure 1 As shown, a rock phase three-dimensional modeling method constrained by deep-sea channel configuration of the present invention comprises the following steps:

[0044] Step S10: obtaining geological background, 3D seismic, well logging curves, drilling coring and other data of the target area;

[0045] Step S20: Based on the core-logging-seismic integrated calibration template, a fine sequence stratigraphic framework of the target area is established to complete the identification and characterization of the third-level sequence interface;

[0046] Step S21, based on the large-scale sandstone-mudstone boundary and cycle changes on the logging curve, establish the third-level sequence interface interpretation framework on the well;

[0047] Specifically, if Figure 2 As shown in the figure, based on the large-scale cyclic variation characteristics of the high-low gamma logging curve on a single well, the location of the third-order sequence interface is identified and the third-order sequence anchor point on the well is established;

[0048] Step S22: Based on the identification results of the third-order sequence interface on the well, the well-seismic joint method is used to identify the position and characteristics of the event axis corresponding to the third-order sequence interface on the seismic section, and based on this, the seismic interpretation of the third-order sequence interface is carried out in the whole area to establish a fine sequence stratigraphic framework in the target area;

[0049] Step S30, taking the third-order sequence interface as a constraint and under the guidance of the deep-sea channel sedimentary architecture model, characterizing and depicting the deep-sea channel architecture interfaces of different orders in the target area;

[0050] Step S31, under the constraint of the third-order sequence interface, based on the guidance of the channel system sedimentary architecture model, a template for well-seismic identification of the channel system is established, and then the identification and characterization of the hierarchical architecture interface of the channel system is carried out in the whole area;

[0051] Specifically, within the third-order sequence interface, the regional distribution of low gamma in the well is used, combined with the distribution boundary of chaotic strong amplitude and weak amplitude on the seismic section, to identify and characterize the boundary of the waterway system architecture, such as Figure 3 As shown, channel system II shows a large "U" shape in cross section with a relatively flat bottom, and channel system III shows an approximate "V" shape in cross section with a relatively protruding bottom;

[0052] Step S32: Under the constraints of the channel system configuration interface, based on the guidance of the composite channel sedimentary configuration model and in combination with the logging and seismic response characteristics, the composite channel hierarchical configuration boundary of the target area is identified by using the vertical staging and lateral delimitation methods;

[0053] Under the constraints of the channel system configuration interface, the composite waterway is divided into stages vertically according to the regular changes of the seismic reflection axis and the cyclic change characteristics of the well logging curve, and then the composite waterway interface of each stage is tracked on the seismic section; on the basis of the vertical stage division, the root mean square attribute plane map within the layer is extracted, and the configuration boundary is characterized by the plane-section interaction method to complete the composite waterway configuration characterization. Figure 4 As shown in the figure, the composite water channels in the target area can be divided into ten stages (CC1~CC10). Vertically, the composite water channels in different composite water channel series have a certain superposition pattern. In the early stage of deposition, the composite water channels are mainly contacted and spliced ​​in the vertical direction, while in the late stage of deposition, the composite water channels are mainly cut and overlapped in the lateral direction.

[0054] Step S33, under the constraint of the composite channel configuration interface, based on the guidance of the single channel sedimentary configuration model and combined with the single channel well-seismic calibration results, the single channel hierarchical configuration boundary of the study area is identified and characterized;

[0055] Under the constraint of the composite channel configuration boundary, the single channel configuration boundary is mainly tracked in the lateral direction according to the change of the logging curve and the sudden change interface of the seismic reflection axis. On the seismic profile, due to the change of sand and mud content at the boundary of the single channel, the seismic amplitude suddenly changes, that is, the strong / weak amplitude alternates, and an obvious sudden change point will appear. Figure 5 As shown, different numbers of single channels are developed in the composite channel of the target area. In the early stage of deposition, the single channels are mainly incised, overlapped and spliced ​​laterally, while in the late stage of deposition, the single channels begin to be incised, overlapped and spliced ​​laterally.

[0056] Step S40: Based on the characterization result of the deep-sea channel configuration interface, the configuration element recognition and coding are carried out on the single well, and the configuration elements interpreted and coded on the single well are discretized so that the coding results are evenly distributed;

[0057] Specifically, if Figure 6 As shown in the figure, the interpreted configuration interfaces of different levels are projected onto a single well, and a corresponding sedimentary configuration code is newly created on the single well. CS (Channel System), CC (Complex Channel) and C (Channel) correspond to the channel system, complex channel and single channel respectively, and discretization is performed on this basis;

[0058] Step S50, based on the deep-sea channel configuration interface representation and configuration element coding results, the structural characteristics of the target area are integrated, and the top and bottom interfaces of the channel and the fault distribution in the target area are used as constraints to establish a structural geological model of the target area;

[0059] Specifically, based on the three-level sequence interface interpreted in the early stage, the structural model function was used to establish the overall structural model of the study area with the top and bottom interfaces as the framework, and the lithology label of deep-sea mudstone was given as the background lithofacies of the target area.

[0060] Step S60: Based on the characterization results of each level of the deep-sea channel and the discretization results of the well configuration elements, different designated color codes are assigned to distinguish the configuration units with the boundaries of different configuration units as constraints, and special labels are assigned to the configuration boundaries to achieve fine characterization of the sand body morphology, spatial distribution and superposition style of different deep-sea channel reservoirs, thereby establishing a three-dimensional geological model of the deep-sea channel configuration;

[0061] Step S61: Based on the structural geological model of the target area, combined with the well seismic characterization and plane attribute extraction results, the boundary of the waterway system is characterized and the corresponding plane boundary polygon is established. Then, through the deterministic assignment method, the waterway system code CS is input under the constraint of the waterway system plane polygon. Then, under the constraint of the profile waterway system interface, the boundary morphology of the waterway system is adjusted manually to obtain the final three-dimensional geological model of the waterway system that conforms to the geological understanding and geological model.

[0062] like Figure 8 As shown in the figure, based on the structural model of the target area, combined with the well-seismic joint water channel system interface characterization results and seismic attribute map, the maximum boundary of the water channel system plane is clarified and a polygon is established. The CS (water channel system) data label is assigned within the polygon range through a deterministic assignment method. Then, in the three directions of I, J, and K, the edge morphology is depicted at the edge of the water channel system with reference to the water channel interface, so that the model is consistent with the actual geological understanding.

[0063] Step S62: Based on the three-dimensional geological model of the water channel system, the well-seismic configuration characterization results are used as hard constraints to determine the plane position projection points of the lateral boundary of the composite water channel. At the same time, the seismic attribute plane map extracted from the sublayer where the composite water channel is located is used as a soft constraint to characterize the plane boundary of the composite water channel and establish the corresponding plane boundary polygon. Through a deterministic assignment method, the composite water channel code CC is input under the constraint of the composite water channel plane polygon; then, under the constraint of the composite water channel profile interface, the composite water channel boundary morphology is manually intervened to characterize the cut-over surface in the composite water channel cut-over area, and the composite water channel boundary code CC-Boundary is input to obtain the final composite water channel three-dimensional geological model that conforms to geological knowledge and geological model;

[0064] like Fig. 9 As shown in the figure, based on the three-dimensional geological model of the waterway system, combined with the characterization results of the composite waterway interface combined with well-seismic data and the seismic attribute map, the maximum boundary of the composite waterway plane is clarified and a polygon is established. The CC (composite waterway) data label is assigned within the polygon range through a deterministic assignment method, and then the edge morphology is depicted at the edge of the waterway system in the three directions of I, J, and K with reference to the waterway interface; at the same time, in the area where some composite waterways overlap each other, the length and shape of the composite waterway overlap surface are depicted with reference to the composite waterway configuration overlap style to show the overlap direction, and the CC-Boundary data label is input to distinguish and separate composite waterways of different periods, so that the model is consistent with the actual geological understanding;

[0065] Step S63, based on the three-dimensional geological model of the composite waterway, taking the well-seismic configuration characterization result as a hard constraint, determining the plane position projection point of the lateral boundary of the single waterway, and taking the plane map of seismic attributes extracted from the sublayer where the single waterway is located as a soft constraint, characterizing the plane boundary of the single waterway and establishing the corresponding plane boundary polygon, and inputting the single waterway code C under the constraint of the single waterway plane polygon through a deterministic assignment method; then, under the constraint of the single waterway profile boundary, manually adjusting the single waterway boundary morphology, characterizing the cut-over surface in the single waterway cut-over area, and inputting the single waterway boundary code C-Boundary, and finally obtaining the single waterway three-dimensional geological model that conforms to the geological understanding and geological model; under the constraint of the structural model of the target area, and based on the structural geological model of the target area, deterministic assignment is performed in layers, so as to construct a three-dimensional geological model of the waterway system;

[0066] like Fig.10As shown in the figure, based on the three-dimensional geological model of the composite waterway, combined with the characterization results of the single waterway interface and the seismic attribute map, the maximum boundary of the single waterway plane is clarified and a polygon is established. The C (single waterway) data label is assigned within the polygon range through a deterministic assignment method, and then the edge morphology is depicted at the edge of the single waterway in the three directions of I, J, and K with reference to the single waterway interface; at the same time, in the area where some single waterways overlap each other, the length and shape of the overlapping surface between the single waterways are depicted with reference to the overlapping style of the single waterway configuration to show the overlapping direction, and the C-Boundary data label is input to distinguish and separate single waterways of different periods, so that the model is consistent with the actual geological understanding;

[0067] Step S70, based on the three-dimensional geological model of the deep-sea channel configuration, the lithofacies distribution characteristics in the waterway are simulated with the lithofacies filling type and proportion inside the waterway identified on the well as constraints to obtain a three-dimensional geological model of the lithofacies of the deep-sea channel;

[0068] Step S71: Based on the type and proportion of lithofacies filling inside the waterway analyzed in the target area, corresponding lithofacies labels are established and discretized;

[0069] Step S72: Under the constraints of a single water channel three-dimensional geological model, starting from the lithofacies types and proportions identified on the well, the distribution ratios of the main body and edge of the water channel are set, and then the truncated Gaussian simulation method is used to simulate the three-dimensional distribution of turbidite channel lithofacies.

[0070] Specifically, under the constraints of a single waterway geological model, combined with logging curves and coring data, the lithofacies genesis type on a single well and its proportion in the waterway, such as debris flow, landslide, Bouma sequence, block transport deposits and abandoned mudstone, etc., are identified, and corresponding labels are established on the well, such as Fig.11 Then, within the scope of a single waterway, based on the lithofacies labels on the wells in the target area, starting from the well area, the sequential Gaussian truncation algorithm is used to simulate the three-dimensional distribution of lithofacies inside a single waterway, as shown in Fig.12 As shown in the figure, from bottom to top of the channel in the direction perpendicular to the provenance, there is a positive cyclic sequence of conglomerate deposition-medium-coarse-grained Bouma sequence Ta / Tb-medium-fine-grained Tc segment-mud debris or fine-grained mud deposition, and the superposition relationship (oblique arrangement, vertical superposition, etc.) between single channels is obvious; from the perspective of parallel channel deposition, the positive grain sequence characteristics inside the channel are obvious, the continuity of the sand body is enhanced, and affected by the degree of channel curvature, a complete single channel morphology can still be seen at the cross-sectional position.

[0071] The above description is not intended to limit the present invention in any form. Although the present invention has been disclosed through the above embodiments, it is not intended to limit the present invention. Any technician familiar with the profession can make some changes or modifications to equivalent embodiments of equivalent changes using the technical contents disclosed above without departing from the scope of the technical solution of the present invention. However, any simple modification, equivalent changes and modifications made to the above embodiments based on the technical essence of the present invention without departing from the content of the technical solution of the present invention are within the scope of the technical solution of the present invention.

Claims

1. A three-dimensional modeling method of rock phases constrained by deep-sea channel configuration, characterized in that: The following steps are involved: Step S10, obtaining geological background, 3D seismic, well logging curves, and drilling coring data of the target area; Step S20, establishing a fine sequence stratigraphic framework for the target area based on the core-logging-seismic integrated calibration template; Step S30, combining the sequence framework interface and deep-sea channel configuration model theory, characterizing and describing the deep-sea channel configuration interfaces of different levels in the target area; Step S40, encoding configuration elements on a single well based on the characterization result of the deep-sea channel configuration interface; Step S50, based on the deep-sea channel configuration interface representation and configuration element coding results, the structural characteristics of the target area are integrated, and the top and bottom interfaces of the channel and the fault distribution in the target area are used as constraints to establish a structural geological model of the target area; Based on the characterization of deep-sea channel configuration interface and the coding results of configuration elements, the isoT0 map of the top and bottom interfaces of the third-level sequence is established, and the structural geological model of the target area is formed with the isoT0 map and faults as the framework; Step S60: On the basis of the structural geological model, a three-dimensional geological model of the deep-sea channel configuration is established with the boundaries of different configuration units as constraints; Step S61: Based on the structural geological model of the target area, combined with the well seismic characterization and plane attribute extraction results, the boundary of the waterway system is characterized and the corresponding plane boundary polygon is established. Then, through the deterministic assignment method, the waterway system code CS is input under the constraint of the waterway system plane polygon. Then, under the constraint of the profile waterway system interface, the boundary morphology of the waterway system is adjusted manually to obtain the final three-dimensional geological model of the waterway system that conforms to the geological understanding and geological model. Step S62: Based on the three-dimensional geological model of the water channel system, the plane position projection point of the lateral boundary of the composite water channel is determined with the well-seismic configuration characterization result as a hard constraint, and the plane map of seismic attributes extracted from the sublayer where the composite water channel is located is used as a soft constraint to characterize the plane boundary of the composite water channel and establish the corresponding plane boundary polygon. Through a deterministic assignment method, the composite water channel code CC is input under the constraint of the composite water channel plane polygon; then, under the constraint of the composite water channel profile interface, the composite water channel boundary morphology is adjusted by manual intervention, the cut-over surface is depicted in the cut-over area of ​​the composite water channel, and the composite water channel boundary code CC-Boundary is input, so as to obtain a final three-dimensional geological model of the composite water channel that conforms to the geological understanding and geological model; Step S63, based on the three-dimensional geological model of the composite waterway, the plane position projection point of the lateral boundary of the single waterway is determined with the result of the well-seismic configuration characterization as a hard constraint, and the plane map of the seismic attributes extracted from the sublayer where the single waterway is located is used as a soft constraint to characterize the plane boundary of the single waterway and establish the corresponding plane boundary polygon, and input the single waterway code C under the constraint of the single waterway plane polygon through a deterministic assignment method; then, under the constraint of the single waterway profile boundary, the single waterway boundary morphology is manually adjusted, the cut-over surface is depicted in the cut-over area of ​​the single waterway, and the single waterway boundary code C-Boundary is input, so as to obtain the final three-dimensional geological model of the single waterway that conforms to the geological understanding and geological model; Step S70: On the basis of the three-dimensional geological model of the deep-sea channel configuration, a three-dimensional geological model of the rock phase of the deep-sea channel in the target area is established with the type and proportion of the well filling rock phase as constraints; Step S71: Based on the type and proportion of lithofacies filling inside the waterway analyzed in the target area, corresponding lithofacies labels are established and discretized; Step S72: Under the constraints of a single water channel three-dimensional geological model, starting from the lithofacies types and proportions identified on the well, set the distribution ratios of the main body and edge of the water channel, use the truncated Gaussian simulation method to simulate the three-dimensional distribution of turbidite channel lithofacies, and establish a three-dimensional geological model of the deep-sea channel lithofacies in the target area.

2. The method for three-dimensional rock facies modeling constrained by deep-sea channel configuration according to claim 1, characterized in that: The specific steps of step S20 are: Step S21, based on the large-scale sandstone-mudstone boundary and cycle changes on the logging curve, establish the third-level sequence interface interpretation framework on the well; Step S22: Based on the interpretation framework of the third-order sequence interface on the well, the position of the third-order sequence interface of the seismic profile is determined by using the well-seismic joint analysis, the sequence framework interpretation of the entire target area is carried out, and a fine sequence stratigraphic framework of the target area is established.

3. The method for three-dimensional rock facies modeling constrained by deep-sea channel configuration according to claim 1, characterized in that: The specific steps of step S30 are: Step S31, under the constraint of the third-order sequence interface, based on the guidance of the channel system sedimentary architecture model, a template for well-seismic identification of the channel system is established, and then the identification and characterization of the hierarchical architecture interface of the channel system is carried out in the whole area; Step S32: Under the constraints of the channel system configuration interface, based on the guidance of the composite channel sedimentary configuration model and in combination with the logging and seismic response characteristics, the composite channel hierarchical configuration boundary of the target area is identified by using the vertical staging and lateral delimitation methods; Step S33: Under the constraint of the composite channel configuration interface, based on the guidance of the single channel sedimentary configuration model and combined with the single channel well-seismic calibration results, the single channel hierarchical configuration boundary identification and characterization are carried out in the study area.

4. The method for three-dimensional rock facies modeling constrained by deep-sea channel configuration according to claim 1, characterized in that: The specific process of step S40 is: based on the characterization result of the deep-sea channel configuration interface, project the corresponding interpretation elements on the single well and establish the corresponding labels, and then discretize the configuration labels on all wells in the target area.

Images