A three-dimensional automatic wiring layout method for pre-embedded cooling water pipes of a wind power foundation

Abstract

The present application relates to the technical field of line layout, in particular to a three-dimensional automatic wiring layout method for a wind power foundation pre-buried cooling water pipe. The present application divides the entity space of the wind power foundation into multiple space units; analyzes the heat accumulation, back-lag and temperature control deviation of the space units during the temperature control period, determines the comprehensive time sequence overheat dose of the space units; compares the comprehensive time sequence overheat dose, distance, obstacle component proportion between the space units and the corresponding passable space units, and the channel continuity between the space units, determines the heat disaster corridor value of the passable space units; combines the heat disaster corridor value and the geometric passing consumption to determine the path unit cost; based on the path unit cost, generates the current cooling water pipe and constantly updates the data to realize the update of the cooling water pipe path. The present application greatly reduces the risk of temperature control exceeding and temperature cracks, and improves the foundation durability and the safety of the wind turbine operation.

Description

Technical Field

[0001] This invention relates to the field of line layout technology, specifically to a three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations. Background Technology

[0002] Onshore wind turbine foundations typically employ large-volume concrete caps, spread foundations, or foundations with anchored cages. During concrete pouring, the heat of hydration from cement releases heat, causing the internal temperature of the foundation to rise. If the peak internal temperature, the temperature difference between the inside and outside, or the cooling rate are not properly controlled, thermal stress can easily form and induce cracks, thereby affecting the foundation's durability, fatigue resistance, and the long-term operational safety of the wind turbine.

[0003] Pre-embedded cooling water pipes are a common practice in the temperature control construction of large-volume concrete foundations for onshore wind power. Before concrete pouring, cooling water pipes are laid in key heat-generating areas inside the foundation. Water circulation removes the heat of hydration, keeping the internal temperature field within a controllable range. However, since cooling water pipes cannot be adjusted after concrete pouring, their design quality directly affects subsequent temperature control effectiveness, construction rework risks, and foundation formation. Furthermore, onshore wind power foundations are characterized by large-volume concrete, complex structural shapes, dense internal steel reinforcement, and the inability to alter the pre-embedded pipes. Therefore, the layout of cooling water pipes must simultaneously meet requirements for spatial obstacle avoidance, cooling efficiency, construction feasibility, and structural safety.

[0004] Existing methods for laying out pre-embedded cooling water pipes in wind turbine foundations often rely on static temperature fields at a specific moment after pouring, empirically determined hot zones, or standardized pipe spacing. This simplifies temperature control risks to the identification of instantaneous high-temperature areas, failing to adequately reflect the dynamic evolution of the hydration heat of large-volume concrete during the critical temperature control period, including its continuous release, conduction, and attenuation. In actual engineering projects, some spatial areas may not be instantaneously at their highest temperature, but due to prolonged periods of high temperature, delayed cooling, or repeated approaches to temperature control boundaries, significant heat damage can accumulate. Conversely, some instantaneous high-temperature areas may have shorter durations, and their actual temperature control priority may not be the highest. If cooling water pipe densification or path optimization is still based on temperature distribution at a single moment, it can easily lead to a mismatch between cooling resource allocation and the actual heat damage accumulation areas, resulting in insufficient local cooling, uneven temperature gradient control, redundant pipe network layout, or fluctuations in temperature control effectiveness. This directly impacts the reliability of subsequent full-cycle temperature control. Summary of the Invention

[0005] To address the technical problem that relying solely on temperature distribution at a single moment for cooling water pipe routing can negatively impact the reliability of temperature control throughout its lifecycle, this invention aims to provide a three-dimensional automatic routing method for pre-embedded cooling water pipes in wind turbine foundations. The specific technical solution adopted is as follows: Acquire comprehensive engineering data for pre-embedded cooling water pipes in wind power foundations, including temperature data, temperature control target data, and obstacle data. The physical space of the wind power foundation is divided into multiple spatial units; the heat accumulation, fallback hysteresis and temperature control deviation of the spatial units during the temperature control period are analyzed to determine the comprehensive time-series overheating dose of the spatial units. By comparing the combined temporal overheating dose, distance, proportion of obstruction components, and channel continuity between spatial units and their corresponding accessible spatial units, the responsive heat hazard corridor value of the accessible spatial unit is determined. Based on the responsive heat hazard corridor value and geometric passage consumption, the path unit cost is determined; based on the path unit cost, the current cooling water pipe path is generated; Based on the position of the passable space unit in the current cooling water pipe path, the responsive heat hazard corridor value of the passable space unit is iteratively updated to obtain the updated responsive heat hazard corridor value; based on the responsive heat hazard corridor value of the next space unit, the path unit cost is updated, and the cooling water pipe path is updated.

