River regulation scheme generation method and system based on river regime evolution map

Abstract

The application discloses a riverway management scheme generation method and system based on a river regime evolution map, belongs to the technical field of computer-aided design and simulation analysis, and comprises the following steps: position registration and time sorting of multi-period riverway observation data to generate a river regime state sequence; identification of a river regime change unit and connection to generate a river regime evolution map; mapping of a management target and extraction of a contradiction subgraph under engineering constraints; embedding of a candidate management unit to generate a response path through counterfactual deduction; combination of identified action convergence, offset and amplification to generate a combined response structure; tailoring to generate an implementable scheme map and analysis to generate a riverway management scheme. The application can uniformly express the causal conduction relationship of the river regime along a river section, quantitatively evaluate the synergy and conflict effect of various management actions, and correct the map according to monitoring data after implementation to continuously iteratively optimize the scheme.

Description

Technical Field

[0001] This invention relates to the field of computer-aided design and simulation analysis technology, and in particular to a method and system for generating river management schemes based on river morphology evolution maps. Background Technology

[0002] River management refers to engineering activities that intervene in the evolution of river morphology to maintain the flood control capacity of rivers and stabilize the shoreline and navigation channel. River morphology refers to the overall changes in the river's planar shape, channel cross-section, and main channel position as a result of water flow and sediment deposition. Its evolution is characterized by segment-by-segment transmission along the direction of water flow, and changes in one river segment often affect the responses of related upstream and downstream segments. With the development of remote sensing observation, underwater topographic surveying, and computer modeling technologies, using computers to process river observation data and assist in the design of management schemes has become an important means of water conservancy engineering planning.

[0003] Existing computer-aided river management schemes mostly employ hydrodynamic and sediment mathematical models to perform forward numerical simulations of given working conditions. By setting boundary conditions, they solve the equations for water flow and sediment transport to obtain the riverbed scouring and deposition distribution. Then, designers manually select management measures based on the simulation results. Another type of scheme uses statistical regression or trend extrapolation of historical observation data to predict changes in the shoreline or riverbed at specific cross-sections, thereby determining whether measures such as bank protection, groynes, or dredging are necessary.

[0004] The above approach treats each section or working condition of the river channel in relative isolation, making it difficult to express the causal transmission relationship between river sections and the changes in river morphology. It also lacks a unified structured expression for the conflicts between the governance objectives and the related river sections. Forward numerical simulation requires resetting the working conditions for each alternative measure and solving it repeatedly, which is computationally expensive. The mutual influence of multiple measures implemented at the same time is difficult to quantify and assess during the design phase. Statistical extrapolation relies on the continuation of existing trends and lacks the ability to extrapolate the evolutionary changes caused by engineering intervention. Furthermore, the actual response after the implementation of the project is difficult to feed back to the original model for continuous correction. Summary of the Invention

[0005] To address the aforementioned issues, this invention provides a method and system for generating river management schemes based on river morphology evolution maps. By constructing a river morphology evolution map with causal transmission edges from multi-period river observation data, mapping management objectives and engineering constraints on the map to extract contradictory subgraphs, embedding candidate management units for counterfactual deduction and combined evaluation, and then correcting the map based on post-implementation monitoring data, this method can uniformly express the causal transmission relationship of the river morphology along the river section, quantitatively evaluate the synergistic and conflicting effects of various management actions, and enable the management scheme to be continuously iterated and optimized with actual response.

[0006] The above objectives can be achieved through the following approach:

[0007] A method for generating river management schemes based on river morphology evolution maps includes: acquiring multi-period river observation data of a target river; performing location registration and temporal sorting on the multi-period river observation data to generate a river morphology state sequence; identifying river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence; connecting the river morphology change units with sequential response relationships to generate a river morphology evolution map; acquiring the management objectives and engineering constraints of the target river; mapping the management objectives and engineering constraints to the river morphology evolution map to extract contradictions and generate a contradiction sub-map; and based on the river morphology change unit transmission method in the contradiction sub-map... Candidate governance units are constructed based on the degree of transmission or propagation. These candidate governance units are then embedded one by one into the contradictory subgraph, and the changes in the graph before and after embedding are compared to generate counterfactual response paths. The counterfactual response paths with spatial and temporal connections are combined, and the convergence, cancellation, and amplification of effects after combination are identified to generate a combined response structure. The engineering constraints are used to trim the candidate governance units and their corresponding connections in the combined response structure that cannot be implemented together, generating a feasible implementation diagram. The effect transmission order and construction dependency order of the candidate governance units in the feasible implementation diagram are analyzed to generate a river management scheme.

[0008] Optionally, the step of acquiring multi-period river observation data of the target river, performing location registration and time sorting on the multi-period river observation data, and generating a river morphology sequence includes: acquiring multi-period river observation data of the target river, extracting shoreline location data, channel elevation data, sandbar boundary data, and main channel location data from the multi-period river observation data, and generating phased river morphology data; correcting the location deviations and observation scale differences between the phased river morphology data, and generating river morphology data at the same scale; dividing the continuous river segments in the river morphology data at the same scale and establishing a time index for each continuous river segment, and generating segmented time-series morphology data; connecting adjacent temporal states belonging to the same continuous river segment in the segmented time-series morphology data, and generating a river morphology sequence.

[0009] Optionally, identifying river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connecting river morphology change units with sequential response relationships, and generating a river morphology evolution map includes: comparing the spatial differences and change directions of adjacent temporal states in the river morphology state sequence to generate river morphology change segments; merging spatially connected change segments with consistent change directions and continuous phases in the river morphology change segments to generate river morphology change units; tracking the occurrence order and response direction of the river morphology change units in the direction of water flow connectivity to generate candidate transitive relationships; eliminating relationships in the candidate transitive relationships that lack continuous temporal responses or have mutually exclusive evolution directions to generate causal transitive edges; and connecting the river morphology change units with the causal transitive edges to generate a river morphology evolution map.

[0010] Optionally, the step of obtaining the governance objectives and engineering constraints of the target river channel, mapping the governance objectives and engineering constraints to the river morphology evolution map for contradiction extraction, and generating a contradiction subgraph includes: converting the governance objectives into the target evolution direction of the river morphology change unit, generating target state markers; converting the engineering constraints into constraint markers that restrict the range and manner of change of the river morphology change unit, generating map constraint boundaries; using the target state markers to trace the causal transmission edges backward from the river morphology evolution map, generating target action chains; using the map constraint boundaries to identify the opposite evolution direction formed when the target action chain continues to propagate to the associated river segment, generating contradiction action pairs; and extracting the river morphology change unit and causal transmission edge that connects the contradiction action pairs and does not contain irrelevant branches, generating a contradiction subgraph.

[0011] Optionally, the step of constructing candidate governance units based on the transmission direction or degree of river morphology change units in the contradictory subgraph, embedding the candidate governance units one by one into the contradictory subgraph, and comparing the changes in the graph before and after embedding to generate a counterfactual response path includes: identifying the spatial location or evolution stage in the contradictory subgraph that can separate the contradictory action pair, generating action separation location; matching governance actions that can change the state of water flow guidance, shoreline constraints, channel flow, or sandbar evolution at the action separation location, generating candidate governance units; embedding the candidate governance units into the action separation location and changing the causal transmission edges connected to the corresponding action separation location, generating an intervention graph branch; comparing the transmission direction and degree of the intervention graph branch with the corresponding node in the contradictory subgraph, generating a response difference chain; and extracting the continuous portion of the response difference chain extending from the candidate governance unit to the target improved state or the unfavorable state of the associated river section, generating a counterfactual response path.

