WO2022228450A1 - 基于非三角洲海侵的大陆架剖面确定方法及装置 - Google Patents
基于非三角洲海侵的大陆架剖面确定方法及装置 Download PDFInfo
- Publication number
- WO2022228450A1 WO2022228450A1 PCT/CN2022/089455 CN2022089455W WO2022228450A1 WO 2022228450 A1 WO2022228450 A1 WO 2022228450A1 CN 2022089455 W CN2022089455 W CN 2022089455W WO 2022228450 A1 WO2022228450 A1 WO 2022228450A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- delta
- transgression
- continental shelf
- river
- profile
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000013049 sediment Substances 0.000 claims abstract description 72
- 230000015654 memory Effects 0.000 claims description 36
- 238000003860 storage Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000011144 upstream manufacturing Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 description 24
- 230000006870 function Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 14
- 238000004891 communication Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 239000000872 buffer Substances 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 3
- 238000005260 corrosion Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000002085 persistent effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Definitions
- the invention relates to the technical field of oil and natural gas exploration and marine geology, in particular to the technical field of stratum superposition pattern analysis, and in particular to a method and device for determining a continental shelf profile based on non-delta transgression.
- the calculation of the continental shelf profile is based on the erosion effect of waves on the seabed.
- the basic principle is that the water depth of the continental shelf gradually becomes shallower from the direction of the deep sea to the direction of the continent. The stronger the effect, the stronger the damage to the sediment. Sediments that have been destroyed in shallow water will be transported to relatively deep water and deposited. This process will eventually reach an equilibrium state, in which case the profile of the continental shelf is also called the equilibrium profile. Therefore, the equilibrium profile of the continental shelf should have the characteristic of steepening towards the land.
- Bruun extended this feature to the background of sea level rise, called Bruun's law (Bruun, 1962).
- the existing technology is applicable to a small time scale, in the range of tens to hundreds of years, belonging to the field of engineering, and it is assumed that external factors such as geological structure and sediment supply do not change. For the time scale of thousands to hundreds of thousands of years, the external conditions are also slowly changing, and the above techniques are difficult to apply.
- the method and device for determining the continental shelf profile based on non-delta transgression provided by the present invention, under the background of non-delta transgression, comprehensively consider the influence of river sediment supply and sea level rise on the continental shelf profile, In this way, the changes of the continental shelf profile in the historical period of the target project area can be more accurately determined.
- the present invention provides a method for determining a continental shelf profile based on non-delta transgression, comprising:
- a two-dimensional model of the continental shelf profile was established according to the amount of sediment supplied by the river in the target project area and the rate of sea level rise;
- the establishment of a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target project area and the rate of sea level rise includes:
- the present invention provides a device for determining a continental shelf profile based on non-delta transgression, the device comprising:
- the model building unit is used to build a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target project area and the rate of sea level rise;
- a model solving unit used for solving the two-dimensional model of the continental profile in the space rectangular coordinate system, so as to generate a relational expression that can characterize the non-delta transgressive slope and the length of the river;
- the model establishment unit includes:
- the river At the onset of sea level rise, the river has an initial length.
- the model solving unit is specifically configured to solve the two-dimensional model of the continental section in the x-z coordinate system to generate the relational expression.
- the present invention provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and running on the processor, when the processor executes the program, a method for determining a continental shelf profile based on non-delta transgression is implemented. step.
- the present invention provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, implements the steps of a method for determining a continental shelf profile based on non-delta transgression.
- the method and device for determining the continental shelf profile based on non-delta transgression provided by the embodiments of the present invention firstly establish a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target work area and the sea level rise rate; Under the Cartesian coordinate system, the two-dimensional model of the continental profile is solved to generate a relational expression that can characterize the non-delta transgressive slope and river length; finally, the continental shelf profile of the target project area is determined according to the relational expression.
- the invention proposes that under the background of non-delta transgression, comprehensively considering the influence of river sediment supply and sea level rise on the continental shelf profile, the variation of the continental shelf profile in the historical period of the target project area can be more accurately determined.
- FIG. 1 is a schematic diagram of a continental shelf balance cross-section based on Brunn's law in the background of the present invention
- FIG. 2 is a schematic flowchart of a method for determining a continental shelf profile based on non-delta transgression in an embodiment of the present invention
- step 100 is a schematic flowchart of step 100 in an embodiment of the present invention.
- FIG. 4 is a schematic diagram of the longitudinal area increment of the alluvial range formed in the time increment ⁇ t in the embodiment of the present invention.
- FIG. 5 is a schematic flowchart of step 200 in an embodiment of the present invention.
- FIG. 7 is a diagram showing the relationship between the continental shelf gradient ( ⁇ ) and the dimensionless river length (L * ) in a specific application example of the present invention.
- Fig. 8 is a schematic diagram of setting the x-z coordinate system in a specific application example of the present invention.
- FIG. 9 is a schematic flowchart of a method for determining a continental shelf profile based on non-delta transgression in a specific application example of the present invention.
- FIG. 10 is a schematic structural diagram of a device for determining a continental shelf profile based on non-delta transgression in an embodiment of the present invention
- FIG. 11 is a schematic diagram of the composition and structure of the model building unit 10 in the device for determining the continental shelf profile based on non-delta transgression according to an embodiment of the present invention
- FIG. 12 is a schematic block diagram of a system configuration of an electronic device 600 according to an embodiment of the present invention.
- embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
- computer-usable storage media including, but not limited to, disk storage, CD-ROM, optical storage, etc.
- An embodiment of the present invention provides a specific implementation of a method for determining a continental shelf profile based on non-delta transgression. Referring to FIG. 2, the method specifically includes the following contents:
- the target work area is in the sedimentary background of non-delta transgression.
- transgression when the sea level retreats towards the land, it is called transgression, and new continental shelves can be generated.
- Transgression can be further divided into two cases: delta transgression and non-delta transgression.
- delta transgression The end of the river with the development of delta
- non-delta transgression which depends on the length of the river.
- L crt critical length
- L ⁇ L crt it is a non-delta transgression.
- Step 200 Under the space rectangular coordinate system, solve the two-dimensional model of the continental profile to generate a relational expression that can characterize the non-delta transgressive slope and the length of the river.
- the continental shelf is placed in the x-z coordinate system (x represents the horizontal direction, z represents the vertical direction), and the slope of the continental shelf profile at a certain position is equivalent to the tangent slope of the continental shelf profile at that position. Therefore, in this x-z coordinate system, the solution of the non-delta transgressive profile can be found.
- Step 300 Determine the continental shelf profile of the target work area according to the relational expression.
- step 300 can represent the non-delta transgressive slope, it can determine the continental shelf profile of the target work area.
- the method for determining the continental shelf profile based on non-delta transgression provided by the embodiment of the present invention firstly establishes a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target work area and the rate of sea level rise; Under the system, the two-dimensional model of the continental profile is solved to generate a relational expression that can characterize the non-delta transgressive slope and the length of the river; finally, the continental shelf profile of the target project area is determined according to the relational expression.
- the invention proposes that under the background of non-delta transgression, comprehensively considering the influence of river sediment supply and sea level rise on the continental shelf profile, the variation of the continental shelf profile in the historical period of the target project area can be more accurately determined.
- the river At the onset of sea level rise, the river has an initial length.
- Equation (9) is the general form of the slope and river length (dimensionless) of any non-delta transgressive profile.
- L has the form:
- Equation (12) has no analytical solution, but a numerical solution can be obtained by numerical simulation.
- Formulas (11)–(13) can also divide the numerator and denominator on the right side of the equal sign by ⁇ 2D to obtain the general expression form of dimensionless transgression profile.
- the method for determining the continental shelf profile based on non-delta transgression provided by the embodiment of the present invention first selects independent variables based on the internal detection data of the pipeline, uses the Lasso algorithm to select the input variables, and finally determines the internal corrosion defects of the pipeline.
- the prediction adopts a generalized linear additive model driven by internal detection data, and the embodiment of the present invention can perform modeling based on the existing two internal corrosion data to estimate the future corrosion depth. Since the analysis results of the model can provide a reference for the dangerous defects of the pipeline, it can be an important part of the integrity management, which can help determine the internal inspection cycle and formulate maintenance plans, which is beneficial to the safe operation of the pipeline.
- the beneficial effects of the embodiments of the present invention are as follows:
- a non-delta transgression-based continental shelf profile determination device capable of realizing a non-delta transgression-based continental shelf profile determination method.
- a non-delta transgression-based continental shelf profile determination device Specifically include the following:
- the model establishment unit 10 is used for establishing a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target project area and the rate of sea level rise;
- the model establishment unit 10 includes:
- the longitudinal increment calculation module 101 is used to calculate the longitudinal increment of the alluvial range in the target work area per unit time
- the river supply volume calculation module 103 is configured to calculate the river supply sediment volume according to the sediment volume per unit time.
- the river At the onset of sea level rise, the river has an initial length.
- the model solving unit 20 is specifically configured to solve the two-dimensional model of the continental section in the x-z coordinate system to generate the relational expression.
- the device for determining the continental shelf profile based on non-delta transgression firstly establishes a two-dimensional model of the continental shelf profile according to the amount of sediment supplied by the river in the target work area and the sea level rise rate; Under the system, the two-dimensional model of the continental profile is solved to generate a relational expression that can characterize the non-delta transgressive slope and the length of the river; finally, the continental shelf profile of the target project area is determined according to the relational expression.
- the present invention proposes that under the background of non-delta transgression, comprehensively considering the influence of river sediment supply and sea level rise on the continental shelf profile, the change of the continental shelf profile in the historical period of the target project area can be more accurately determined.
- FIG. 12 is a schematic block diagram of a system configuration of an electronic device 600 according to an embodiment of the present invention.
- the electronic device 600 may include a central processing unit 100 and a memory 140 ; the memory 140 is coupled to the central processing unit 100 .
- this figure is exemplary; other types of structures may be used in addition to or in place of this structure to implement telecommunication functions or other functions.
- the non-delta transgression-based continental shelf profile determination method may be integrated into the central processing unit 100 .
- the central processing unit 100 can be configured to perform the following controls:
- the continental shelf profile of the target work area is determined according to the relationship.
- the apparatus for determining the continental shelf profile based on non-delta transgression may be configured separately from the central processing unit 100, for example, the apparatus for determining the continental shelf profile based on non-delta transgression may be configured as a chip connected to the central processing unit 100, The control function of the continental shelf profile determination method based on non-delta transgression is realized through the control of the central processor.
- the central processing unit 100 may include a microprocessor or other processor device and/or logic device, and the central processing unit 100 receives input and controls various aspects of the electronic device 600 component operation.
- the memory 140 may be, for example, one or more of a cache, a flash memory, a hard drive, a removable medium, a volatile memory, a non-volatile memory or other suitable devices.
- the above-mentioned information related to the failure can be stored, and a program executing the related information can also be stored.
- the central processing unit 100 can execute the program stored in the memory 140 to realize information storage or processing.
- the input unit 120 provides input to the central processing unit 100 .
- the input unit 120 is, for example, a key or a touch input device.
- the power supply 170 is used to provide power to the electronic device 600 .
- the display 160 is used for displaying display objects such as images and characters.
- the display can be, for example, but not limited to, an LCD display.
- the memory 140 may be solid state memory, eg, read only memory (ROM), random access memory (RAM), SIM card, and the like. There may also be memories that retain information even when powered off, selectively erased and provided with more data, examples of which are sometimes referred to as EPROMs or the like. Memory 140 may also be some other type of device. Memory 140 includes buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application/function storage part 142 for storing application programs and function programs or for performing operations of the electronic device 600 through the central processing unit 100 .
- the memory 140 may also include a data store 143 for storing data such as contacts, digital data, pictures, sounds and/or any other data used by the electronic device.
- the driver storage section 144 of the memory 140 may include various drivers of the electronic device for communication functions and/or for executing other functions of the electronic device (eg, a messaging application, a contact book application, etc.).
- the communication module 110 is the transmitter/receiver 110 that transmits and receives signals via the antenna 111 .
- a communication module (transmitter/receiver) 110 is coupled to the central processing unit 100 to provide input signals and receive output signals, as may be the case with conventional mobile communication terminals.
- multiple communication modules 110 may be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and/or a wireless local area network module.
- Communication module (transmitter/receiver) 110 is also coupled to speaker 131 and microphone 132 via audio processor 130 to provide audio output via speaker 131 and to receive audio input from microphone 132 for general telecommunication functions.
- Audio processor 130 may include any suitable buffers, decoders, amplifiers, and the like.
- the audio processor 130 is also coupled to the central processing unit 100 so as to enable recording on the local unit through the microphone 132 and to play back the sound stored on the local unit through the speaker 131 .
- each module may be implemented in the same one or more software and/or hardware, and the modules implementing the same function may be implemented by a combination of multiple sub-modules or sub-units.
- the apparatus embodiments described above are only illustrative.
- the division of the units is only a logical function division.
- multiple units or components may be combined or integrated. to another system, or some features can be ignored, or not implemented.
- the shown or discussed mutual coupling or direct coupling or communication connection may be through some interfaces, indirect coupling or communication connection of devices or units, and may be in electrical, mechanical or other forms.
- the controller in addition to implementing the controller in the form of pure computer-readable program code, the controller can be implemented as logic gates, switches, application-specific integrated circuits, programmable logic controllers and embedded devices by logically programming the method steps.
- the same function can be realized in the form of a microcontroller, etc. Therefore, this kind of controller can be regarded as a kind of hardware component, and the devices included in it for realizing various functions can also be regarded as a structure in the hardware component. Or even, the means for implementing various functions can be regarded as both a software module implementing a method and a structure within a hardware component.
- Memory may include forms of non-persistent memory, random access memory (RAM) and/or non-volatile memory in computer readable media, such as read only memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
- RAM random access memory
- ROM read only memory
- flash RAM flash memory
- Computer-readable media includes both persistent and non-permanent, removable and non-removable media, and storage of information may be implemented by any method or technology.
- Information may be computer readable instructions, data structures, modules of programs, or other data.
- Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read only memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), Flash Memory or other memory technology, Compact Disc Read Only Memory (CD-ROM), Digital Versatile Disc (DVD) or other optical storage, Magnetic tape cartridges, magnetic tape disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information that can be accessed by a computing device.
- computer-readable media does not include transitory computer-readable media, such as modulated data signals and carrier waves.
- Embodiments of this specification may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer.
- program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types.
- Embodiments of the description may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network.
- program modules may be located in both local and remote computer storage media including storage devices.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于非三角洲海侵的大陆架剖面确定方法及装置，基于非三角洲海侵的大陆架剖面确定方法包括：根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；根据所述关系式确定所述目的工区的大陆架剖面。本发明提供的基于非三角洲海侵的大陆架剖面确定方法及装置，在非三角洲海侵背景下，综合考虑到河流沉积物供给以及海平面上升对大陆架剖面的影响，从而更加精准的确定目的工区历史时期中大陆架剖面的变化情况。
Description
本发明涉及石油天然气勘探技术领域与海洋地质领域，尤其是地层叠加样式分析技术领域，具体涉及一种基于非三角洲海侵的大陆架剖面确定方法及装置。
现有技术中，大陆架剖面计算是基于波浪对海底的侵蚀作用，具体地，其基本原理为，大陆架从深海方向向大陆方向水深逐渐变浅，而随着水深变浅，波浪对海底沉积物的作用加强，沉积物遭受破坏越强烈。在浅水区被破坏的沉积物将被搬运到相对深水位置沉积下来。这个过程最终将达到一个平衡状态，在这种情况下大陆架的剖面也被叫做平衡剖面。因此大陆架的平衡剖面应该具有向陆地方向变陡的特征。Bruun将这一特征扩展到了海平面上升背景，称为布鲁恩法则(Bruun,1962)。其基本原理为，由于海平面上升，原为陆地的部分被淹没，成为浅水区，因此早期的陆地沉积物遭受剥蚀，被搬运至深水区，形成新的平衡剖面。因此，在海平面上升背景下，新生成的大陆架仍然表现为向陆地方向变陡的特征(图1)。上述大陆架剖面计算方法存在以下缺点：
(1)现有技术适用于较小的时间尺度，在数十至数百年间，属于工程学的领域，并且假设地质构造、沉积物供给等外部因素不发生变化。而对于数千至数十万年的时间尺度，外部条件也在缓慢地变化过程中，上述技术很难适用。
(2)现有技术没有考虑陆地方向的沉积物供给作用。大陆架的产生通常是海平面上升所致。海平面上升使得陆地淹没，形成大陆架。而在海岸线附近，除了海平面上升趋向于使岸线向陆地方向迁移之外，常常同时河流作用趋向于使岸线向海洋方向迁移，现有的方案并没有考虑后者的作用。
发明内容
针对现有技术中的问题，本发明提供的基于非三角洲海侵的大陆架剖面确定方法及装置，在非三角洲海侵背景下，综合考虑到河流沉积物供给以及海平面上升对大陆架剖面的影响，从而更加精准的确定目的工区历史时期中大陆架剖面的变化情况。
第一方面，本发明提供一种基于非三角洲海侵的大陆架剖面确定方法，包括：
根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；
在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；
根据所述关系式确定所述目的工区的大陆架剖面。
一实施例中，所述根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型，包括：
计算所述目的工区单位时间内的冲积范围的纵向增量；
根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；
根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
一实施例中，在所述大陆架剖面二维模型中：
内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及
在海平面上升开始，河流的具有初始长度。
一实施例中，所述在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式，包括：
在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
第二方面，本发明提供一种基于非三角洲海侵的大陆架剖面确定装置，该装置包括：
模型建立单元，用于根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；
模型求解单元，用于在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；
剖面确定单元，用于根据所述关系式确定所述目的工区的大陆架剖面。
一实施例中，所述模型建立单元包括：
纵向增量计算模块，用于计算所述目的工区单位时间内的冲积范围的纵向增量；
沉积物体量计算模块，用于根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；
河流供给体量计算模块，用于根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
一实施例中，在所述大陆架剖面二维模型中：
内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及
在海平面上升开始，河流的具有初始长度。
一实施例中，所述模型求解单元具体用于在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现基于非三角洲海侵的大陆架剖面确定方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现基于非三角洲海侵的大陆架剖面确定方法的步骤。
从上述描述可知，本发明实施例提供的基于非三角洲海侵的大陆架剖面确定方法及装置，首先根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；接着，在空间直角坐标系下，求解大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；最后根据关系式确定目的工区的大陆架剖面。本发明提在非三角洲海侵背景下，综合考虑到河流沉积物供给以及海平面上升对大陆架剖面的影响，从而更加精准的确定目的工区历史时期中大陆架剖面的变化情况。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明背景技术中基于布鲁恩法则的大陆架平衡剖面示意图；
图2为本发明的实施例中的基于非三角洲海侵的大陆架剖面确定方法的流程示意图；
图3为本发明的实施例中步骤100的流程示意图；
图4为本发明的实施例中时间增量Δt内形成的冲积范围的纵向面积增量示意图；
图5为本发明的实施例中步骤200的流程示意图；
图6为本发明的具体应用实例中非三角洲海侵剖面求解示意图；
图7为本发明的具体应用实例中大陆架坡度(φ)与无量纲的河流长度(L
*)关系图；
图8为本发明的具体应用实例中x-z坐标体系设置示意图；
图9为本发明的具体应用实例中基于非三角洲海侵的大陆架剖面确定方法的流程示意图；
图10为本发明的实施例中的基于非三角洲海侵的大陆架剖面确定装置的结构示意图；
图11为本发明的实施例中基于非三角洲海侵的大陆架剖面确定装置中模型建立单元10组成结构示意图；
图12为本发明实施例的电子设备600的系统构成的示意框图。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明的实施例提供一种基于非三角洲海侵的大陆架剖面确定方法的具体实施方式，参见图2，该方法具体包括如下内容：
步骤100：根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型。
可以理解的是，目的工区处于是非三角洲海侵的沉积背景下，另外，需要说明的是，当海平面向陆地方向后退，才称之为海侵，才能够产生新的大陆架。海侵可以进一步分为两种情况：三角洲海侵与非三角洲海侵。河流的末端伴随三角洲的发育称之为三角洲海侵；河流的末端不发育三角洲，为非三角洲海侵，这取决于河流的长度。当冲积河流的长度(L)大于临界长度L
crt(L>L
crt)时为非三角洲海侵；反之，若L<L
crt，为三角洲海侵。
步骤200：在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非 三角洲海侵坡度以及河流长度的关系式。
具体地，将大陆架置于x-z坐标系中(x代表水平方向，z代表垂直方向)，大陆架剖面在某一位置的坡度，相当于大陆架剖面在该位置处的切线斜率。因此，在该x-z坐标系中，可以找到非三角洲海侵剖面的解决方案。
步骤300：根据所述关系式确定所述目的工区的大陆架剖面。
因为步骤300中的关系式可以表征非三角洲海侵坡度，故其可以确定目的工区的大陆架剖面。
从上述描述可知，本发明实施例提供的基于非三角洲海侵的大陆架剖面确定方法，首先根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；接着，在空间直角坐标系下，求解大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；最后根据关系式确定目的工区的大陆架剖面。本发明提在非三角洲海侵背景下，综合考虑到河流沉积物供给以及海平面上升对大陆架剖面的影响，从而更加精准的确定目的工区历史时期中大陆架剖面的变化情况。
一实施例中，在所述大陆架剖面二维模型中：
内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及
在海平面上升开始，河流的具有初始长度。
可以理解的是，为了找到非三角洲海侵坡度(即大陆架坡度φ)与河流长度(L)之间的定量关系，使用二维几何模型来阐明随着海平面上升，河流长度变化(L>>L
crt)与大陆架坡度变化的过程。为了简化，在建立大陆架剖面二维模型的过程中，采用以下条件：
(1)内陆基底坡度(γ)和冲积坡度(α)保持恒定。
(2)上游的沉积物供给速率(q
s，二维空间中为单位宽度内的沉积物体积供给速率)保持恒定(即q
s＝定值)。
(3)在海平面开始上升至结束的过程中，河流的长度(L)比始终大于临界值(即L>>L
crt)，即海侵始终保持非三角洲海侵。
(4)海平面上升的速率保持恒定(R
sl＝定值>0)。
(5)在海平面上升开始，河流的具有初始长度(L
0；L
0>>L
crt)。
一实施例中，参见图3，步骤100进一步包括：
步骤101：计算所述目的工区单位时间内的冲积范围的纵向增量；
参见图4，计算在时间增量Δt内形成的冲积范围的纵向增量。通过近似处理，可以将空间增量视为梯形形状，可以得出梯形区域S
1的表达式：
步骤102：根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；
具体地，计算Δt期间内供给的沉积物数量，用S
2来表示，即：
S
2＝q
sΔt， (2)
步骤103：根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
一实施例中，参见图5，步骤200进一步包括：
步骤201：在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
可以理解的是，空间任意选定一点O,过点O作三条互相垂直的数轴x，y，z，它们都以O为原点且具有相同的长度单位。这三条轴分别称作横轴，纵轴，竖轴，该坐标系称为空间直角坐标系。
为进一步地说明本方案，本发明还提供基于非三角洲海侵的大陆架剖面确定方法的具体应用实例，该具体应用实例具体包括如下内容，参见图9。
为克服现有技术中的各种问题，在本具体应用实例中，同时考虑海平面上升与河流作用，即在两者的共同作用下，确定大陆架产生的条件及理论剖面形态。
可以理解的是，在海平面上升过程中，海岸线既可以向海洋方向推进，也可以向陆地方向后退，还可以在原位置保持不动。这一过程取决于河流供给沉积物与海平面上升两者的相对强弱。海平面上升将产生多余的可供沉积物充填的空间，当河流搬运来的沉积物足够多，在填充了这部分空间之外还有剩余的沉积物，则剩余的沉积物将进行造陆，海岸线向海洋方向推进；反之，当河流搬运来的沉积物不足以填充由海平面上升增加的空间，则海平面将向陆地方向后退；如果二者相当，则海岸线在原位置固定不变。只有当海平面向陆地方向后退，才称之为海侵，才能够产生新的大陆架。海侵可以进一步分为两种情况：三角洲海侵与非三角洲海侵。河流的末端伴随三角洲的发育称之为三角洲海侵；河流的末端不发育三角洲，为非三角洲海侵，这取决于河流的长度。研究表明，当冲积河流的长度(L)大于临界长度L
crt(L>L
crt)时为非三角洲海侵；反之，若L<L
crt，为三角洲海侵(Tomer et al.,2011；Wang et al.,2019)。可以想象，对于一个大型河流-三角洲体系(L>L
crt)，在海平面上升一般应表现为非三角洲海侵。随着海侵的继续与河流长度的缩短，当达到L<L
crt时，三角洲将重新形成(Tomer et al.,2011)。
在本具体应用实例中，仅针对非三角洲海侵的情况(图6)。图6是非三角洲海侵剖面求解示意。R
sl为海平面变化速率(>0为上升)，q
s为沉积物供给速率。L为河流的长度，随海平面上升将变短，φ为与L相应的大陆架坡度，随着L的变短将变大。α为河流的坡度，γ为陆地基底的坡度。ΔH为海平面上升的微分单元。
对于非三角洲海侵，所有的沉积物全部沉积在海平面之上。假定沉积物供给速率恒定，则沉积物在河流表面的加积速率与河流的长度成反比。如果河流长度非常长(L>>L
crt)，则沿冲积剖面的加积速率将非常低，海岸线附近的沉积物聚积受到限制，海岸线能够更“自由”地向陆地方向迁移。相比之下，如果河流长度比较短(但仍然L>L
crt)，海岸线附近的沉积物聚积速率变大，导致陆架更陡峭，海岸线向陆地方向的迁移也会受到限制(变缓)。因此，河流在向陆退积和长度(L)缩短之后，非三角洲海侵剖面的形态基本上会在上游变陡(如图6所示)。
S1：建立大陆架剖面二维模型。
在公式(2)的基础上，当Δt间隔内由海平面上升形成的新增大陆架可以线性近似(即新增的大陆架坡度φ＝定值)时，存在S
1≈S
2的关系。此时可以得到：
公式(3)两侧同时除以R
blrq
sΔt，可得：
其中，令：
Λ
2D为二维空间中海平面升降背景下的空间尺度度量，具有长度单位。
则：
接着，利用微分方法求取海侵剖面坡度。实际上，仅当Δt无限接近0时，S
1≈S
2才成立。通过这样做，可以将以上方程式缩写为：
或者：
以上等式是φ的解析解：它取决于L和Λ
2D，适用于非三角洲海侵(即L>L
crt)。
将公式(8)等号右侧分子与分母同时除以Λ
2D，即将河流长度L进行无量纲化，得到公式(9)。公式(9)为任意非三角洲海侵剖面坡度与河流长度(无量纲化)的一般形式。
其中上标*表示无量纲的量。
参见图7，φ与L
*相关关系图，从图中可以看出，φ与L
*呈反相关关系。从图7中可以看出，
①如果L(或L
*)非常大，φ趋近于α。大陆架剖面和河流剖面一致。
②随着L(或L
*)的缩短，陆架坡度(φ)变陡，并且陆架剖面在海侵过程中逐渐远离原有的河流剖面。
以上等式为物理解释提供了一种数学方案：因为河流加积速率与河长之间呈反比关系，因此，较长的冲积河中的海平面上升将无法有效地将洪泛面从已有的河流剖面上移开。即：随着海平面上升，形成的大陆架表现为向陆地一侧变陡，向海洋一侧变缓，且向海洋一侧无线接近于早期暴露的河流的坡度。
S2：求解大陆架剖面二维模型。
首先对海侵剖面坡度表达式进行形式转换。从数学意义上讲，如果将大陆架置于x-z坐标系中(x代表水平方向，z代表垂直方向)，大陆架剖面在某一位置的坡度，相当于大陆架剖面在该位置处的切线斜率(即φ)。因此，在该x-z坐标系中，可能找到非三角洲海侵剖面的解决方案。
设定x–z坐标系的原点是在t＝0(海平面刚刚开始上升)时穿过陆架边缘点的水平线与内陆剖面的交点(图8)。在此指定的坐标系内，L具有以下形式：
将公式(10)代入公式(8)，可以将φ由x–z坐标以及α,γ和Λ
2D表示：
图7是适用于本发明的x-z坐标体系设置(x代表水平方向，向海洋方向为正；z代表垂直方向，向上为正)示意图。注意：坐标原点是在t＝0(海平面刚刚开始上升)时 穿过陆架边缘点的水平线与内陆剖面的交点。
接着，确定海侵剖面的表达形式。由于φ代表了大陆架剖面在任意位置处的切线斜率，公式(11)代表了海侵剖面的导函数。即海侵剖面应该为公式(11)所示的原函数。假设η(x,z)是原函数，则η(x,z)可以表示为公式(11)自初始位置(x
0,z)至(x,z)的积分：
其中(x
0，0)和(x
1，z
1)分别表示初始和最终海岸线位置。具体来说，在最初的海岸线处，有以下关系式：
其中L
0为t＝0(海平面刚刚开始上升)时的河流长度。在最后的海岸线(x
1，z
1)处，仍然满足不等式：L≥L
crt.。公式(12)没有解析解，但可以通过数值模拟求得数值解。
坐标系与海侵剖面表达式的无量纲化：公式(11)–(13)也可以将等号右侧分子与分母同除以Λ
2D，从而得到无量纲化海侵剖面的一般表达形式。
其中上标*表示无量纲的量。
从上述描述可知，本发明实施例提供的基于非三角洲海侵的大陆架剖面确定方法，首先基于管道的内检测数据选取自变量，并使用Lasso算法对输入变量进行选择，最后对管道内腐蚀缺陷预测采用以内检测数据为驱动的广义线性可加模型，本发明实施例可基于现有的两次内腐蚀数据进行建模对未来的腐蚀深度进行估算。由于模型的分析结果可以为管道的危险缺陷提供参考依据，因此可以是完整性管理的重要部分，可以帮助确定内检测周期以及制定维修计划，有利于管道安全运行。具体地，本发明实施例的有益效果如下：
基于同一发明构思，本申请实施例还提供了基于非三角洲海侵的大陆架剖面确定装 置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于基于非三角洲海侵的大陆架剖面确定装置解决问题的原理与基于非三角洲海侵的大陆架剖面确定方法相似，因此基于非三角洲海侵的大陆架剖面确定装置的实施可以参见基于非三角洲海侵的大陆架剖面确定方法实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
本发明的实施例提供一种能够实现基于非三角洲海侵的大陆架剖面确定方法的基于非三角洲海侵的大陆架剖面确定装置的具体实施方式，参见图10，基于非三角洲海侵的大陆架剖面确定装置具体包括如下内容：
模型建立单元10，用于根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；
模型求解单元20，用于在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；
剖面确定单元30，用于根据所述关系式确定所述目的工区的大陆架剖面。
一实施例中，参见图11，所述模型建立单元10包括：
纵向增量计算模块101，用于计算所述目的工区单位时间内的冲积范围的纵向增量；
沉积物体量计算模块102，用于根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；
河流供给体量计算模块103，用于根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
一实施例中，在所述大陆架剖面二维模型中：
内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及
在海平面上升开始，河流的具有初始长度。
一实施例中，所述模型求解单元20具体用于在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
从上述描述可知，本发明实施例提供的基于非三角洲海侵的大陆架剖面确定装置，首先根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；接着，在空间直角坐标系下，求解大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；最后根据关系式确定目的工区的大陆架剖面。本发明提在非 三角洲海侵背景下，综合考虑到河流沉积物供给以及海平面上升对大陆架剖面的影响，从而更加精准的确定目的工区历史时期中大陆架剖面的变化情况。
本申请的实施例还提供一种电子设备，该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照上述实施例的方法的实施及上述实施例所述的装置，其内容被合并于此，重复之处不再赘述。
图12为本发明实施例的电子设备600的系统构成的示意框图。如图12所示，该电子设备600可以包括中央处理器100和存储器140；存储器140耦合到中央处理器100。值得注意的是，该图是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
一实施例中，基于非三角洲海侵的大陆架剖面确定方法可以被集成到中央处理器100中。其中，中央处理器100可以被配置为进行如下控制：
根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；
在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；
根据所述关系式确定所述目的工区的大陆架剖面。
在另一个实施方式中，基于非三角洲海侵的大陆架剖面确定装置可以与中央处理器100分开配置，例如可以将基于非三角洲海侵的大陆架剖面确定装置配置为与中央处理器100连接的芯片，通过中央处理器的控制来实现基于非三角洲海侵的大陆架剖面确定方法的管控功能。
如图12所示，该电子设备600还可以包括：通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是，电子设备600也并不是必须要包括图12中所示的所有部件；此外，电子设备600还可以包括图12中没有示出的部件，可以参考现有技术。
如图12所示，中央处理器100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中，存储器140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序，以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142，该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143，该数据存储部143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132，以经由扬声器131提供音频输出，并接收来自麦克风132的音频输入，从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器130还耦合到中央处理器100，从而使得可以通过麦克风132能够在本机上录音，且使得可以通过扬声器131来播放本机上存储的声音。
本申请的实施例还提供能够实现上述实施例中的基于非三角洲海侵的大陆架剖面确定方法中全部步骤的一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的基于非三角洲海侵的大陆架剖面确定方法的全部步骤，例如，处理器执行计算机程序时实现下述步骤：
根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；
在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；
根据所述关系式确定所述目的工区的大陆架剖面。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接 可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。 计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的 精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (10)
- 一种基于非三角洲海侵的大陆架剖面确定方法，其特征在于，包括：根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；根据所述关系式确定所述目的工区的大陆架剖面。
- 根据权利要求1所述的基于非三角洲海侵的大陆架剖面确定方法，其特征在于，所述根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型，包括：计算所述目的工区单位时间内的冲积范围的纵向增量；根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
- 根据权利要求2所述的基于非三角洲海侵的大陆架剖面确定方法，其特征在于，在所述大陆架剖面二维模型中：内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及在海平面上升开始，河流的具有初始长度。
- 根据权利要求1所述的基于非三角洲海侵的大陆架剖面确定方法，其特征在于，所述在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式，包括：在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
- 一种基于非三角洲海侵的大陆架剖面确定装置，其特征在于，包括：模型建立单元，用于根据目的工区的河流供给沉积物体量以及海平面上升速率建立大陆架剖面二维模型；模型求解单元，用于在空间直角坐标系下，求解所述大陆剖面二维模型，以生成可以表征非三角洲海侵坡度以及河流长度的关系式；剖面确定单元，用于根据所述关系式确定所述目的工区的大陆架剖面。
- 根据权利要求5所述的基于非三角洲海侵的大陆架剖面确定装置，其特征在于，所述模型建立单元包括：纵向增量计算模块，用于计算所述目的工区单位时间内的冲积范围的纵向增量；沉积物体量计算模块，用于根据所述纵向增量以及沉积物体积供给速率计算单位时间内的沉积物体量；河流供给体量计算模块，用于根据单位时间内的沉积物体量计算所述河流供给沉积物体量。
- 根据权利要求6所述的基于非三角洲海侵的大陆架剖面确定装置，其特征在于，在所述大陆架剖面二维模型中：内陆基底坡度、冲积坡度、河流上游的沉积物供给速率以及海平面上升速率均为一恒定值；以及在海平面上升开始，河流的具有初始长度。
- 根据权利要求5所述的基于非三角洲海侵的大陆架剖面确定装置，其特征在于，所述模型求解单元具体用于在x-z坐标系中，对所述大陆剖面二维模进行求解，以生成所述关系式。
- 一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述基于非三角洲海侵的大陆架剖面确定方法的步骤。
- 一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至4任一项所述基于非三角洲海侵的大陆架剖面确定方法的步骤。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110458579.6 | 2021-04-27 | ||
| CN202110458579.6A CN113094929B (zh) | 2021-04-27 | 2021-04-27 | 基于非三角洲海侵的大陆架剖面确定方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2022228450A1 true WO2022228450A1 (zh) | 2022-11-03 |
Family
ID=76680148
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2022/089455 WO2022228450A1 (zh) | 2021-04-27 | 2022-04-27 | 基于非三角洲海侵的大陆架剖面确定方法及装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN113094929B (zh) |
| WO (1) | WO2022228450A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113094929B (zh) * | 2021-04-27 | 2024-02-06 | 中国石油大学(北京) | 基于非三角洲海侵的大陆架剖面确定方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109190538A (zh) * | 2018-08-24 | 2019-01-11 | 华北水利水电大学 | 一种基于遥感技术的多泥沙河流三角洲海岸带演化分析方法 |
| CN113094929A (zh) * | 2021-04-27 | 2021-07-09 | 中国石油大学(北京) | 基于非三角洲海侵的大陆架剖面确定方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108629130B (zh) * | 2018-05-10 | 2019-09-24 | 山东科技大学 | 确定控制陆源分散有机质分布三角洲体系外缘坡度的方法 |
-
2021
- 2021-04-27 CN CN202110458579.6A patent/CN113094929B/zh active Active
-
2022
- 2022-04-27 WO PCT/CN2022/089455 patent/WO2022228450A1/zh active Application Filing
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109190538A (zh) * | 2018-08-24 | 2019-01-11 | 华北水利水电大学 | 一种基于遥感技术的多泥沙河流三角洲海岸带演化分析方法 |
| CN113094929A (zh) * | 2021-04-27 | 2021-07-09 | 中国石油大学(北京) | 基于非三角洲海侵的大陆架剖面确定方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| LIU JIANPING; XIAN BENZHONG; JI YOULIANG; GONG CHENGLIN; WANG JUNHUI; WANG ZHEN; CHEN PENG; SONG DELIN; WEI WENZONG; ZHANG XIMENG;: "Alternating of aggradation and progradation dominated clinothems and its implications for sediment delivery to deep lake: The Eocene Dongying Depression, Bohai Bay Basin, east China", MARINE AND PETROLEUM GEOLOGY, ELSEVIER, AMSTERDAM, NL, vol. 114, 19 December 2019 (2019-12-19), AMSTERDAM, NL , XP086013734, ISSN: 0264-8172, DOI: 10.1016/j.marpetgeo.2019.104197 * |
| WANG JUNHUI, TETSUJI MUTO: "Autostratigraphic modelling of the growth of alluvial-shelf systems during steady base-level cycles: Two-dimensional tank experiments", SEDIMENTOLOGY, vol. 68, no. 1, 9 August 2020 (2020-08-09), pages 135 - 167, XP055981987, DOI: 10.1111/sed.12789 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113094929A (zh) | 2021-07-09 |
| CN113094929B (zh) | 2024-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Baldock et al. | Morphological hysteresis in the evolution of beach profiles under sequences of wave climates-Part 1; observations | |
| Yoshida et al. | Changes in depositional processes—an ingredient in a new generation of sequence-stratigraphic models | |
| WO2022228450A1 (zh) | 基于非三角洲海侵的大陆架剖面确定方法及装置 | |
| Myrow et al. | Bedform disequilibrium | |
| Chadwick et al. | Origin of a preferential avulsion node on lowland river deltas | |
| Hutton et al. | Sedflux 2.0: An advanced process-response model that generates three-dimensional stratigraphy | |
| WO2017069860A1 (en) | System and method for modeling coupled systems of hydrodynamics and sediment transport | |
| Carlson et al. | Control of basin water depth on channel morphology and autogenic timescales in deltaic systems | |
| Hampel et al. | Response of the tectonically erosive south Peruvian forearc to subduction of the Nazca Ridge: Analysis of three‐dimensional analogue experiments | |
| Hsu et al. | An experimental and numerical investigation on wave‐mud interactions | |
| Wu et al. | Impacts of backwater hydrodynamics on fluvial–deltaic stratigraphy | |
| Muto et al. | Planform evolution of deltas with graded alluvial topsets: Insights from three‐dimensional tank experiments, geometric considerations and field applications | |
| Anderson et al. | A geomorphic enthalpy method: Description and application to the evolution of fluvial-deltas under sea-level cycles | |
| CN107656312B (zh) | 基于分方位角叠加的河道砂体预测方法及装置 | |
| CN115480296A (zh) | 天然气水合物稳定带底界预测方法及装置 | |
| Burgess et al. | How to interpret, understand, and predict stratal geometries using stratal-control spaces and stratal-control-space trajectories | |
| CN113433587B (zh) | 岩石物理反演方法及系统 | |
| MUTO et al. | The autostratigraphic view of responses of river deltas to external forcing: A review of the concepts | |
| Olariu et al. | Ubiquity of wave-dominated deltas in outer-shelf growth-faulted compartments | |
| Matsumoto et al. | Modelling the development of varied shore profile geometry on rocky coasts | |
| CN112926805B (zh) | 基于深度强化学习的直井试井智能解释方法及装置 | |
| CN113281825A (zh) | 岩石物理模型构建方法及装置 | |
| CN109063228B (zh) | 一种水驱油藏含水上升率变化的确定方法及装置 | |
| CA3158470A1 (en) | Method and apparatus for predicting hydrocarbon reservoir with two-phase medium | |
| CN112596102B (zh) | 油藏渗流特征参数的地球物理预测方法、装置及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22794920 Country of ref document: EP Kind code of ref document: A1 |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 18288541 Country of ref document: US |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 22794920 Country of ref document: EP Kind code of ref document: A1 |