[0006] The embodiments of the present invention have at least the following beneficial effects: This invention first acquires comprehensive engineering data required for the layout of pre-embedded cooling water pipes in wind turbine foundations, and constructs a unified spatial unit adapted for thermal state representation and pipeline access judgment. By transforming instantaneous temperature deviations into heat hazard accumulation results with duration, residual fall phases, and repeated deviation characteristics, it avoids cooling water pipes responding only to a single momentary high temperature point, providing a quantitative basis for subsequent identification of responsive heat hazard corridors and quantifying the comprehensive temporal overheating dose of each spatial unit throughout the entire temperature control cycle. It identifies traversable heat hazard response corridors that can intercept heat hazard transmission, generates the current cooling water pipe path, and continuously updates the data to update the cooling water pipe path. This invention accurately quantifies latent heat hazards through temporal overheating dose, achieving a high degree of matching between cooling resources and actual heat hazard accumulation areas, significantly reducing the risk of temperature control exceeding limits and temperature cracks, and improving foundation durability and wind turbine operational safety. Attached Figure Description

[0007] Figure 1 This is a flowchart illustrating a three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations, provided as an embodiment of the present invention. Detailed Implementation

[0008] The following describes in detail, with reference to the accompanying drawings, a specific scheme for a three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations provided by the present invention.

[0009] Please see Figure 1The diagram illustrates a three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to an embodiment of the present invention. The method includes the following steps: Step S100: Obtain full-dimensional engineering data of the pre-embedded cooling water pipes of the wind power foundation. The full-dimensional engineering data includes: temperature data, temperature control target data, and obstacle data.

[0010] In onshore wind power foundations with central anchor cages, the central anchorage zone typically features concentrated heat, dense reinforcement, a large space occupied by the anchor cage, and difficulty in directly penetrating pipelines. The appropriate location for cooling water pipes is not necessarily the point of highest overheating, but rather a corridor through which heat damage is transferred from the high-heat zone to the workable space.

[0011] First, obtain full-dimensional engineering data for the pre-embedded cooling water pipes of the wind turbine foundation. The full-dimensional engineering data includes: temperature data, temperature control target data, obstacle data, foundation three-dimensional spatial data, and cooling water pipe process constraint data.

[0012] The temperature data is obtained through hydration heat temperature field or temperature field inversion by thermometry, yielding temperature values ​​at different times at different spatial sampling locations. These temperature values ​​constitute the temperature data. The temperature data can be used to reflect the temperature rise, maintenance, and fall process after pouring.

[0013] The temperature control target data is obtained from the project's temperature control design documents or temperature control special plan. The temperature control target data is the ideal temperature data at different sampling locations in the space.

[0014] The obstacle data includes the smallest rectangular outer envelope of obstacle components such as steel bars, anchor bolt cages, and embedded parts. Here, the smallest rectangular outer envelope of the obstacle data is directly used as the occupied space of the obstacle data.

[0015] Cooling water pipe constraint data includes pipe outer diameter, design bending radius, installation tolerance, protective layer requirements, candidate inlet and candidate outlet, which are used to limit the conditions for pipeline passage envelope and continuous routing.

[0016] Step S200: Divide the physical space of the wind power foundation into multiple spatial units; analyze the heat accumulation, hysteresis, and temperature control deviation of the spatial units during the temperature control period, and determine the comprehensive time-series overheating dose of the spatial units.

[0017] Since temperature data typically comes from thermal analysis networks or temperature inversion points, and cooling water pipes must be routed within the actual space constrained by reinforcing bars, anchor cages, embedded parts, and protective layers, directly using the location corresponding to the temperature data as the pipe routing node may result in correct temperature judgments but pipelines that cannot pass through; if only geometric channels are used as the calculation object, information on the evolution of local thermal hazards will be lost.

[0018] Meanwhile, the release and conduction of hydration heat in large-volume concrete is a three-dimensional unsteady-state heat transfer process in a continuous solid medium; the temperature control function of the cooling water pipe is to form a continuous temperature influence field centered on the pipe through heat exchange between the pipe wall and the surrounding concrete, rather than a point-to-point effect on discrete temperature nodes.

[0019] In onshore wind power foundations with central anchor cages, the core anchoring zone is the area with concentrated hydration heat and the highest priority for temperature control. However, due to the dense anchor cages, steel mesh, and embedded parts, it is fragmented into irregular, only locally continuous passageways for pipelines, resulting in a significant misalignment and disconnect between high-risk heat hazard zones and available pipeline space. Therefore, it is necessary to first construct spatial units that simultaneously accommodate thermal state representation and pipeline passage assessment.

[0020] First, the physical space of the wind turbine foundation is divided into multiple spatial units of the same size. It should be noted that the size of the spatial unit is determined by the implementer based on the original resolution of the temperature field, the outer diameter of the cooling water pipes, the design bending radius, and the local clearance formed by the reinforcing steel and anchor bolt cages.

[0021] For each spatial unit, obtain the center coordinates of the spatial unit. Effective volume Adjacent spatial unit set The minimum distance to the smallest rectangular outer envelope of the nearest obstacle component. And whether the conditions for pipeline passage are met.

[0022] For each spatial unit, temperature data at the sampling locations within the unit are acquired, and the average temperature data within the unit is used as the actual temperature data for that unit; based on the sampling time... Spatial Unit For example, spatial units At sampling time The actual temperature data below is denoted as .

[0023] Similarly, temperature control target data is obtained at sampling locations within the space unit. Specifically, the average temperature data within the space unit is used as the temperature control target data for the space unit; based on the sampling time... Spatial Unit For example, spatial units At sampling time The target temperature control data is denoted as follows: .