[0012] Optionally, the combination of counterfactual response paths with spatial and temporal connections, identifying the convergence, cancellation, and amplification of effects after combination, and generating a combined response structure includes: arranging the order in which the counterfactual response paths arrive at the same river morphology change unit to generate a path temporal relationship; comparing the evolution direction of arrival at the same river morphology change unit based on the path temporal relationship to generate a path interaction relationship; identifying path combinations with the same evolution direction and increasing transmission degree, path combinations with opposite evolution directions and decreasing transmission degree, and path combinations ending transmission at the same node in the path interaction relationship to generate amplification, cancellation, and convergence relationships; connecting candidate governance units with the amplification, cancellation, or convergence relationship, and retaining the corresponding path temporal relationship to generate a combined response structure.

[0013] Optionally, the step of using the engineering constraints to trim the candidate governance units and corresponding connection relationships in the combined response structure that cannot be implemented together, and generating a feasible implementation diagram, includes: mapping the engineering constraints to the candidate governance units and path timing relationships in the combined response structure to generate governance unit constraint relationships; determining whether the governance unit constraint relationships have spatial occupancy conflicts, implementation timing conflicts, or situations where the effect is amplified and crosses the map constraint boundary, and generating conflict governance combinations; removing candidate governance units in the conflict governance combinations that cause the adverse state of the associated river segment to continue to propagate and cannot be offset by other counterfactual response paths, and generating constraint-satisfied structures; connecting candidate governance units in the constraint-satisfied structures that jointly reduce the contradictory action pairs and maintain spatial connection relationships, and generating collaborative governance combinations; retaining the collaborative governance combinations and their corresponding action transmission order and construction dependency order, and generating a feasible implementation diagram.

[0014] Optionally, the step of parsing the action transmission order and construction dependency order of the candidate treatment units in the feasible implementation diagram to generate a river treatment plan includes: reading the action transmission order in the feasible implementation diagram to generate a treatment unit action sequence; adjusting the candidate treatment units in the treatment unit action sequence that cannot be directly implemented according to the action sequence based on the construction dependency order to generate a constructable treatment sequence; extracting the action separation position, change method, and connection object corresponding to each candidate treatment unit in the constructable treatment sequence to generate treatment unit implementation information; and combining the treatment unit implementation information with the constructable treatment sequence to generate a river treatment plan.

[0015] Optionally, the method further includes: acquiring river monitoring data generated after the implementation of the river management plan, and generating a river state sequence after implementation; comparing the river state sequence after implementation with the river state change units corresponding to the counterfactual response path, and generating path response differences; locating the causal transmission edge where the path response difference first appears in the river evolution map, and generating a deviation transmission position; correcting the transmission direction or degree of the corresponding causal transmission edge according to the deviation transmission position, and generating a corrected river evolution map; and using the corrected river evolution map to correct the river management plan.

[0016] Based on the same inventive concept, this invention also provides a river management scheme generation system based on river morphology evolution maps. The system includes: a state sequence generation module, used to acquire multi-period river observation data of the target river, perform location registration and time sorting on the multi-period river observation data, and generate a river morphology state sequence; an evolution map generation module, used to identify river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connect the river morphology change units with sequential response relationships, and generate a river morphology evolution map; a contradiction location module, used to acquire the management objectives and engineering constraints of the target river, map the management objectives and engineering constraints to the river morphology evolution map for contradiction extraction, and generate a contradiction sub-graph; and a counterfactual processing module, used for... Candidate governance units are constructed based on the transmission direction or degree of river morphology change units in the contradictory subgraph. These candidate governance units are then embedded one by one into the contradictory subgraph, and the changes in the graph before and after embedding are compared to generate counterfactual response paths. A combined response generation module is used to combine the counterfactual response paths with spatial and temporal connections, identify the convergence, cancellation, and amplification of effects after combination, and generate a combined response structure. A scheme diagram generation module is used to use the engineering constraints to trim the candidate governance units and their corresponding connections in the combined response structure that cannot be implemented together, generating a feasible scheme diagram. A scheme output module is used to parse the effect transmission order and construction dependency order of the candidate governance units in the feasible scheme diagram to generate a river management scheme.

[0017] Compared with the prior art, the present invention has the following advantages:

[0018] This invention converts multi-period river observation data into a river morphology evolution map with causal transmission relationships. It unifies the originally scattered observation records of shorelines, channels, sandbars, and main streams onto the same spatial coordinates and time axis and organizes them into continuous river segments. This allows the evolution of river morphology over time to be expressed in a calculable and structured form, enabling those skilled in the art to systematically model and quantitatively analyze complex river morphology changes.

[0019] This invention maps governance objectives and engineering constraints on a graph and extracts contradictory subgraphs by tracing causal transmission edges in reverse. It centrally presents the conflicts between governance objectives and natural evolution, as well as between the current river segment and related river segments. This avoids repeated calculations for the entire river channel and allows the contradictions between governance needs and river evolution patterns to be located, reducing the computational scope and manual judgment workload required for scheme design.

[0020] This invention utilizes a counterfactual reasoning mechanism that embeds candidate governance units one by one into a contradictory subgraph and compares the changes in the graph before and after embedding. This allows for the prediction of the response effect of each governance action after transmission along the river section before the actual implementation of the project. Furthermore, it identifies the convergence, cancellation, and amplification relationships when multiple actions are implemented simultaneously, enabling the synergy and conflict between governance actions to be quantitatively evaluated during the design phase, thereby improving the scientific nature and feasibility of the governance scheme.

[0021] This invention utilizes the comparison between post-implementation monitoring data and counterfactual response paths to locate the deviation transmission position and correct the causal transmission edge, enabling the river morphology evolution map to be continuously corrected based on the actual engineering response. The scheme is then iteratively updated, allowing the governance scheme to continuously approach the true river morphology response in long-term operation, thereby enhancing the method's adaptability to different river conditions. Attached Figure Description

[0022] Figure 1 This is a flowchart illustrating the method for generating river management schemes based on river morphology evolution maps according to an embodiment of the present invention.

[0023] Figure 2 This is a schematic diagram of continuous river segment division according to an embodiment of the present invention.

[0024] Figure 3 This is a schematic diagram of the construction of the river state sequence according to an embodiment of the present invention.

[0025] Figure 4 This is a schematic diagram of the river morphology evolution map structure according to an embodiment of the present invention.

[0026] Figure 5 This is a schematic diagram of contradictory subgraphs and counterfactual deductions in an embodiment of the present invention.

[0027] Figure 6 This is a schematic diagram of the combined response structure according to an embodiment of the present invention.

[0028] Figure 7 This is a schematic diagram of an embodiment of the present invention.

[0029] Figure 8 This is a schematic diagram of the structure of the river management scheme generation system based on river morphology evolution map according to an embodiment of the present invention. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0031] Reference Figure 1 One embodiment of the present invention proposes a method for generating river management schemes based on river morphology evolution maps. By constructing a river morphology evolution map with causal transmission edges from multi-period river observation data, mapping management objectives and engineering constraints on the map to extract contradictory subgraphs, embedding candidate management units for counterfactual inference and combination evaluation, and then correcting the map based on post-implementation monitoring data, the method can uniformly express the causal transmission relationship of the river morphology along the river section, quantitatively evaluate the synergistic and conflicting effects of various management actions, and enable the management scheme to be continuously iterated and optimized with actual response.

[0032] The method described in this embodiment specifically includes:

[0033] S1. Obtain multi-period river observation data of the target river channel, perform location registration and time sorting on the multi-period river observation data, and generate a river state sequence.

[0034] In one embodiment of the present invention, step S1 includes the following steps:

[0035] Acquire multi-period river observation data of the target river channel, and extract shoreline location data, channel elevation data, sandbar boundary data and main channel location data from the multi-period river observation data to generate phased river channel status data;