[0024] Based on the actual temperature data and temperature control target data of each space unit, the temperature control deviation of the space unit is determined. Specifically, for each space unit, the difference between the actual temperature data and the temperature control target data is calculated as the temperature control deviation of the space unit. If the difference between the actual temperature data and the temperature control target data is negative, the temperature control deviation of the space unit is directly set to 0.

[0025] The difference between the actual temperature data and the temperature control target data characterizes the deviation of the actual thermal state from the engineering temperature control target. Truncation of the difference between the actual temperature data and the temperature control target data to 0 is to exclude normal conditions where no overheating risk has occurred, ensuring that subsequent calculations only accumulate the heat damage portion requiring a cooling water pipe response.

[0026] As a specific example, taking the sampling time Spatial Unit For example, spatial unit At sampling time Temperature control deviation The calculation formula is: Where max is the function for maximizing the value. Temperature control deviation represents the positive deviation of the current temperature data from the target temperature control data.

[0027] Temperature control deviation only indicates whether a space unit deviates from the temperature control target at a certain sampling moment. However, the arrangement of cooling water pipes needs to determine whether the space unit continuously requires a cooling response throughout the entire temperature control period. If only the maximum temperature or the maximum temperature control deviation is read, it will still be difficult to identify areas of continuous overheating, delayed cooling, and repeated deviations.

[0028] Therefore, after analyzing the temperature control deviation and obtaining the temperature control deviation amount, the heat accumulation and fallback hysteresis of the space unit during the temperature control period are further analyzed in combination with the temperature control deviation amount in order to determine the comprehensive time-series overheating dose of the space unit.

[0029] The temperature control risk of large-volume concrete for wind power foundations does not depend solely on the instantaneous temperature peak. It is directly related to the temporal thermal hazards throughout the entire hydration heat cycle, including the duration of heat retention, local temperature lag during the cooling and falling phase, and the cumulative effect of temperatures approaching the temperature control threshold at multiple time periods.

[0030] Especially in the foundation center anchor cage area, the concrete volume is concentrated and the heat dissipation path is long. This area may not be the instantaneous highest temperature point during the temperature control period, but it is very easy for heat retention phenomena to deviate from the temperature control target data for a long time. Conventional pipe layout designs that only target short-term high temperature areas are very likely to miss the core control area.

[0031] In this embodiment of the invention, the spatial distribution and severity of this type of time-series thermal hazard are quantitatively identified by full-cycle temperature data, and the core control area for long-term heat retention is accurately delineated, providing a direct quantitative basis for the targeted optimization of the subsequent cooling water pipe layout.

[0032] By taking into account the average magnitude of temperature control deviations at adjacent sampling times during the combined temperature control period, the basic overheat accumulation of the space unit is determined, specifically: First, calculate the time interval between adjacent sampling times, using the sampling time... For example, sampling time Corresponding time interval for: ; Sampling time The next sampling time. Sampling time. This can be referred to as the k-th sampling time. This can be referred to as the (k+1)th sampling moment.

[0033] Then, for each space cell, the basic superheat accumulation of the space cell is determined.

[0034] As a concrete example, the i-th spatial unit Basic superheat accumulation The calculation formula is: Where m represents the number of sampling moments during the temperature control period; Sampling time Temperature control deviation; Sampling time Temperature control deviation; Sampling time The corresponding time interval.

[0035] In the formula for calculating the basic superheat accumulation, This represents the average deviation over the time interval, multiplied by This represents the contribution of thermal hazards caused by sustained deviations within a time interval. This calculation incorporates both the magnitude and duration of the deviation into the pipe layout criteria. The basic superheat accumulation can characterize the spatial unit. The area of ​​deviation time during the temperature control period.

[0036] Furthermore, based on the temperature change trend of the spatial unit at different sampling times, the starting time of temperature drop is determined, specifically: By combining the actual temperature data of all space units, the overall temperature change is determined. Based on the trend of the overall temperature change at different times, the reverse fluctuation amount at different times is obtained; the time corresponding to the minimum value of the reverse fluctuation amount is taken as the starting time of the temperature drop. More specifically: First, calculate the actual temperature data of all spatial cells at each sampling time. The median value is denoted as the overall temperature. Where k = 1, 2, ..., m, and m represents the total number of recording times within the temperature control period. Further calculation is performed on the overall temperature change between adjacent sampling times, which is then used as the overall temperature change. Specifically: calculate the difference between the overall temperature at the next sampling time and the previous sampling time, and use it as the overall temperature change.

[0037] Among them, when the overall temperature change This indicates an overall increase in median temperature, when the overall temperature change is... The time indicates a decrease in the overall median temperature, representing the overall temperature change. The time indicates that the overall median temperature remains constant.

[0038] Then, based on the trend of the overall temperature change at different times, the reverse fluctuation at different times is obtained. Specifically: at each sampling time... As candidate split points, the reverse fluctuation corresponding to each candidate split point The calculation formula is: Among them, the first item Indicates candidate split points The previously observed decrease was a reverse fluctuation during the warming phase; the second item Indicates candidate split points The subsequent increase in volume is a reverse fluctuation during the pullback phase.

[0039] Reverse volatility The smaller the value, the lower the value at the sampling time. When used as a dividing point, the period before the dividing point is more consistent with the overall warming process, while the period after the dividing point is more consistent with the overall cooling process.

[0040] The moment corresponding to the minimum value of the reverse fluctuation is taken as the starting moment of the temperature drop. If there are multiple moments corresponding to the same minimum value, the overall temperature among them is selected. The moment of maximum temperature is taken as the starting point of temperature decline.

[0041] After obtaining the start time of temperature drop, the drop hysteresis of the spatial unit is obtained by combining the average amplitude of temperature control deviation within adjacent sampling times.

[0042] As a specific example of the present invention, spatial unit retracement The calculation formula is: Where s represents the starting time of the overall temperature drop. The sequence number in the temperature time series.

[0043] pullback lag Characterizes spatial units The larger the accumulated overheating after the overall decline phase, the greater the stagnation of heat release at that location.

[0044] Then, based on the temperature control deviation, the deviation period in which the deviation occurred is selected; and by combining the duration and location of the deviation period with the temperature control deviation, the deviation contribution is determined.

[0045] spatial units Continuously meet temperature control deviation Time periods with values ​​greater than zero are identified as deviation periods; for example, the q-th deviation period is denoted as... The duration of the deviation period is recorded as . The average temperature control deviation during the deviation period was The total number of deviation periods is In this deviation period, the temperature control deviation is greater than zero at all times. That is, the deviation period is composed of consecutive times where the temperature control deviation is greater than zero. Therefore, the duration of the deviation period is also the number of sampling times within the deviation period.

[0046] It should be noted that when no deviation from the time period is detected, the spatial unit is directly set. Deviation contribution It is zero.

[0047] When a deviation period is detected, the contribution of the deviation is determined by combining the duration and location of the deviation period with the amount of temperature control deviation.

[0048] As a specific example of the present invention, spatial unit Deviation contribution The calculation formula is: ;in, This indicates the relative sequential position of the deviation time period within the deviation sequence of the same spatial unit.

[0049] Finally, by combining the basic superheat accumulation, the fallback hysteresis, and the deviation contribution, the comprehensive time-series superheat dose of the space unit is determined. Specifically, based on the interquartile range of the basic superheat accumulation, the interquartile range of the fallback hysteresis, and the interquartile range of the deviation contribution of all space units, the contribution ratios of the basic superheat accumulation, the fallback hysteresis, and the deviation contribution are determined respectively. The comprehensive time-series superheat dose of each space unit is obtained by weighted summing of the basic superheat accumulation, the fallback hysteresis, and the deviation contribution using the contribution ratios.

[0050] Among them, the integrated temporal overheating dose of the space unit The calculation formula is: ; ; ; ; Basic superheat accumulation Interquartile difference; For the lag in the decline Interquartile difference; For deviation contribution The interquartile difference is used to determine the contribution ratios of basic overheat accumulation, fallback hysteresis, and deviation contribution based on spatial discrimination capability. Basic superheat accumulation The proportion of contribution; The lag in the decline The proportion of contribution; For deviation contribution The proportion of contribution.

[0051] Integrated time-series overheating dose Characterizes spatial units During the temperature control period, the comprehensive heat damage accumulation intensity required for the cooling water pipes to respond changes the pipe layout basis from static high temperature points to full-time heat damage accumulation areas.

[0052] Step S300: Compare the integrated time-series overheating dose, distance, proportion of obstacle components, and channel continuity between spatial units and corresponding accessible spatial units to determine the responsive heat hazard corridor value of the accessible spatial unit.

[0053] The composite time-series overheating dose can indicate which spatial units have a high accumulation of heat hazards, but spatial units with a high composite time-series overheating dose are not necessarily directly traversable by cooling water pipes. In wind turbine foundations with central anchor cages, high heat hazard areas often overlap with the anchor cages, dense reinforcement, or embedded parts. Directly connecting to areas with high composite time-series overheating doses can easily result in impassable paths; while simply implementing ordinary obstacle avoidance may keep pipelines away from the heat hazard source. Therefore, it is necessary to convert high heat hazard areas into corridors within the constructable space that can intercept the transfer of heat hazards.

[0054] The cooling water pipes are not required to pass through the center of the heat source, but rather to be placed along the path of heat transfer from the high-heat-hazard area outwards, so that the heat is carried away by the circulating water during conduction. For the central anchor cage area, a reasonable path is usually located outside the anchor cage or in the continuous area of ​​the rebar clearance, rather than penetrating the anchor core.

[0055] First, for spatial units Select the corresponding passable spatial units from the other spatial units, specifically: Spatial Unit The passable space unit must meet two conditions: first, the passable space unit is within the neighborhood of a space unit; second, the sum of the total time-series superheat dose of the passable space unit minus the sum of the total time-series superheat dose of the space unit is a positive value. More specifically, condition one requires that the distance between the center coordinates of the space unit and the passable space unit be less than a preset spatial distance. This preset spatial distance is 3 meters, but in other embodiments, the implementer can adjust this preset spatial distance according to the actual situation.