[0036] Correct the positional deviations and observation scale differences between the phased river state data to generate river state data of the same scale;

[0037] Divide the continuous river segments in the same-scale river state data and establish a time index for each continuous river segment to generate segmented time-series state data;

[0038] Connect adjacent time-phase states belonging to the same continuous river segment in the segmented time-series state data to generate a river morphology state sequence.

[0039] Specifically, river observation data comes from remote sensing images, underwater topographic survey records, and hydrological cross-section survey records collected for the same river section in different years or at different hydrological periods. Each collection constitutes the raw material for phased river state data. The river state sequence is a continuous set of states obtained by arranging multiple observation results on a unified spatial coordinate and a unified time axis, used to reflect the changes in river morphology over time.

[0040] Four types of data were extracted from each observation period: shoreline location data, channel elevation data, sandbar boundary data, and main stream location data. Shoreline location data refers to the set of planar coordinate points at the water-land interface, reflecting the planar orientation of the river channel boundaries. Channel elevation data refers to the vertical height of each measuring point at the bottom of the riverbed relative to a unified reference surface, reflecting the depth distribution of the cross-section. Sandbar boundary data refers to the set of coordinate points on the outer edges of sandbars or beaches exposed above the water surface in the middle of the river channel, reflecting the location and shape of water-blocking bodies within the channel. Main stream location data refers to the line connecting the centers of the areas with the highest flow velocity in the cross-section, reflecting the planar path of the main flow force. These four types of data collected in the same period were combined to form phased river channel state data, denoted as the [number missing]. The phased river channel status data are as follows:

[0041] ,

[0042] in Indicates the first Phased river channel status data, Indicates the first Coastline location data for the period Indicates the first Riverbed elevation data for the period Indicates the first Qizhou Beach boundary data, Indicates the first mainstream location data, This indicates the collection sequence number, with a value range of [value range missing]. to , This indicates the total number of observation periods.

[0043] Different acquisition periods will result in deviations in planar position and observation scale due to equipment positioning errors, coordinate system differences, and imaging resolution variations. Correction for these positional deviations employs rigid body transformations involving translation and rotation to unify the coordinates of each period to the reference period coordinate system. For the [missing information - likely a specific data point or data point]... Any coordinate point in the middle Corrected coordinates satisfy:

[0044] ,

[0045] in Indicates the first Issue No. A planar vector of the original coordinate points. Represents the corrected planar vector. Indicates the first The rotation matrix of the period relative to the reference period is a 2x2 orthogonal matrix. Indicates the first The translation vector of the period relative to the reference period. This represents the coordinate point number. The rotation matrix and translation vector are obtained by matching stable control points shared by both periods, minimizing the sum of squared positional errors before and after control point correction. The error objective is:

[0046] ,

[0047] in Indicates the first Control point matching error during the period Indicates the first Issue No. Coordinates of control points This indicates the coordinates of the control points corresponding to the reference period. Indicates the first The number of control points participating in the matching process. This represents the Euclidean length of the vector. The dimension of this error is the square of the length, consistent on both sides. Stable ground control points are selected from bridge pier foundations, fixed hydrological markers, or bedrock outcrops that remain unchanged across time periods. The number is determined based on identifiable ground features on site, with no fewer than [number missing] points selected for each period in actual measurements. Several control points are used to ensure stable solution.

[0048] Observation scale differences refer to the inconsistency in spatial resolution between different periods; for example, the difference in spatial resolution between each pixel in a certain period. Another phase corresponds to Mi Meters. All periods are uniformly resampled to the baseline scale, with a resampling ratio of:

[0049] ,

[0050] in Indicates the first Scaling ratio, Indicates the first The ground length corresponding to a single pixel in the period. This represents the ground length corresponding to a single pixel in the baseline period. Both dimensions are length, while the scale itself is dimensionless. After location and scale correction, elements from each period are combined to form river channel state data at the same scale.

[0051] River channel condition data at the same scale is still a holistic record of the entire river and needs to be spatially divided into continuous river segments for segment-by-segment analysis. For example... Figure 2 As shown, equally spaced cross-sections are arranged along the direction of the main channel location data. The river sections between adjacent cross-sections form a continuous river segment. A continuous river segment is denoted as A time index is added to each period of the state of each continuous river segment. The time index is taken from the actual duration of the observation in that period, resulting in segmented time-series state data:

[0052] ,

[0053] in Indicates the first River section Segmented time-series state data for the period, Indicates the first The spatial range of a continuous river segment Indicates the first The time coordinates corresponding to the period of observation Indicates falling into the first River channel condition data at the same scale within the river section. Indicates the river section number.

[0054] Reference Figure 3 By connecting the states of the same continuous river segment in adjacent time phases sequentially from earliest to latest according to the time coordinates, the state evolution chain of that river segment is obtained. The sum of the chains of all river segments constitutes the river morphology state sequence:

[0055] ,

[0056] in Represents a sequence of river state conditions. This represents a pair of adjacent states in the same river segment, subject to time constraints. Ensure the connection direction is in the positive time direction.

[0057] For example, data was collected from a certain section of the river. Periodic data, the ground length of a single pixel is respectively rice, rice, rice, Rice, take the first The period is the base period. Scaling ratio according to Calculate, put the first The length of each pixel in the periodic image is calculated according to... Double the conversion, making Rice pixels are unified to Meter scale. Position correction during the first... Expected to obtain There are several bridge pier control points, among which the coordinates of one control point before correction are: The coordinates of the corresponding control points in the reference period are: Suppose that the rotation matrix is ​​close to the identity matrix and the translation vector is obtained. Then the corrected coordinates are as follows have to The error term at this point coincides with the base period. Square meters. After the stationing was laid out along the main stream, this section of the river was designated as... Section, No. Duan Di Period and the The time coordinate of the period is , ,satisfy Therefore, and By connecting them sequentially, and connecting all adjacent periods of the river segment in this way, a sequence of river morphology states for that segment can be obtained. .

[0058] S2. Identify river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connect the river morphology change units with sequential response relationship, and generate a river morphology evolution map.

[0059] In one embodiment of the present invention, step S2 includes the following steps:

[0060] By comparing the spatial differences and directions of change of adjacent time phases in the river morphology state sequence, river morphology change segments are generated.

[0061] Merge the spatially connected, consistent-direction, and continuous-stage river morphology change segments to generate river morphology change units.

[0062] The occurrence order and response direction of the river morphology change units in the direction of water flow connectivity are tracked to generate candidate transmission relationships;

[0063] Eliminate candidate transitive relationships that lack continuous temporal responses or whose evolution directions are mutually exclusive, and generate causal transitive edges;

[0064] Connect the river morphology change unit with the causal transmission edge to generate a river morphology evolution map.

[0065] Specifically, a river morphology change unit refers to a region of river morphology change that is spatially connected and maintains a consistent direction of evolution. A river morphology evolution map is a directed structure obtained by connecting these change units according to their causal relationships, and is used to express the path of change transmission from the upstream section to the downstream section.

[0066] By comparing the spatial differences between adjacent time phases in the river morphology sequence, the changes in four elements—shoreline, channel, sandbars, and main stream—are calculated point by point. For the first... River section Period to the Period in cross section The change at the location is taken as follows:

[0067] ,

[0068] in Indicates the first The river section in Period to the Expected at cross-section The annual rate of change at that location Indicates the first River section Period in cross section The state value at that point, This represents the time interval between two periods. The dimension of the change in state value is length; dividing by time yields the dimension of length per time, representing the rate. The direction of change is determined by the sign of the change: a positive value indicates shoreline expansion or riverbed silting, while a negative value indicates shoreline retreat or riverbed erosion. Spatially adjacent cross-sections with the same sign of the rate of change are grouped together to obtain a segment of river morphology change.