[0056] By comparing the combined temporal overheating dose, distance, and proportion of obstruction components between a space unit and its corresponding accessible space unit, the thermal hazard response value of the accessible space unit is determined; specifically: By comparing the combined temporal superheat dose of a space cell with that of its corresponding traversable space cell, the superheat dose difference between the space cells is obtained. Specifically: for spatial units and the corresponding passable space unit Calculate the superheated dose difference If the overheating dose difference A value greater than zero indicates that from a spatial unit To accessible space unit There are candidate directions for release from high heat damage accumulation areas to low heat damage accumulation areas.

[0057] According to spatial units and the corresponding passable space unit Distance between and the proportion of obstacle components in the spatial unit Determine the equivalent transmission distance The corresponding calculation formula is: ;in, The size and distance of the current spatial unit The ratio is used to avoid situations where the denominator is zero.

[0058] Among them, the proportion of obstacle components is spatial unit. and the corresponding passable space unit The proportion of the obstacle's width to the total length of the line connecting the center points. Obstacle component proportions. The higher the value, the greater the equivalent transfer distance, and the less effective the characterization of heat damage response by the cooling water pipe in that direction.

[0059] Based on the trend of combined temporal superheat dose variation between spatial units and their corresponding traversable spatial units, the dose reduction continuity of spatial units is determined. , specifically: For spatial units To accessible space unit The sequence of spatial units traversed sequentially along the line connecting their center coordinates is read, and the combined temporal superheat dose of each spatial unit in the sequence is recorded; let the total number of adjacent spatial unit pairs in the spatial unit sequence be... The statistical analysis identifies adjacent spatial unit pairs that satisfy a decreasing continuous relationship in terms of the overall time-series superheat dose. Specifically, if the overall time-series superheat dose of the preceding spatial unit is not lower than that of the following spatial unit, then the number of all adjacent spatial unit pairs satisfying this condition is: .

[0060] To mitigate the impact of local fluctuations in the discrete temperature field, the interquartile range of the integrated time-series superheat dose along this path is used. As a natural fluctuation tolerance; when the combined time-series superheat dose of the subsequent space unit is higher than that of the previous space unit, but does not exceed the combined time-series superheat dose value of the previous space unit and the natural fluctuation tolerance. When the sum is equal to the sum, the decreasing continuous relationship is still included; when it exceeds this range, it is considered that there is a significant increase in reverse heat damage in that direction, and the decreasing continuous relationship is not included.

[0061] As a specific example of the present invention, spatial unit To accessible space unit Dose reduction continuity The calculation formula is: .

[0062] Dosage Decrease Continuity Characterized from to Continuity of dose reduction in direction. The larger the value, the more the overheating dose in that direction conforms to the characteristic of decreasing heat transfer from the high heat hazard area outwards. The more suitable the location is as a passageway for intercepting heat transfer in cooling water pipes.

[0063] By combining the equivalent transmission distance, the dose reduction continuity, and the comprehensive time-series overheating dose, the thermal hazard response value of the passable space unit is determined.

[0064] Using any accessible space unit as the target accessible space unit; obtaining the comprehensive time-series overheat dose of the accessible space unit corresponding to the target accessible space unit; obtaining the equivalent transfer distance and dose reduction continuity of the accessible space unit corresponding to the target accessible space unit; determining the thermal hazard response value of the target accessible space unit based on the equivalent transfer distance, dose reduction continuity, and comprehensive time-series overheat dose of the target accessible space unit and its corresponding accessible space unit; wherein, the comprehensive time-series overheat dose and dose reduction continuity are positively correlated with the thermal hazard response value, and the equivalent transfer distance is negatively correlated with the thermal hazard response value.

[0065] As a specific example of the present invention, a passable space unit Taking the target accessible space unit as an example, calculate the accessible space unit. Thermal damage response value : ;in, Indicates a passable space unit The corresponding number of passable units; As a passable space unit The integrated temporal overheating dose of the corresponding m-th traversable spatial unit; As a passable space unit The dose reduction continuity up to the corresponding m-th traversable spatial unit; As a passable space unit The equivalent transmission distance to its corresponding m-th passable spatial unit.

[0066] Accessible space unit The thermal damage response value no longer represents ordinary spatial proximity, but is used to quantify traversable spatial units. The degree of interception and adaptation for the transfer of high-time-series overheating doses to the workable space. Thermal damage response value. The larger the value, the more it indicates At the same time, it is closer to the source of high-sequence overheating, has a stronger continuity of overheating dose reduction in its direction, and has a smaller equivalent transmission distance after being blocked by obstacles. Therefore, it has a higher thermal hazard response value as a candidate passage location for cooling water pipes.

[0067] The continuity of passageways between the corresponding accessible spatial units is then calculated. For passable space units Read its adjacent spatial units and filter out those that can be accessed by passing through. Adjacent passable space units are formed for pipeline splicing. The pipeline splicing must simultaneously meet the following requirements: the pipelines between the two lines do not spatially conflict with the boundaries of reinforcing bars, anchor cages, embedded parts, and protective layers; the minimum clearance in the splicing direction can accommodate the outer diameter of the cooling water pipe, installation tolerance, and protective layer requirements; and the design bending radius requirement of the cooling water pipe can be met when switching from the previous splicing direction to this splicing direction.

[0068] As an example of the present invention, a passable space unit Channel continuity The calculation formula is: ;in, Indicates the ability to interact with accessible space units The number of adjacent spatial units that form an effective pipeline connection; express Total number of adjacent accessible space units.

[0069] Indicates a passable space unit The degree of fit as a component of a continuous cabling channel; The larger, the more it means The more readily available the surrounding area is for continuous pipeline extension, the more suitable it is as a location for cooling water pipe paths.

[0070] Finally, by combining the thermal hazard response value of the accessible space unit with the corresponding passage continuity of the accessible space unit, the responsive thermal hazard corridor value is determined. (The thermal hazard response value is missing from the original text.) and channel continuity Based on this, a responsive heat hazard corridor value is formed. : Heat hazard corridor value Characterized accessible spatial units The overall fit strength as a candidate path location for cooling water pipes. If Higher but A lower value indicates that while the location has thermal hazard response significance, it is difficult to form a suitable location for pipeline construction; if Higher but A lower value indicates that while the location is easy to route, its contribution to preventing heat transfer is limited. Only when both are high... This will improve the priority of pipeline deployment, thereby obtaining a higher priority in subsequent path cost calculations.

[0071] Step S400: Combine the responsive heat hazard corridor value and geometric passage consumption to determine the path unit cost; based on the path unit cost, generate the current cooling water pipe path.

[0072] Once cooling water pipes are placed near a heat hazard corridor, they will cool the surrounding concrete, reducing the necessity of laying new pipelines in that area. Conversely, areas farther from existing pipelines that still have high heat hazard response value should receive higher priority in subsequent routing. Without updating remaining heat hazard corridors, the algorithm will continue to be attracted by initial high-value areas, leading to localized pipeline density and uneven overall cooling coverage. Therefore, if all cooling water pipe paths are generated at once based solely on the initial heat hazard corridor values, the paths are prone to repeatedly clustering near the same high-heat-hazard corridor, while other persistent heat hazard areas not yet covered by cooling water pipes may still lack response.

[0073] This invention employs a progressive processing method of path generation, reduction of heat hazard corridors, and regeneration of remaining corridors, so that the pipeline layout gradually converges as the heat hazard response requirements change.

[0074] First, the initial heat hazard corridor value is used. As the initial residual heat hazard corridor value ; For passable space units The path unit cost is determined by combining the responsive heat hazard corridor value and geometric passage consumption.

[0075] As an example, a passable space unit Path unit cost The calculation formula is: ;in, This indicates that the cooling water pipe passes through a passable space unit. The overall cost; Indicates a passable space unit The geometrical travel cost is determined at least by the spatial clearance, turning feasibility, and path detour length; This represents the initial residual heat hazard corridor value; Indicates the spatial unit in the current state The remaining heat hazard corridor value that has not yet been responded to by cooling water pipes is used to guide paths to prioritize corridor areas that still have high heat hazard response requirements. Geometric passage consumption. This is used to prevent the path from entering spaces that are difficult to construct.

[0076] As a specific example of the present invention, the method for obtaining geometrical passage consumption is as follows: The geometrical passage consumption of the corresponding currently passable spatial unit is calculated by weighted summing the normalized value of the reciprocal of the obstacle clearance distance of the passable spatial unit, the normalized value of the length of the line connecting the center of the passable spatial unit to the previous corresponding spatial unit, and the normalized value of the turning angle in the current expansion direction. Each weight is pre-set according to engineering objectives such as construction obstacle avoidance priority, single-loop length control requirements, and cooling water pipe bending radius control requirements. As an example, the reciprocal of the obstacle clearance distance, the length of the line connecting the center, and the turning angle in the current expansion direction can also be assigned the same weight, for example, 1 / 3, and the sum of the weights of these three features should be 1. For spatial units or expansion directions that do not meet the requirements for cooling water pipe outer diameter, protective layer thickness, or minimum bending radius, their corresponding geometrical passage consumption is assigned an infinite cost. The normalization method used in calculating geometric passage consumption is the maximum-minimum normalization method. The maximum value is the maximum value of the corresponding data in the historical database, and the minimum value is the minimum value of the corresponding data in the historical database. For example, when performing maximum-minimum normalization on the inverse of the obstacle net distance of the passable space unit, the maximum value is the maximum value of the inverse of the obstacle net distance of the passable space unit in the historical database, and the minimum value is the minimum value of the inverse of the obstacle net distance of the passable space unit in the historical database.

[0077] More specifically: the reciprocal of the obstacle clearance is the reciprocal of the minimum distance from the current passable space cell to the nearest obstacle surface, and the turning angle in the current expansion direction is the deflection angle of the line connecting the center of the current passable space cell and the previous space cell relative to the horizontal direction.