[0069] When merging segments of river course changes, the segments must be spatially connected end-to-end, have the same direction of change, and be sequentially connected in terms of temporal duration. Assume that two segments satisfy the condition of spatial connection. If the conditions for interval connection, consistent direction, identical signs, and phase connection are met, then the intervals are merged into the same river morphology change unit. Its spatial extent is taken as the union of two spatial extents, and its duration is taken as the union of two time intervals. The intensity of a river morphology change unit is measured by the spatial average of the rate of change within its range.

[0070] ,

[0071] in Indicates the first The average intensity of change in each river morphology unit This represents the set of cross sections covered by the element. Indicates the number of cross sections within the set. Indicates the number of the changed unit. Dimensions and Consistent, with length per time.

[0072] The order of occurrence of each river morphology change unit along the direction of water flow connectivity is tracked, determined by the main current location data from upstream to downstream. If the starting point of the duration of an upstream unit is earlier than that of a downstream unit, and the two are spatially connected along the main current, then a candidate transmission relationship is considered to exist between them. The degree of consistency in the transmission direction is determined by the product of the directional signs of the change intensities of the two units.

[0073] ,

[0074] in Indicates upstream unit With downstream units Direction sign product, Indicates the sign function, , Representing units respectively , The average intensity of change, Dimensionless.

[0075] Candidate transit relationships were eliminated if they lacked continuous temporal response or exhibited mutually exclusive evolutionary directions. Lack of continuous temporal response refers to time gaps and lack of connection between the continuous stages of upstream and downstream units; mutually exclusive evolutionary directions refer to combinations that cannot be physically explained by the same kinetic chain, such as upstream scouring corresponding to downstream synchronous sedimentation. This contradicts the actual direction of sediment transport. The remaining relationships form causal transitive edges, and the degree of transitivity of the edges is taken as:

[0076] ,

[0077] in Indicates from unit To unit Causal transitive edge weights This indicates the distance between two units along the main flow. This represents the spacing attenuation coefficient, whose value is determined based on the attenuation pattern along the path observed in historical observations. To make the weights dimensionless, and Divide to eliminate the dimension of length per time, in the exponential term Take the reciprocal of the length and... The product is dimensionless, and the entire weight is dimensionless. (See reference...) Figure 4 By treating all river morphology change units as nodes and all causal transitive edges as directed connections, a river morphology evolution map is obtained:

[0078] ,

[0079] in A diagram showing the evolution of river course. Represents a set of river morphology change units. This represents a set of weighted causal transitive edges.

[0080] For example, in a certain section of the river... Period to the Period in adjacent The changes in shoreline at each cross-section are as follows: rice, rice, rice, rice, meters, time interval Year. Before. The rate of change of each cross section is as follows , , Meters are calculated annually, and fragments with the same symbol are grouped into collapse segments; later Cross-sectional rate , Each year, the data is categorized as an outward expansion segment. The average change intensity of the collapse segment is calculated according to... Calculated per meter per year. Assume the strength of the other downstream collapse unit. meters per year, product of direction signs The directions are consistent, and the two durations are connected sequentially, thus retaining them as causal propagation edges. The distance along the main current is taken. meters, attenuation coefficient For each meter, the edge weight is calculated as follows: The causal transit edge is calculated and combined with all units and edges to form a river morphology evolution map. .

[0081] S3. Obtain the governance objectives and engineering constraints of the target river channel, map the governance objectives and engineering constraints to the river morphology evolution map to extract contradictions, and generate a contradiction sub-map;

[0082] In one embodiment of the present invention, step S3 includes the following steps:

[0083] The governance objective is transformed into the target evolution direction of the river morphology change unit, and a target state marker is generated;

[0084] The engineering constraints are converted into constraint markers that limit the range and manner of change of river morphology units, and a map constraint boundary is generated.

[0085] Using the target state marker, the causal propagation edges are traced backward from the river morphology evolution map to generate the target action chain;

[0086] Using the map constraint boundary, the opposite evolution direction formed when the target action chain continues to propagate to the associated river segment is identified, generating contradictory action pairs;

[0087] Extract the river morphology change units and causal propagation edges that connect the contradictory action pairs and do not contain irrelevant branches to generate a contradictory subgraph.

[0088] Specifically, refer to Figure 5 The governance objective refers to the desired river state for a specific river section, such as curbing shoreline erosion or maintaining the stability of the main channel. Engineering constraints refer to the limitations imposed on governance actions by construction conditions, such as prohibiting the occupation of waterways, restricting construction during the flood season, or limiting the scope of engineering projects. The conflict sub-graph is a local structure extracted from the river evolution map that centrally reflects the conflict between the governance objective and natural evolution or with related river sections.

[0089] The governance objectives are transformed into the target evolution direction of river morphology change units, resulting in target state labels. For the units involved in the objectives... The expected average change intensity is to tend towards a set threshold, and the target state is marked as follows:

[0090] ,

[0091] in Representation unit The target evolution direction marker, Representation unit Current average change intensity This indicates the desired intensity of the governance objectives. Indicates the sign function, Dimensionless, positive values ​​indicate that adjustments need to be made in the direction of siltation or outward expansion, while negative values ​​indicate that adjustments need to be made in the direction of scouring or collapse.

[0092] The engineering constraints are converted into constraint markers that limit the range and manner of change of river morphology elements, resulting in the atlas constraint boundaries. The constraint markers are represented by the upper limit of the allowed intensity adjustment for each element.

[0093] ,

[0094] in Representation unit Allowable range of intensity values Represents the element given by engineering constraints. The intensity is adjustable, and its dimension is the same as that of intensity, which is length per time. The allowable engineering volume and the river section's carrying capacity are used for calibration. The sum of the allowable intervals of all units constitutes the map constraint boundary.

[0095] By using target state markers to trace upstream along the causal transmission edges in the river morphology evolution map, the chain of upstream changing units that can influence the target unit is identified, thus obtaining the target action chain. The reverse tracing starts from the target unit and proceeds upstream along the incoming edges, accumulating the degree of transmission:

[0096] ,

[0097] in Indicates from source unit To the target unit The cumulative transmission degree of the chain, Indicates from arrive The set of causal transit edges that pass through sequentially along the transit path. The corresponding edge weights are dimensionless, and the larger the value, the stronger the influence of the source cell on the target cell.

[0098] By using map-constrained boundaries to identify the opposing evolutionary directions formed as the target action chain continues to propagate to the associated river segment, contradictory action pairs are obtained. When the target requires a unit to adjust towards sedimentation, but the same action chain, when propagated to a downstream associated unit, causes that downstream unit to change towards scour beyond its allowable range due to physical continuity, the two units constitute a contradictory action pair. The determination criteria are:

[0099] ,

[0100] in Representing associated units Evolutionary direction markings resulting from action chain propagation. Represents the unit after conduction The predicted intensity Representation unit The allowable intensity range is defined as follows: a product of directions less than zero indicates opposite directions; a predicted intensity falling outside the allowable range indicates an out-of-bounds violation. Conditions are met. and They constitute a pair of contradictory forces.

[0101] Extract the river morphology change unit connecting the contradictory action pairs and the causal transmission edge, and remove the side branches unrelated to the contradiction to obtain the contradictory subgraph:

[0102] ,

[0103] in Represents a contradictory subgraph. Indicates the connection target unit Associated Units The transmission path.

[0104] For example, a governance objective in one area is to eliminate the collapsed units. The intensity from the current Rice increased to the target every year Meters per year, target markers according to The calculation indicates that adjustments are needed in the direction of siltation. Downstream related units. Current intensity Rice per year, permissible range meters per year, allowable intervals according to Calculation. The action chain transmits upstream sedimentation adjustments to downstream areas for prediction. Each year, this value does not fall below [a certain amount]. Within the interval, simultaneously Directional product according to Calculate if both conditions are met simultaneously. and This constitutes a contradictory pair. The weights along each edge of the chain are... and Cumulative transmission degree according to Calculations, and extraction accordingly , By passing edges and removing side branches, we obtain a contradictory subgraph. .

[0105] S4. Based on the transmission direction or degree of the river morphology change unit in the contradictory subgraph, construct candidate governance units, embed the candidate governance units one by one into the contradictory subgraph, compare the changes in the graph before and after embedding, and generate counterfactual response paths.

[0106] In one embodiment of the present invention, step S4 includes the following steps:

[0107] The spatial location or evolution stage of the contradictory action pair can be separated in the contradictory subgraph, generating the action separation location;

[0108] Matching governance actions that can change the state of water flow guidance, shoreline constraints, channel flow, or beach evolution at the separation point of action, and generating candidate governance units;

[0109] The candidate governance unit is embedded into the action separation position and the causal transmission edge connected to the corresponding action separation position is changed to generate an intervention graph branch;

[0110] By comparing the transmission direction and degree of the corresponding nodes in the intervention graph branch and the contradictory subgraph, a response difference chain is generated;

[0111] Extract the continuous portion of the response difference chain that extends from the candidate governance unit to the target improved state or the unfavorable state of the associated river segment, and generate a counterfactual response path.

[0112] Specifically, refer to Figure 5 Candidate governance units refer to engineering actions taken at specific locations to change river morphology, such as bank protection, groynes, dredging, or sandbar protection; counterfactual response paths refer to the transmission routes of the difference in river morphology after the deployment of a candidate governance unit compared to when it is not deployed, used to compare the difference in effects of having and not having governance.

[0113] Identifying the spatial location or evolutionary stage of contradictory action pairs within a contradictory subgraph reveals the separation point. The separation point is the key node that severs the causal transmission between contradictory action pairs; its determination is measured by the decrease in the cumulative degree of transmission between the contradictory action pairs after removing a node.

[0114] ,

[0115] in Indicates the removal of a node. The subsequent transmission reduction rate This indicates the cumulative degree of transmission between contradictory pairs before removal. Indicates the removal of a node. The cumulative degree of transmission afterward, Dimensionless, the closer the value is to This indicates that the node is more critical. The node position where the reduction rate exceeds the set threshold is the effective separation position.

[0116] Candidate treatment units are obtained by matching treatment actions that can alter the flow guidance, shoreline constraints, channel flow, or sandbar evolution at the point of separation. The change in the intensity of each treatment action on the unit is denoted as the action effect.

[0117] ,

[0118] in Indicates the first The change in unit strength caused by a governance action, measured in units of length per time. Indicates the first The scale of engineering deployment for this type of action, This represents the intensity change coefficient resulting from a unit deployment scale. Based on historical deployments of similar projects and measured river response, the values ​​are selected to make... and The multiplication result in a dimension of length per time. Bank protection actions primarily reduce the intensity of erosion, groynes diverting currents primarily change the position of the main stream and adjust the flow of branch channels, dredging actions primarily increase the flow area of ​​the river channel, and sandbar protection actions primarily stabilize the boundaries of sandbars.

[0119] By embedding candidate governance units into the action separation location and altering the causal transitive edges connected to that location, an intervention graph branch is obtained. After embedding, the intensity of the affected unit is the sum of its original intensity and the action effect.

[0120] ,

[0121] in Indicates the unit after embedding governance actions The strength, This represents the strength before embedding. After the strength of the affected unit changes, the weights of its downstream causal propagation edges are updated according to the weight definition in step S2, forming an intervention graph branch.

[0122] By comparing the transmission direction and degree of corresponding nodes in the intervention map branches and contradictory subgraphs, the response difference chain is obtained. The intensity difference before and after intervention is calculated for each affected unit.

[0123] ,

[0124] in Representation unit The intensity difference caused by the governance action is equal to the action effect. Units with non-zero intensity differences that are continuous along the transmission direction are connected in series to form a response difference chain.

[0125] By extracting the continuous portion of the response difference chain that extends from the candidate governance unit to the target improved state or to the unfavorable state of the associated river segment, the counterfactual response path is obtained:

[0126] ,

[0127] in Indicates the first Counterfactual response paths corresponding to candidate governance units Indicates the unit where the governance action takes place. Indicates the target improvement unit. This indicates an unfavorable unit in the associated river segment. The path record shows the complete process of an action propagating from the placement location along the causal propagation edge to the target unit or unfavorable unit.

[0128] For example, the intermediate unit located between contradictory action pairs in the contradictory subgraph The cumulative degree of transmission before and after removal is , Reduction rate according to Calculation, exceeding the threshold , The location is defined as the action separation point. At this location, the revetment action is matched, and the intensity change coefficient is taken. Rice per meter per year, scale of deployment Linear meter, motion effect according to Calculated per meter per year. Strength of the collapsed unit before embedding. Rice per year, after embedding according to Calculated per meter per year, strength difference according to The annual yield is [amount missing]. The downstream connected units of this unit have the same and continuous strength difference, extending from the revetment unit all the way to the target unit. By cascading consecutive units, a counterfactual response path is obtained. .

[0129] S5. Combine the counterfactual response paths that have spatial and temporal connections, identify the convergence, cancellation, and amplification of effects that occur after combination, and generate a combined response structure.

[0130] In one embodiment of the present invention, step S5 includes the following steps:

[0131] Arrange the counterfactual response paths in the order in which they arrive at the same river morphology change unit to generate a path temporal relationship;

[0132] By combining the path temporal relationships to compare the evolution directions reaching the same river morphology change unit, path interaction relationships are generated;

[0133] Identify path combinations with the same evolution direction and increasing transmission degree, path combinations with opposite evolution direction and decreasing transmission degree, and path combinations that end transmission at the same node in the path interaction relationship, and generate amplification relationship, cancellation relationship and convergence relationship;

[0134] Connect candidate governance units that have the amplification relationship, the cancellation relationship, or the convergence relationship, and retain the corresponding path timing relationship to generate a combined response structure.

[0135] Specifically, refer to Figure 6 A combined response structure refers to the overall action structure formed when multiple counterfactual response paths are implemented simultaneously, interconnected, and mutually influential. It is used to express the synergistic or antagonistic relationship between governance actions. Action convergence refers to multiple paths ending their propagation within the same unit; action cancellation refers to paths in opposite directions causing a decrease in the overall intensity change within the same unit; and action amplification refers to paths in the same direction causing an increase in the overall intensity change within the same unit.

[0136] The timing relationship of multiple counterfactual response paths arriving at the same river morphology unit is obtained by arranging the order in which they arrive. The arrival time is determined by adding the propagation delay along the path to the completion time of the action.

[0137] ,

[0138] in Representing a path The function reaches the unit At that moment, Representing a path The corresponding completion time of the remediation actions. Indicates the position of the action to the unit The spacing along the main stream, Indicates the speed at which changes in river morphology are transmitted along the course of the river. The dimension is length divided by length per time, which is time, consistent with the dimension of the moment construction is completed. Transmission velocity. Based on historical observations of the distance and time of the advancing front, the path time sequence is obtained by sorting the arrival times from earliest to latest.

[0139] By comparing the evolution directions of paths reaching the same river morphology unit based on their temporal relationships, the path interaction relationships are obtained. For two paths reaching the same unit, the directional relationship is determined by the sign of their respective intensity changes, and the combined intensity change after the interaction is superimposed is:

[0140] ,