[0078] Multiple paths are obtained by connecting accessible spatial units. More specifically: Accessible spatial units are connected to their corresponding accessible spatial units, and then, based on the already connected accessible spatial units, their corresponding accessible spatial units are connected again to obtain multiple paths. The path with the lowest cumulative cost is selected as the current cooling water pipe path, and the inlet candidate region and outlet candidate region are used as the path's start and end constraints. Here, the inlet candidate region serves as the initial spatial unit, and subsequent accessible spatial units are developed based on this spatial unit. It should be noted that a single generated cooling water pipe path is strictly restricted from repeatedly passing through the same accessible spatial unit.

[0079] Step S500: Based on the position of the passable space unit in the current cooling water pipe path, iteratively update the responsive heat hazard corridor value of the passable space unit to obtain the updated responsive heat hazard corridor value; based on the responsive heat hazard corridor value of the next space unit, update the path unit cost and update the cooling water pipe path.

[0080] After generating the current cooling water pipe path, the thermal response reduction ratio is first determined based on the position of the passable space unit in the current cooling water pipe path.

[0081] Calculate accessible space units The shortest distance to the center point of the current cooling water pipe path Based on this, the current cooling water pipe path is determined for accessible space units. Thermal response reduction ratio The thermal response reduction ratio serves only to characterize the degree to which existing pipelines reduce the demand for surrounding heat hazard corridors. As a preferred embodiment of the invention, the thermal response reduction ratio... The calculation formula is: ,in, The effective cooling influence radius of the cooling water pipes is pre-set by the implementer through analysis and observation of historical data.

[0082] The respondable heat hazard corridor value is iteratively updated by using the heat response reduction ratio to obtain the updated respondable heat hazard corridor value.

[0083] Updated responsive heat hazard corridor value The calculation formula is: ; in, This indicates the passable spatial unit after the current cooling water pipe path is generated. The residual heat hazard corridor value; This represents the remaining heat hazard corridor value before the update, which is also the responsive heat hazard corridor value before the update. This indicates the current cooling water pipe path to this accessible space unit. The reduction ratio of heat hazard corridor demand. The closer a spatial unit is to an existing pipeline, the more its heat hazard corridor demand is reduced, and the weaker its attraction in subsequent path searches; areas that are not yet effectively covered have a higher responsive heat hazard corridor value, thus continuing to guide subsequent paths to extend into that area.

[0084] After updating the remaining heat hazard corridor values, the path unit cost is recalculated, and subsequent pipeline paths are generated. This continues until the main remaining heat hazard corridor areas are covered, or the effect of subsequent newly added paths on reducing remaining heat hazard corridors is no longer significant. For the generated discrete paths, their closed-loop hydraulic verification, collision verification, minimum bending radius verification, continuous curve generation, and construction drawing refinement can be completed using existing pipeline integrated design and construction refinement technologies.

[0085] The responsive heat hazard corridor values ​​for each accessible space unit are obtained. Updated responsive heat hazard corridor values Path unit cost Based on the three-dimensional cooling water pipe routing results generated by updating the residual heat hazard corridor, a three-dimensional automatic routing layout scheme for pre-embedded cooling water pipes in wind power foundations is formed.

[0086] The three-dimensional automatic routing layout scheme for pre-embedded cooling water pipes in wind turbine foundations includes at least: the three-dimensional centerline of the cooling water pipes, the connection relationship between the inlet and outlet, the spatial unit sequence corresponding to each pipe path, the responsive heat hazard corridor area covered by the path, the remaining heat hazard corridor value after path generation, and pipe positioning data for construction detailing. The three-dimensional centerline is determined by the path unit cost. The generated area has a relatively high residual heat hazard corridor value. Therefore, it can prioritize responding to persistent heat-damaged areas that have not yet been adequately cooled and covered, rather than repeatedly focusing on a single instantaneous high-temperature location.

Claims

1. A three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind turbine foundations, characterized in that, The method includes the following steps: Acquire comprehensive engineering data for pre-embedded cooling water pipes in wind power foundations, including temperature data, temperature control target data, and obstacle data. The physical space of the wind power foundation is divided into multiple spatial units; the heat accumulation, fallback hysteresis and temperature control deviation of the spatial units during the temperature control period are analyzed to determine the comprehensive time-series overheating dose of the spatial units. By comparing the combined temporal overheating dose, distance, proportion of obstruction components, and channel continuity between spatial units and their corresponding accessible spatial units, the responsive heat hazard corridor value of the accessible spatial unit is determined. Based on the responsive heat hazard corridor value and geometric passage consumption, the path unit cost is determined; based on the path unit cost, the current cooling water pipe path is generated; Based on the position of the passable space unit in the current cooling water pipe path, the responsive heat hazard corridor value of the passable space unit is iteratively updated to obtain the updated responsive heat hazard corridor value; based on the responsive heat hazard corridor value of the next space unit, the path unit cost is updated, and the cooling water pipe path is updated.

2. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 1, characterized in that, The analysis examines the heat accumulation, hysteresis, and temperature control deviation of the space unit during the temperature control period to determine the comprehensive time-series superheat dose of the space unit, including: Based on the actual temperature data and temperature control target data of each space unit, determine the temperature control deviation of the space unit; The basic overheat accumulation of a space unit is determined by the average magnitude of the temperature control deviation at adjacent sampling times during the combined temperature control period. The starting time of temperature drop was determined based on the temperature change trend of the spatial unit at different sampling times. Starting from the initial moment of temperature drop, the average magnitude of temperature control deviation within adjacent sampling moments is used to obtain the drop hysteresis of the spatial unit; Based on the temperature control deviation, the deviation periods in which the deviation occurred are selected; and the deviation contribution is determined by combining the duration and location of the deviation periods with the temperature control deviation. The overall time-series overheating dose of the space unit is determined by combining the basic overheating accumulation, the fallback hysteresis, and the deviation contribution.

3. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 2, characterized in that, The determination of the temperature drop start time based on the temperature change trend of the spatial unit at different sampling times includes: By combining the actual temperature data of all space units, the overall temperature change is determined; based on the trend of the overall temperature change at different times, the reverse fluctuation at different times is obtained; the time corresponding to the minimum value of the reverse fluctuation is taken as the starting time of temperature drop.

4. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 2, characterized in that, The determination of the comprehensive temporal superheat dose of the space unit by combining the basic superheat accumulation, the fallback hysteresis, and the deviation contribution includes: Based on the interquartile range of the basic superheat accumulation, the interquartile range of the fallback hysteresis, and the interquartile range of the deviation contribution of all space units, the contribution ratios of the basic superheat accumulation, the fallback hysteresis, and the deviation contribution are determined respectively. The comprehensive time-series overheating dose for each space unit is obtained by weighted summing of the basic overheating accumulation, the fallback hysteresis, and the deviation contribution for each space unit using the contribution ratio.

5. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 1, characterized in that, The comparison of the spatial unit and the corresponding accessible spatial unit, considering the comprehensive temporal overheating dose, distance, proportion of obstruction components, and channel continuity between spatial units, determines the responsive heat hazard corridor value of the accessible spatial unit, including: By comparing the combined time-series overheating dose, distance, and proportion of obstruction components between the space unit and the corresponding accessible space unit, the thermal hazard response value of the accessible space unit is determined. Determine the continuity of passageways between accessible spatial units corresponding to spatial units. ; By combining the thermal hazard response value of the accessible space unit with the corresponding passage continuity of the accessible space unit, the responsive thermal hazard corridor value is determined.

6. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 5, characterized in that, The comparison of the combined time-series overheating dose, distance, and proportion of obstruction components between the comparative spatial unit and the corresponding accessible spatial unit determines the thermal hazard response value of the accessible spatial unit, including: The equivalent transmission distance is determined based on the distance between the spatial unit and the corresponding passable spatial unit, as well as the proportion of obstacle components in the spatial unit. Based on the trend of the combined time-series overheat dose between the spatial unit and the corresponding traversable spatial unit, the dose reduction continuity between the spatial unit and the corresponding traversable spatial unit is determined. By combining the equivalent transmission distance, the dose reduction continuity, and the comprehensive time-series overheating dose, the thermal hazard response value of the passable space unit is determined.

7. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 6, characterized in that, The determination of dose reduction continuity between a spatial unit and its corresponding accessible spatial unit based on the trend of changes in the comprehensive temporal superheat dose between the spatial unit and the corresponding accessible spatial unit includes: For spatial units to passable spatial units, read the sequence of spatial units traversed sequentially along the line connecting their center coordinates; The proportion of adjacent spatial unit pairs in the spatial unit sequence whose integrated temporal overheating dose satisfies a decreasing continuity relationship is calculated as the dose decreasing continuity between spatial units and traversable spatial units.

8. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 6, characterized in that, The determination of the thermal hazard response value of a passable space unit by combining the equivalent transfer distance, the dose reduction continuity, and the comprehensive time-series overheat dose includes: Use any passable spatial unit as the target passable spatial unit; Obtain the comprehensive temporal overheating dose of the accessible space unit corresponding to the target accessible space unit; Obtain the equivalent delivery distance and dose reduction continuity of the accessible space unit corresponding to the target accessible space unit; Based on the equivalent transfer distance, dose reduction continuity, and comprehensive temporal overheat dose of the target accessible space unit and its corresponding accessible space unit, the thermal hazard response value of the target accessible space unit is determined. Among them, the comprehensive temporal overheat dose and dose reduction continuity are positively correlated with the thermal hazard response value, while the equivalent transfer distance is negatively correlated with the thermal hazard response value.

9. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 1, characterized in that, The step of generating the current cooling water pipe path based on the path unit cost includes: Multiple paths are obtained by connecting accessible space units; the path with the lowest cumulative cost of the path units is selected as the current cooling water pipe path, and the inlet candidate area and outlet candidate area are used as the path start and end constraints.

10. The three-dimensional automatic routing layout method for pre-embedded cooling water pipes in wind power foundations according to claim 1, characterized in that, The step of iteratively updating the responsive heat hazard corridor value of the accessible space unit based on its position in the current cooling water pipe path to obtain the updated responsive heat hazard corridor value includes: The thermal response reduction ratio is determined based on the location of the accessible space unit in the current cooling water pipe path; The respondable heat hazard corridor value is iteratively updated by using the heat response reduction ratio to obtain the updated respondable heat hazard corridor value.

Images