[0141] in Representation unit The intensity change resulting from the synthesis of multiple paths. Indicates arrival unit The set of paths Representing a path In unit The change in intensity at a given point is measured in units of length per time.

[0142] Path combinations with the same evolution direction and increased transmission after synthesis are identified as amplification relationships; path combinations with opposite evolution directions and decreased transmission after synthesis are identified as cancellation relationships; and path combinations ending in the same unit are identified as convergence relationships. These three types of relationships are determined by comparing the synthesis strength with the individual path strength.

[0143] ,

[0144] in Representation unit The superposition coefficient at the point is dimensionless. Greater than It is judged to be an amplification relationship. Less than If multiple paths terminate in the same unit, they are considered to be in a merging relationship.

[0145] Connect candidate governance units that have amplification, cancellation, or convergence relationships, and preserve the corresponding path temporal relationships to obtain the combined response structure:

[0146] ,

[0147] in Represents the combined response structure. Represents the set of candidate governance units. This represents the set of governance units connected by the superposition coefficients and arrival times.

[0148] For example, path Corresponding to the completion time of the revetment construction Spacing to target unit meters, conduction speed Rice per year, arrival time according to Calculation. Path The arrival time is calculated in the same way as the groyne's diversion action. Both paths arrive at the same target cell sequentially. The intensity changes of both paths within that cell are positive, respectively. , Rice per year, synthetic intensity according to Calculated annually in meters, with the superposition factor based on... Calculation, greater than This is determined to be an amplification relationship. Based on this, the revetment unit and the groyne unit are connected, and the arrival time sequence is preserved. and The combined response structure is obtained. .

[0149] S6. Using the engineering constraints, trim the candidate governance units and corresponding connection relationships that cannot be implemented together in the combined response structure to generate an implementation scheme diagram;

[0150] In one embodiment of the present invention, step S6 includes the following steps:

[0151] The engineering constraints are mapped to the candidate governance units and path timing relationships in the combined response structure to generate governance unit constraint relationships;

[0152] Determine whether the constraint relationship of the governance unit has spatial occupancy conflict, implementation timing conflict, or crosses the graph constraint boundary after the effect is amplified, and generate a conflict governance combination;

[0153] Remove candidate governance units from the conflict governance combination that cause the unfavorable state of the associated river segment to continue to propagate and cannot be offset by other counterfactual response paths, and generate a constraint satisfaction structure.

[0154] Connect the candidate governance units in the constraint-satisfying structure that jointly reduce the contradictory effects and maintain spatial connectivity to generate a collaborative governance combination;

[0155] Retaining the aforementioned collaborative governance combination and its corresponding action transmission sequence and construction dependency sequence, a feasible implementation diagram is generated.

[0156] Specifically, refer to Figure 7 The feasible implementation diagram refers to the structure that removes the governance actions that cannot be carried out simultaneously, retains the governance actions that can synergistically reduce conflicts, and their connection sequence, and is used to provide a feasible governance combination.

[0157] By mapping engineering constraints to candidate governance units and path sequence relationships in the combined response structure, the governance unit constraint relationships are obtained. The constraint mapping appends the spatial occupancy range, construction time window, and deployment scale upper limit of each governance action to the corresponding unit.

[0158] ,

[0159] in Indicates governance actions The constraint relationship, This indicates the spatial range occupied by the action along the path. Indicates the time window during which construction is permitted. This indicates the upper limit of the deployment scale.

[0160] To determine whether there are spatial occupancy conflicts, implementation timing conflicts, or situations where the effect amplifies and exceeds the map constraint boundary, a conflict governance combination is obtained. Spatial occupancy conflicts refer to the overlap of spatial intervals between two actions, and the criterion is:

[0161] ,

[0162] Implementation sequence conflict refers to two mutually exclusive construction time windows that are required to be sequentially connected; effect amplification exceeding the limit refers to the unit strength falling outside the allowable range after the amplification relationship is synthesized. The criteria are as follows:

[0163] ,

[0164] The meanings of the symbols are consistent with those in the preceding text. Action combinations that satisfy any one of the criteria are categorized into conflict management combinations.

[0165] Remove candidate governance units from the conflict governance combination that cause the adverse state of the associated river segment to continue to propagate, and whose adverse propagation cannot be offset by other counterfactual response paths, to obtain a constraint-satisfied structure. The offsetting capability is determined by the change in residual strength at the adverse unit.

[0166] ,

[0167] in Indicates unfavorable unit The remaining intensity change after all preserved paths are superimposed, if it still makes If the allowed range is exceeded and there is no opposite path to offset it, then the action that caused the transmission is removed.

[0168] Candidate governance units that satisfy the connection constraints to jointly reduce contradictory effects in the structure and maintain spatial connectivity are obtained as collaborative governance combinations. Then, by retaining the action transmission order and construction dependency order corresponding to the collaborative governance combinations, a feasible implementation diagram is obtained:

[0169] ,

[0170] in This represents a feasible implementation diagram. This represents the set of candidate governance units to be retained. Indicates the arrival time of the action. Indicates the construction dependency sequence. Record the sequential relationship where the completion of one action depends on the completion of another action.

[0171] For example, the operational space range of the revetment is: Meters, the dredging operation space range is Rice, intersection Neither is empty, and the two constitute a space occupation conflict. Each year, the intensity of the amplified relationship changes. Rice per year, accumulated according to Calculated in meters per year, Falling In addition, it constitutes an amplified and excessive action. Therefore, dredging operations and some revetment deployments are included in the conflict management framework, removing any related units. Dredging operations that are detrimental to transmission and have no counter-pathway for offsetting should be retained, along with bank protection operations with reduced deployment scale. Returning to the interval, connecting it with the groyne diversion action forms a collaborative governance combination. Retaining the arrival time sequence and construction dependency sequence, a feasible implementation diagram is obtained. .

[0172] S7. Analyze the action transmission order and construction dependency order of the candidate treatment units in the feasible implementation diagram to generate a river treatment plan.

[0173] In one embodiment of the present invention, step S7 includes the following steps:

[0174] Read the action transmission sequence in the proposed implementation diagram to generate the action sequence of the governance unit;

[0175] Based on the construction dependency sequence, candidate treatment units that cannot be directly implemented in the order of action of the treatment units are adjusted to generate a constructable treatment sequence.

[0176] Extract the separation location, change method, and connection object of each candidate treatment unit in the constructable treatment sequence to generate treatment unit implementation information;

[0177] By combining the implementation information of the governance units with the feasible governance sequence, a river governance plan is generated.

[0178] Specifically, a river management plan refers to a complete management arrangement arranged in the order of feasibility and including implementation information for each management action, for use on-site.

[0179] By reading the action propagation order in the feasible solution diagram, and arranging the candidate governance units in ascending order of action arrival time, we obtain the governance unit action sequence:

[0180] ,

[0181] in Indicates the sequence of actions of governance units. This indicates that the order of arrival of the function is ranked as follows: Candidate governance units, Indicates the time when its effect arrives. This indicates the total number of governance actions that are retained.

[0182] By adjusting the sequence of remediation units in the action sequence that cannot be directly implemented according to the order of actions, a constructable remediation sequence is obtained. When an action should be implemented first according to the action sequence, but its construction depends on the completion of another action, the construction order is exchanged according to the dependency relationship. The adjustment rule is as follows:

[0183] ,

[0184] in Indicates action Construction-dependent actions Complete it first. This indicates the order of construction. After rearranging according to the dependencies, a workable treatment sequence is obtained. .

[0185] Extract the action separation location, change method, and connection object corresponding to each candidate treatment unit in the constructable treatment sequence to obtain the treatment unit implementation information:

[0186] ,

[0187] in Indicates governance actions Implementation information, The coordinates along the path indicate the location of the separation of actions; the method of representation can be changed to indicate the specific actions taken in bank protection, groynes diversion, dredging, or sandbar protection. Indicates the scale of deployment. This indicates the downstream unit in the graph where the action is connected.

[0188] By combining the implementation information of the treatment units with the feasible treatment sequence, a river treatment plan is obtained:

[0189] ,

[0190] in Indicates the river management plan, Indicates the sequence of feasible treatments. This represents the collection of information on all governance actions implemented.

[0191] For example, the implementation diagram can retain two actions: the revetment and the groynes diverting the current, with the arrival times of their effects being respectively... and Arranged in order of function bank protection, groynes and diversion channels The construction of groynes requires the completion of the revetment as a prerequisite, and is dependent on this condition. Therefore, the construction order remains that the revetment precedes the groynes. Maintain this order. Information on the implementation of revetment actions is obtained by taking the coordinates of the separation position. Rice, changing the method to bank protection, and the scale of deployment Linear meters, connecting downstream units The information on the implementation of the groyne action was extracted using the same method, and the information on the implementation of the two actions was combined with... By combining the results, a river management plan can be obtained. .

[0192] S8. Obtain river monitoring data generated after the implementation of the river management plan, and generate a river state sequence after implementation;

[0193] Compare the implemented river state sequence with the river state change units corresponding to the counterfactual response path to generate path response differences;

[0194] Locate the causal propagation edge where the path response difference first appears in the river morphology evolution map, and generate the deviation propagation location;

[0195] Based on the deviation transmission location, the transmission direction or degree of the corresponding causal transmission edge is corrected to generate a corrected river morphology evolution map.

[0196] The revised river management plan was revised using the revised river morphology evolution map.

[0197] Specifically, the post-implementation river state sequence refers to the state sequence generated by re-collecting observation data and following the same procedure as step S1 after the completion of the river management project according to the river management plan. It is used to verify the management effect and, in turn, correct the evolution map and management plan.

[0198] Acquire river monitoring data generated after the implementation of the river management plan, and generate a river morphology sequence after implementation by registering the data by location and sorting it by time. Its generation method is consistent with step S1, ensuring that it can be directly compared with the sequence before treatment.

[0199] By comparing the implemented river morphology state sequence with the river morphology change units corresponding to the counterfactual response path, the path response differences are obtained. For each unit covered by the path, the difference between the measured intensity and the counterfactual predicted intensity is calculated:

[0200] ,

[0201] in Representation unit Differences in path response, This represents the measured element strength after implementation. This represents the predicted cell strength after the governance action is embedded in step S4, with the dimension being length per time. The larger the absolute value of the difference, the more significant the deviation between the prediction and the actual value.

[0202] The location of the deviation propagation is determined by identifying the first causal propagation edge appearing in the river morphology evolution map to pinpoint the path response difference. Following the counterfactual response path downstream from the unit where the governance action was implemented, each edge is examined, and the first propagation edge whose absolute difference value exceeds a set threshold is identified.

[0203] ,

[0204] in The causal transitive edge represents the first occurrence of the deviation. This represents the threshold for determining the difference, measured in units of length per time, and is set based on the level of observation error. This represents the absolute value of the response difference between the edge endpoint cells.

[0205] Based on the location of the deviation propagation, the propagation direction or degree of the corresponding causal propagation edge is corrected to obtain the corrected river morphology evolution map. The correction adjusts the weight of this edge according to the measured difference ratio:

[0206] ,

[0207] in This represents the corrected edge weights. This indicates that the weights have been corrected. The numerator is the measured upstream-downstream intensity difference, and the denominator is the predicted upstream-downstream intensity difference. The ratio is dimensionless, ensuring that the corrected weights remain dimensionless. If the sign of the measured difference is opposite to the predicted sign, the transmission direction of that edge is simultaneously reversed. Replacing the corrected edge back with the map yields the corrected river morphology evolution map. .

[0208] By re-executing steps S3 to S7 using the revised river morphology evolution map, a combination of governance actions that better matches the actual river morphology response is obtained, thus completing the revision of the river management plan.

[0209] For example, the counterfactual prediction strength of a certain unit Meters per year, actual measured intensity after implementation Meters per year, path response differences according to Meters calculated annually, absolute value Exceeding the threshold The unit's ingress edge is defined as the location of deviation propagation, measured in meters per year. Let the predicted upstream-downstream strength difference at this edge be... The annual average intensity difference between upstream and downstream areas is [amount missing]. Rice per year, original edge weight The corrected weights are based on Calculate. Transfer the weights from... Updated to Replacing the map with the corrected river course evolution map yields the revised map. Then, by retracing the process of contradiction extraction and solution generation, a revised river management solution was obtained.

[0210] Based on the same inventive concept, such as Figure 8 As shown, the present invention also provides a river management scheme generation system based on river morphology evolution maps, the system comprising:

[0211] The state sequence generation module is used to acquire multi-period river observation data of the target river, perform location registration and time sorting on the multi-period river observation data, and generate a river state sequence.

[0212] The evolution map generation module is used to identify river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connect the river morphology change units with sequential response relationship, and generate a river morphology evolution map.

[0213] The conflict localization module is used to obtain the governance objectives and engineering constraints of the target river channel, map the governance objectives and engineering constraints to the river evolution map to extract conflicts, and generate a conflict sub-map;

[0214] The counterfact processing module is used to construct candidate governance units based on the transmission direction or degree of transmission of river morphology change units in the contradictory subgraph, embed the candidate governance units one by one into the contradictory subgraph and compare the changes in the graph before and after embedding, and generate counterfact response paths.

[0215] The combined response generation module is used to combine the counterfactual response paths that have spatial and temporal connections, identify the convergence, cancellation, and amplification of effects that occur after combination, and generate a combined response structure.

[0216] The scheme diagram generation module is used to use the engineering constraints to trim the candidate governance units and corresponding connection relationships that cannot be implemented together in the combined response structure, and generate a scheme diagram.

[0217] The scheme output module is used to parse the action transmission order and construction dependency order of the candidate treatment units in the feasible scheme diagram and generate a river treatment scheme.

[0218] All equivalent changes and modifications made in accordance with the teachings of this invention are still within the scope of this invention. Those skilled in the art will readily conceive of other embodiments of this invention upon considering the specification and the disclosure of practical truth. This application is intended to cover any variations, uses, or adaptations of this invention that follow the general principles of this invention and include common knowledge or conventional techniques in the art not described herein.

Claims

1. A method for generating river management schemes based on river morphology evolution maps, characterized in that, The method includes: Acquire multi-period river observation data of the target river channel, perform location registration and time sorting on the multi-period river observation data, and generate a river state sequence; Identify river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connect the river morphology change units with sequential response relationship, and generate a river morphology evolution map; Obtain the governance objectives and engineering constraints of the target river channel, map the governance objectives and engineering constraints to the river morphology evolution map to extract contradictions, and generate a contradiction sub-map; Candidate governance units are constructed based on the transmission direction or degree of river morphology change units in the contradictory subgraph. The candidate governance units are then embedded into the contradictory subgraph one by one, and the changes in the graph before and after embedding are compared to generate counterfactual response paths. The counterfactual response paths with spatial and temporal connections are combined to identify the convergence, cancellation, and amplification of effects that occur after combination, and to generate a combined response structure. By utilizing the engineering constraints, candidate governance units and their corresponding connections that cannot be implemented together in the combined response structure are trimmed to generate a feasible implementation diagram; The river management scheme is generated by analyzing the action transmission order and construction dependency order of the candidate management units in the feasible implementation diagram.

2. The method for generating river management schemes based on river morphology evolution maps according to claim 1, characterized in that, The process of acquiring multi-period river observation data of the target river, performing location registration and time sorting on the multi-period river observation data, and generating a river state sequence includes: Acquire multi-period river observation data of the target river channel, and extract shoreline location data, channel elevation data, sandbar boundary data and main channel location data from the multi-period river observation data to generate phased river channel status data; Correct the positional deviations and observation scale differences between the phased river state data to generate river state data of the same scale; Divide the continuous river segments in the same-scale river state data and establish a time index for each continuous river segment to generate segmented time-series state data. Connect adjacent time-phase states belonging to the same continuous river segment in the segmented time-series state data to generate a river morphology state sequence.

3. The method for generating river management schemes based on river morphology evolution maps according to claim 2, characterized in that, The process of identifying river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connecting the river morphology change units with sequential response relationships, and generating a river morphology evolution map includes: By comparing the spatial differences and directions of change of adjacent time phases in the river morphology state sequence, river morphology change segments are generated. Merge the spatially connected, consistent-direction, and continuous-stage river morphology change segments to generate river morphology change units. The occurrence order and response direction of the river morphology change units in the direction of water flow connectivity are tracked to generate candidate transmission relationships; Eliminate candidate transitive relationships that lack continuous temporal responses or whose evolution directions are mutually exclusive, and generate causal transitive edges; Connect the river morphology change unit with the causal transmission edge to generate a river morphology evolution map.

4. The method for generating river management schemes based on river morphology evolution maps according to claim 3, characterized in that, The process of obtaining the governance objectives and engineering constraints of the target river channel, mapping the governance objectives and engineering constraints to the river morphology evolution map for contradiction extraction, and generating a contradiction sub-map includes: The governance objective is transformed into the target evolution direction of the river morphology change unit, and a target state marker is generated; The engineering constraints are converted into constraint markers that limit the range and manner of change of river morphology units, and a map constraint boundary is generated. Using the target state marker, the causal propagation edges are traced backward from the river morphology evolution map to generate the target action chain; Using the map constraint boundary, the opposite evolution direction formed when the target action chain continues to propagate to the associated river segment is identified, generating contradictory action pairs; Extract the river morphology change units and causal propagation edges that connect the contradictory action pairs and do not contain irrelevant branches to generate a contradictory subgraph.

5. The method for generating river management schemes based on river morphology evolution maps according to claim 4, characterized in that, The process of constructing candidate governance units based on the transmission direction or degree of river morphology change units in the contradictory subgraph, embedding each candidate governance unit into the contradictory subgraph one by one, and comparing the changes in the graph before and after embedding to generate a counterfactual response path includes: The spatial location or evolution stage of the contradictory action pair can be separated in the contradictory subgraph, generating the action separation location; Matching governance actions that can change the state of water flow guidance, shoreline constraints, channel flow, or beach evolution at the separation point of action, and generating candidate governance units; The candidate governance unit is embedded into the action separation position and the causal transmission edge connected to the corresponding action separation position is changed to generate an intervention graph branch; By comparing the transmission direction and degree of the corresponding nodes in the intervention graph branch and the contradictory subgraph, a response difference chain is generated; Extract the continuous portion of the response difference chain that extends from the candidate governance unit to the target improved state or the unfavorable state of the associated river segment, and generate a counterfactual response path.

6. The method for generating river management schemes based on river morphology evolution maps according to claim 5, characterized in that, The combined counterfactual response paths, which have spatial and temporal connections, identify the convergence, cancellation, and amplification of effects that occur after combination, generating a combined response structure including: Arrange the counterfactual response paths in the order in which they arrive at the same river morphology change unit to generate a path temporal relationship; By combining the path temporal relationships to compare the evolution directions reaching the same river morphology change unit, path interaction relationships are generated; Identify path combinations with the same evolution direction and increasing transmission degree, path combinations with opposite evolution direction and decreasing transmission degree, and path combinations that end transmission at the same node in the path interaction relationship, and generate amplification relationship, cancellation relationship and convergence relationship; Connect candidate governance units that have the amplification relationship, the cancellation relationship, or the convergence relationship, and retain the corresponding path timing relationship to generate a combined response structure.

7. The method for generating river management schemes based on river morphology evolution maps according to claim 6, characterized in that, The step of using the engineering constraints to trim the candidate governance units and corresponding connections in the combined response structure that cannot be jointly implemented, and generating an implementation diagram, includes: The engineering constraints are mapped to the candidate governance units and path timing relationships in the combined response structure to generate governance unit constraint relationships; Determine whether the constraint relationship of the governance unit has spatial occupancy conflict, implementation timing conflict, or crosses the graph constraint boundary after the effect is amplified, and generate a conflict governance combination; Remove candidate governance units from the conflict governance combination that cause the unfavorable state of the associated river segment to continue to propagate and cannot be offset by other counterfactual response paths, and generate a constraint satisfaction structure. Connect the candidate governance units in the constraint-satisfying structure that jointly reduce the contradictory effects and maintain spatial connectivity to generate a collaborative governance combination; Retaining the aforementioned collaborative governance combination and its corresponding action transmission sequence and construction dependency sequence, a feasible implementation diagram is generated.

8. The method for generating river management schemes based on river morphology evolution maps according to claim 7, characterized in that, The step of analyzing the action transmission order and construction dependency order of the candidate treatment units in the feasible implementation diagram to generate a river treatment scheme includes: Read the action transmission sequence in the proposed implementation diagram to generate the action sequence of the governance unit; Based on the construction dependency sequence, candidate treatment units that cannot be directly implemented in the order of action of the treatment units are adjusted to generate a constructable treatment sequence. Extract the separation location, change method, and connection object of each candidate treatment unit in the constructable treatment sequence to generate treatment unit implementation information; By combining the implementation information of the governance units with the feasible governance sequence, a river governance plan is generated.

9. The method for generating river management schemes based on river morphology evolution maps according to claim 8, characterized in that, The method further includes: Obtain river monitoring data generated after the implementation of the river management plan, and generate a river state sequence after implementation; Compare the implemented river state sequence with the river state change units corresponding to the counterfactual response path to generate path response differences; Locate the causal propagation edge where the path response difference first appears in the river morphology evolution map, and generate the deviation propagation location; Based on the deviation transmission location, the transmission direction or degree of the corresponding causal transmission edge is corrected to generate a corrected river morphology evolution map. The revised river management plan was revised using the revised river morphology evolution map.

10. A river management scheme generation system based on river morphology evolution maps, characterized in that, The system includes: The state sequence generation module is used to acquire multi-period river observation data of the target river, perform location registration and time sorting on the multi-period river observation data, and generate a river state sequence. The evolution map generation module is used to identify river morphology change units with continuous spatial range and consistent evolution direction in the river morphology state sequence, connect the river morphology change units with sequential response relationship, and generate a river morphology evolution map. The conflict localization module is used to obtain the governance objectives and engineering constraints of the target river channel, map the governance objectives and engineering constraints to the river morphology evolution map to extract conflicts, and generate a conflict sub-map; The counterfactual processing module is used to construct candidate governance units based on the transmission direction or degree of river morphology change units in the contradictory subgraph, embed the candidate governance units one by one into the contradictory subgraph and compare the changes in the graph before and after embedding, and generate counterfactual response paths. The combined response generation module is used to combine the counterfactual response paths that have spatial and temporal connections, identify the convergence, cancellation, and amplification of effects that occur after combination, and generate a combined response structure. The scheme diagram generation module is used to use the engineering constraints to trim the candidate governance units and corresponding connection relationships that cannot be implemented together in the combined response structure, and generate a scheme diagram. The scheme output module is used to parse the action transmission order and construction dependency order of the candidate treatment units in the feasible scheme diagram and generate a river treatment scheme.

Images