Digital twinborn body construction method for urban planning

Abstract

The invention relates to the technical field of data processing, in particular to a digital twinborn body construction method for urban planning, which comprises the following steps: acquiring the daily yield of each industrial element; obtaining a time delay consistency coefficient according to the yield change consistency condition of the industrial element and other industrial elements and the time delay difference under the yield consistency condition; obtaining all upstream industrial elements of each downstream industrial element based on the time delay consistency coefficient; obtaining the yield change consistency according to the covariant relationship of the yields under the condition that the upstream industrial elements are used as the same link raw materials; obtaining a plurality of industrial system elements according to the yield change consistency coefficient and the time delay consistency coefficient; and constructing an urban industrial system network according to all industrial system elements, and obtaining a plurality of digital twin bodies for urban planning. The invention aims to solve the problem that multiple links and products in the urban industrial system are mutually associated, so that the divided urban industrial systems do not interfere with each other.

Description

Technical Field

[0001] The present invention relates to the field of data processing technology, and in particular to a method for constructing a digital twin for urban planning. Background Art

[0002] Digital twin cities refer to the construction of city information models by combining the physical and digital spaces of cities through information technology and intelligent equipment. Through digital twin technology, urban planners can more accurately simulate the future development of cities, predict possible problems, and develop more scientific planning schemes. Digital twin technology is mainly used in urban planning to simulate the production capacity and consumption of cities, and then predict the future direction of cities. Among them, the economic lifeline of cities is usually affected by multiple industrial chains. By analyzing and predicting the production capacity of each industrial chain in the future, the future direction of cities can be obtained, and the planning of future cities can be adjusted.

[0003] The industrial system of a city has multiple complete industrial chains, from raw materials to semi-finished products, from semi-finished products to finished products, from semi-finished products to new semi-finished products, and from raw materials to finished products, which are all related in sequence. Therefore, when analyzing the future direction of the city, due to the existence of local industries, outbound industries, processing industries, etc. in the city, it is necessary to separate different industrial chains and then construct a digital twin model. However, in the process of separation, it is affected by the time mismatch of different industrial chains and the use of one product in multiple industrial system links, resulting in the upstream and downstream relationships of some industrial system links being hidden, which in turn reduces the accuracy of the constructed digital twin model in predicting the city's future. Summary of the invention

[0004] The present invention provides a digital twin construction method for urban planning to solve the existing problems.

[0005] A digital twin construction method for urban planning of the present invention adopts the following technical solution: One embodiment of the present invention provides a method for constructing a digital twin for urban planning, the method comprising the following steps: Get the daily output of each industrial element in the preset date sequence; According to the consistency of the output change of each industrial element with other industrial elements and the difference of time lag in the output under the consistency condition, the time lag consistency coefficient of each industrial element with other industrial elements is obtained; Based on the time-lag consistency coefficient, the downstream industrial elements of each industrial element are screened to obtain all the upstream industrial elements of each downstream industrial element; According to the covariation relationship of the output between all the upstream industrial elements of each downstream industrial element as the raw materials of the same link, the consistency of the output change of each upstream industrial element under each downstream industrial element is obtained; According to the consistency of the output change of each upstream industrial element under the downstream industrial element, combined with the time-lag consistency coefficient of the corresponding industrial element, several industrial system elements are obtained; An urban industrial system network is constructed based on all industrial system elements. After dividing the urban industrial system network into communities, several industrial communities are obtained. Based on each industrial community, several digital twins for urban planning are obtained.

[0006] Preferably, the specific steps of obtaining the time-lag consistency coefficient include: For the i-th industrial element, obtain all production date segments of the i-th industrial element and all matching production date segments in the m-th suspected downstream industrial element of the i-th industrial element; The time-lag consistency coefficient between the i-th industrial element and the m-th suspected downstream industrial element The calculation method is: in, represents the variance function; is an exponential function with a natural constant as base; is the median date of the k-th production date segment of the i-th industrial element, It is the median date of the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

[0007] Preferably, the specific steps of obtaining the production date segment include: Use the dynamic time warping algorithm to dynamically match the daily output of the i-th industrial element with the j-th industrial element except the i-th industrial element, and obtain the daily output matching day of the i-th industrial element on the j-th industrial element; If the date sequence value of the second day's output matching day of the i-th industrial element is greater than or equal to 2, the j-th industrial element is recorded as the suspected downstream industrial element of the i-th industrial element; Obtain all the minimum values ​​of the sequence composed of the daily output of the ith industrial element, and record the date sequence composed of all dates between two adjacent minimum values ​​and including the first minimum value of the two minimum values ​​as a production date segment of the ith industrial element.

[0008] Preferably, the specific steps of obtaining the matching production date segment include: The first date in the k-th production date segment of the i-th industrial element, the production matching date in the date sequence of the m-th suspected downstream industrial element of the i-th industrial element, is recorded as the first day of the k-th matching production date segment in the m-th suspected downstream industrial element of the i-th industrial element; The last date in the k-th production date segment, the production matching date in the date sequence of the m-th suspected downstream industrial element of the i-th industrial element, is recorded as the last day of the k-th matching production date segment in the m-th suspected downstream industrial element of the i-th industrial element; The sequence consisting of all dates between the first day and the last day of the kth matching production date segment in the mth suspected downstream industrial element is recorded as the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

[0009] Preferably, the specific steps of obtaining the upstream industrial elements include: A time lag consistency threshold is preset, and when the time lag consistency coefficient between the ith industrial element and the mth suspected downstream industrial element is greater than or equal to the time lag consistency threshold, the mth suspected downstream industrial element of the ith industrial element is recorded as the downstream industrial element of the ith industrial element; For the nth downstream industrial element of the ith industrial element, the ith industrial element is recorded as the nth downstream industrial element. An upstream industrial element.

[0010] Preferably, the specific steps of obtaining the consistency of the output change include: The time-lag consistency coefficient between the ith industrial element and the nth downstream industrial element of the ith industrial element is recorded as the time-lag consistency coefficient between the nth downstream industrial element and the nth downstream industrial element. The time-lag consistency coefficient of the upstream industrial elements; In the nth downstream industrial element In the date sequence of each covariant contributing industrial element of the upstream industrial element, obtain the date sequence of the first The sequence segment with the same production date segment date as the kth production date segment of the nth downstream industrial element is recorded as the The k-th production date segment under each covariate contribution industrial element of the upstream industrial elements; The nth downstream industrial element Consistency of production changes of upstream industrial elements The calculation method is: in, To find the function of Pearson correlation coefficient; The nth downstream industrial element The number of covariant contributing industrial elements of the upstream industrial elements, The nth downstream industrial element The number of production date segments of each upstream industrial element; The nth downstream industrial element The sequence of the production of the kth production date segment of the upstream industrial elements, The nth downstream industrial element The sequence of the covariant contributions of the upstream industrial elements and the production of the industrial elements in the k-th production date segment.

[0011] Preferably, the specific steps of obtaining the covariant contributing industrial elements include: Among all the upstream industrial elements of the nth downstream industrial element, except the The other upstream industrial elements except the first upstream industrial element are recorded as the nth downstream industrial element. The covariant contribution of upstream industrial elements is the industrial element.

[0012] Preferably, the specific steps of obtaining the industrial system elements include: Obtain the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element, perform linear normalization on the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element to obtain the normalized value of the industrial system correlation coefficient of each upstream industrial element under the nth downstream industrial element, use the Otsu threshold algorithm to perform maximum difference segmentation on the normalized value of the industrial system correlation coefficient, and use all upstream industrial elements in the largest part of the two after the segmentation as several upwardly associated industrial elements of the nth downstream industrial element; The industrial system correlation coefficient of the upstream industrial element corresponding to each upwardly associated industrial element of the nth downstream industrial element is used as the industrial system correlation coefficient of each upwardly associated industrial element of the nth downstream industrial element; each downstream industrial element and each of its upwardly associated industrial elements constitute an industrial system element.

[0013] Preferably, the specific steps of obtaining the industrial system correlation coefficient include: The nth downstream industrial element Industrial system correlation coefficient of upstream industrial elements The calculation method is: in, The relationship between the nth downstream industrial element and the nth downstream industrial element The time-lag consistency coefficient of the upstream industrial elements, For the nth downstream industrial element The consistency of output changes of upstream industrial elements.

[0014] Preferably, the step of constructing an urban industrial system network according to all industrial system elements, dividing the urban industrial system network into communities to obtain a number of industrial communities, and obtaining a number of digital twins for urban planning based on each industrial community includes: Constructing an urban industrial system network based on a relational network structure, wherein the urban industrial system network is a topological structure, wherein each node in the urban industrial system network represents the first industrial element, and if two nodes can constitute an industrial system element, the industrial system correlation coefficient between the downstream industrial element and the upwardly associated industrial element in the industrial system element is used as the degree of the urban industrial system network; The urban industrial system network is divided into communities through the Leuven algorithm to obtain several industrial communities. The output of industrial elements in each community on all dates is used as the data basis for generating a digital twin, thus obtaining several digital twins for urban planning.

[0015] The beneficial effects of the technical solution of the present invention are as follows: the present application obtains the daily output of each industrial element on a preset date sequence; obtains the time lag consistency coefficient between each industrial element and other industrial elements according to the consistency of the output changes between each industrial element and other industrial elements and the difference in time lag when the output is consistent; highlights the consistency of output changes of different industrial elements and the lag effect of output changes in different batches of production between upstream and downstream links in the industrial system through the time lag consistency coefficient; screens the downstream industrial elements of each industrial element based on the time lag consistency coefficient to obtain all upstream industrial elements of each downstream industrial element; obtains the consistency of output changes of each upstream industrial element under each downstream industrial element according to the covariation relationship of the output between all upstream industrial elements of each downstream industrial element when they are raw materials of the same link; and analyzes all upstream industrial elements of the same downstream industrial element. The simultaneous changes between industrial elements are used to analyze the consistency of output changes between different upstream industrial elements in the same upstream link, thereby avoiding the different degrees of influence of different proportions of different industrial elements in the same upstream link on the downstream industrial link, resulting in a decrease in the time-lag consistency coefficient; according to the consistency of output changes of each upstream industrial element under the downstream industrial element, combined with the time-lag consistency coefficient of the corresponding industrial element, a number of industrial system elements are obtained; an urban industrial system network is constructed based on all industrial system elements, and a number of industrial communities are obtained after the urban industrial system network is divided into communities, and a number of digital twins for urban planning are obtained based on each industrial community; different industrial elements are divided into systems, so that the different industrial systems after the division are unrelated to each other, so that in urban planning, only all industrial elements under one industrial system need to be analyzed to avoid interference from other industrial elements. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 A flowchart of the steps of a method for constructing a digital twin for urban planning according to the present invention; Figure 2 It is a schematic diagram of matching the daily output of different industrial elements in one embodiment of the present invention. DETAILED DESCRIPTION

[0018] In order to further explain the technical means and effects adopted by the present invention to achieve the predetermined invention purpose, the following is a detailed description of the digital twin construction method for urban planning proposed by the present invention, its specific implementation method, structure, features and effects, in combination with the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" does not necessarily refer to the same embodiment. In addition, specific features, structures or characteristics in one or more embodiments may be combined in any suitable form.

[0019] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.

[0020] The following is a detailed description of a method for constructing a digital twin for urban planning provided by the present invention in conjunction with the accompanying drawings.

[0021] See also Figure 1 , which shows a flowchart of a method for constructing a digital twin for urban planning provided by an embodiment of the present invention, the method comprising the following steps: Step S001: Obtain the daily output of each industrial element in a preset date sequence.

[0022] It should be noted that the city contains multiple complete industrial system chains. The purpose of the mathematical twin constructed in this embodiment is to accurately divide the multiple industrial system chains, and then use each industrial system chain as the data basis of a mathematical twin, so that the constructed mathematical twins do not affect each other, and each mathematical twin only contains all industrial elements of all links in one industrial system chain.

[0023] It should be further explained that an industrial system chain contains multiple links, and the previous links will have an impact on the subsequent links. The main manifestation of the impact is reflected in the consistency of the output changes of each link in different industrial system chains. Therefore, this embodiment first collects the output of all links in the city, and regards different links of different industrial systems as an industrial element to obtain the daily output of each link.

[0024] Preferably, the specific steps of obtaining the daily output of each industrial element are: The length of the date sequence is preset, and the daily production quantity of each industrial element in the date sequence is obtained according to the public data of the city's industrial data statistics department. The production quantity is then normalized to the maximum and minimum values ​​to obtain the daily output of each industrial element. The set composed of all industrial elements is the city's industrial system set.

[0025] It should be noted that the date sequence length segment described in this embodiment is described by taking one year as an example, wherein the industrial element is the product produced by each factory in the city, and the production quantity is the total amount of the same product produced in the city; as an example, the industrial element of this embodiment is described by taking raw materials such as ore and agricultural products, semi-finished raw materials such as pig iron and kerosene after processing of raw materials, as well as power generation, industrial water output, sewage discharge, etc. as examples; It is particularly important to note that the public data of the city's industrial data statistics department in this embodiment is collected by taking the public data of the Industrial Statistics Division under the Municipal Statistics Bureau as an example. The Municipal Statistics Bureau organizes and implements national economic accounting and input-output surveys, and provides national economic accounting data after compiling the city's gross domestic product.

[0026] Step S002: Obtain the time lag consistency coefficient between each industrial element and other industrial elements according to the consistency of the output changes of each industrial element and other industrial elements and the difference in time lag when the outputs are consistent.

[0027] It should be noted that since industrial production needs to ensure that costs are as low as possible, placing the same industrial system in one city can minimize transportation costs. The conditions for the city where the industrial system is placed are usually the origin of raw materials or low energy consumption costs. Therefore, cities with raw material origins or low energy consumption costs can attract more industries. That is, one raw material or one energy consumption can affect multiple industrial chains. Therefore, when constructing a digital twin, because it is limited by the influence of the same upstream links, these multiple industrial chains need to be combined for analysis.

[0028] It should be further explained that different links in the same industrial system have different time effects, so changes in the upstream link will affect the next link after a certain period of time. Therefore, it is necessary to first conduct similarity analysis and time lag analysis on the current industrial elements and other industrial elements to obtain the time lag consistency coefficient between the current industrial elements and other industrial elements.

[0029] Preferably, according to the time lag difference between the current industrial element and other industrial elements, the specific steps of obtaining the time lag consistency coefficient between the current industrial element and each other industrial element are: It should be noted that, since the previous links in different links of an industrial system will affect the subsequent links, the entire complete industrial system is manifested as a tree structure as the links converge, and in the tree structure, the output changes of the industrial elements corresponding to the previous links will have an impact on the industrial elements corresponding to the subsequent links after a certain period of time, and the impact difference is small and the time of impact feedback is similar. Therefore, this embodiment matches the industrial elements through dynamic time warping, and then analyzes the time lag consistency and change similarity of different impacts, and obtains the time lag consistency coefficient of each industrial element with other industrial elements.

[0030] Specifically, a dynamic time warping algorithm is used to dynamically and regularly match the daily output of the ith industrial element with the jth industrial element except the ith industrial element, so as to obtain the output matching day of the ith industrial element on the jth industrial element every day, and the output matching day is the date with the shortest interval and the least output difference with the date of each day of the ith industrial element.

[0031] It should be noted that the dynamic time warping algorithm is a well-known technology and will not be described in detail in this embodiment. Figure 2 As shown, Figure 2 It is a schematic diagram of the matching of the daily output of different industrial elements, specifically, it is a schematic diagram of the daily output of the i-th industrial element and the j-th industrial element after being matched through the dynamic time warping algorithm.

[0032] Furthermore, if the date sequence value of the output matching day of the second day of the i-th industrial element is greater than or equal to 2, the j-th industrial element is recorded as a suspected downstream industrial element of the i-th industrial element.

[0033] It should be noted that the dynamic time warping algorithm will match elements with the same values ​​in two sequences. If the output is similar, they will be matched into a pair. If the i-th industrial element is the upstream industrial element and the j-th industrial element is the downstream industrial element, the i-th industrial element will have a hysteresis effect on the output of the j-th industrial element. Therefore, the dynamic time warping algorithm will use the first several days of the j-th industrial element as the first day of the i-th industrial element. The first day of the output matching day, at this time, it is considered that the j-th industrial element may be the downstream industrial element of the i-th industrial element.

[0034] Furthermore, all suspected downstream industrial elements of the i-th industrial element are obtained in the same way.

[0035] It should be noted that the above matching through the dynamic time warping algorithm has preliminarily obtained the suspected downstream industrial elements affected by the i-th industrial element, and has not taken into account the correlation strength of the influence. Therefore, it is necessary to conduct a consistency analysis of the influence time lag between each industrial element and its suspected downstream industrial element, and obtain the time lag consistency coefficient of each suspected downstream industrial element of the industrial element, and then screen the suspected downstream industrial elements to obtain the downstream industrial elements of each industrial element.

[0036] Preferably, according to the daily production matching relationship between the industrial element and its suspected downstream industrial element, the specific steps of obtaining the time lag consistency coefficient of each suspected downstream industrial element of the industrial element are: It should be noted that the amount of output ultimately required by the downstream of the industrial system is reflected in the upstream link, resulting in changes in the daily output of the upstream link, and the peak values ​​in the change can better reflect the degree of demand in the downstream link. Therefore, in order to highlight the time-lag relationship between the impact of industrial elements and their suspected downstream industrial elements, this embodiment first segments the date sequence of industrial elements according to the fluctuations in output data, and then matches the date sequences of industrial elements and their suspected downstream industrial elements based on the daily output matching days of industrial elements in suspected downstream industrial elements, thereby obtaining several output date segments of industrial elements and matching output date segments of their suspected downstream industrial elements.

[0037] Specifically, all minimum values ​​of the sequence formed by the daily output of the i-th industrial element are obtained, and the date sequence formed by all dates between two adjacent minimum values ​​and including the first minimum value of the two minimum values ​​is recorded as an output date segment of the i-th industrial element. Similarly, all output date segments of the i-th industrial element are obtained.

[0038] Further, for the kth production date segment of the ith industrial element, the first date in the kth production date segment, the production matching day in the date sequence of the mth suspected downstream industrial element of the ith industrial element, is recorded as the first day of the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element; the last date in the kth production date segment, the production matching day in the date sequence of the mth suspected downstream industrial element of the ith industrial element, is recorded as the last day of the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element. The sequence consisting of all dates between the first day and the last day of the kth matching production date segment in the mth suspected downstream industrial element is recorded as the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

[0039] Similarly, each production date segment of the i-th industrial element and the corresponding matching production date segment in each suspected downstream industrial element of the i-th industrial element are obtained.

[0040] It should be noted that if the mth suspected downstream industrial element is the downstream industrial element of the ith industrial element, then when the output of the ith industrial element changes, the time lag relationship of the change of the mth suspected downstream industrial element is consistent. The time lag is caused by the need for transportation, transit and other operations in the process of industrial elements from upstream links to downstream links. Therefore, this embodiment obtains the time lag consistency coefficient between each industrial element and its suspected downstream industrial element according to the consistency of the time difference between the output date segment and the corresponding matching output date segment.

[0041] Preferably, the time-lag consistency coefficient between the i-th industrial element and the m-th suspected downstream industrial element The calculation method is: in, represents the variance function; is an exponential function with a natural constant as base; is the median date of the k-th production date segment of the i-th industrial element, It is the median date of the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

[0042] It should be noted that the variance can reflect the size of the value, so this embodiment uses It reflects the consistency of the time difference between the production date segment and the corresponding matching production date segment in the suspected downstream industrial element. The more consistent the value is, the more consistent the time lag effect of the ith industrial element on the mth suspected downstream industrial element is, and the mth suspected downstream industrial element is more likely to be the downstream industrial element of the ith industrial element.

[0043] Similarly, the time-lagged consistency coefficient between each industrial element and other industrial elements is obtained.

[0044] Step S003: Filter the downstream industrial elements of each industrial element based on the time-lag consistency coefficient to obtain all the upstream industrial elements of each downstream industrial element; obtain the consistency of the output change of each upstream industrial element under each downstream industrial element according to the covariation relationship of the output between all the upstream industrial elements of each downstream industrial element when they are raw materials of the same link; obtain a number of industrial system elements according to the consistency of the output change of each upstream industrial element under the downstream industrial element in combination with the time-lag consistency coefficient of the corresponding industrial element.

[0045] A time lag consistency threshold is preset, and this embodiment is described by taking 0.4 as an example. When the time lag consistency coefficient between the ith industrial element and the mth suspected downstream industrial element is greater than or equal to the time lag consistency threshold, the mth suspected downstream industrial element of the ith industrial element is recorded as the downstream industrial element of the ith industrial element.

[0046] It should be noted that each link in the city's industrial system presents a tree structure, that is, the production of an industrial element may be affected by multiple industrial elements, which will lead to the existence of mutual correlation between multiple industrial elements that produce the same industrial element. Therefore, when the downstream industrial element changes, the changes in the upstream industrial elements are not obvious, or because the industrial element accounts for a small proportion in the generation of the next industrial element, it leads to a mismatch in output, which in turn leads to errors in the regularization due to the dissimilarity of changes during dynamic time regularization, which makes the accuracy of the time lag consistency coefficient low. Therefore, the above-mentioned time lag consistency threshold value is preset to be small.

[0047] It should be further explained that, since the proportion of industrial elements used as raw materials is different when producing the same industrial element, the time-lag consistency coefficient of the industrial elements with a larger proportion is more accurate than the numerical consistency coefficient of the industrial elements with a smaller proportion. Since the demand for industrial elements used as raw materials is similar, the change range of the output of industrial elements used as raw materials is similar. Therefore, this embodiment obtains the industrial system correlation coefficient of each downstream industrial element with the upstream industrial element based on the consistency of the output changes of the upstream industrial elements of the downstream industrial elements and the time-lag consistency coefficient.

[0048] Preferably, for the nth downstream industrial element of the ith industrial element, the ith industrial element is recorded as the nth downstream industrial element. upstream industrial element, the time-lag consistency coefficient between the ith industrial element and the nth downstream industrial element of the ith industrial element is recorded as the time-lag consistency coefficient between the nth downstream industrial element and the nth downstream industrial element. The time-lag consistency coefficient of the upstream industrial elements.

[0049] Similarly, the time-lag consistency coefficients of all upstream industrial elements of the nth downstream industrial element and between the nth downstream industrial element and all upstream industrial elements are obtained.

[0050] Furthermore, for the nth downstream industrial element The kth production date segment of the upstream industrial element, the The k-th production date segment of the n-th upstream industrial element is the same as the k-th production date segment of the ith upstream industrial element; among all the upstream industrial elements of the n-th downstream industrial element, except the The other upstream industrial elements except the first upstream industrial element are recorded as the first downstream industrial element. The covariant contribution of upstream industrial elements is the industrial element.

[0051] Further, in the nth downstream industrial element In the date sequence of each covariant contributing industrial element of the upstream industrial element, obtain the date sequence of the first The sequence segment with the same production date segment date as the kth production date segment of the nth downstream industrial element is recorded as the Each covariate contribution of the upstream industrial elements to the k-th production date segment under the industrial element.

[0052] Further, according to the nth downstream industrial element The k-th production date segment of the upstream industrial element and the n-th downstream industrial element The kth covariant contribution date segment of each covariant contribution industrial element under the nth upstream industrial element is combined with the nth downstream industrial element and the nth downstream industrial element. The time-lag consistency coefficient of the nth upstream industrial element is obtained, and the time-lag consistency coefficient of the nth downstream industrial element is obtained. The specific steps for calculating the industrial system correlation coefficient of each upstream industrial element are as follows: The nth downstream industrial element Consistency of production changes of upstream industrial elements The calculation method is: in, To find the function of Pearson correlation coefficient; The nth downstream industrial element The number of covariant contributing industrial elements of the upstream industrial elements, The nth downstream industrial element The number of production date segments of each upstream industrial element; The nth downstream industrial element The sequence of the production of the kth production date segment of the upstream industrial elements, The nth downstream industrial element The sequence of the covariant contributions of the upstream industrial elements and the production of the industrial elements in the k-th production date segment.

[0053] It is used to measure the consistency of the output changes of the upstream industrial element and its covariant contributing industrial element at the same time. The larger the value, the more likely it is that the element and its covariant contributing industrial element are both used as the raw materials of the downstream industrial element. The more upstream industrial elements there are, the more likely they are to form a link in an industrial system with the nth downstream industrial element.

[0054] Further, based on the nth downstream industrial element The consistency of the output change of the nth upstream industrial element and the consistency of the output change of the nth downstream industrial element The time-lag consistency coefficient of the upstream industrial element is corrected to obtain the nth downstream industrial element Industrial system correlation coefficient of upstream industrial elements The calculation method is: in, The relationship between the nth downstream industrial element and the nth downstream industrial element The time-lag consistency coefficient of the upstream industrial elements, For the nth downstream industrial element The consistency of output changes of upstream industrial elements.

[0055] Furthermore, the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element are obtained, and the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element are linearly normalized to obtain the normalized value of the industrial system correlation coefficient of each upstream industrial element under the nth downstream industrial element. The Otsu threshold algorithm is used to perform maximum difference segmentation on the normalized value of the industrial system correlation coefficient, and all upstream industrial elements in the largest part of the two after the segmentation are used as several upward-associated industrial elements of the nth downstream industrial element, and the industrial system correlation coefficients of the upstream industrial elements corresponding to each upward-associated industrial element of the nth downstream industrial element are used as the industrial system correlation coefficients of each upward-associated industrial element of the nth downstream industrial element; each downstream industrial element and each of its upward-associated industrial elements constitute an industrial system element.

[0056] It should be noted that this embodiment uses the maximum and minimum value normalization algorithm to achieve the normalization of the industrial system correlation coefficient; among which the Otsu threshold algorithm is an existing well-known technology and will not be described in detail in this embodiment.

[0057] Step S004: construct an urban industrial system network based on all industrial system elements, divide the urban industrial system network into communities to obtain several industrial communities, and obtain several digital twins for urban planning based on each industrial community.

[0058] It should be noted that after obtaining the industrial system elements composed of all downstream industrial elements and each of their upward associated industrial elements, an industrial system element represents the relationship between two adjacent links of an industrial system in the city, and all industrial systems in the city are connected due to this relationship. Therefore, by constructing an urban industrial system network based on a relational network structure for all industrial system elements, and then dividing the urban industrial system network, we can obtain several industrial communities, and then construct a digital twin for each industrial community to obtain several digital twins for urban planning.

[0059] Specifically, an urban industrial system network based on a relational network structure is constructed, wherein the urban industrial system network is a topological structure, and each node in the urban industrial system network represents the first industrial element. If two nodes can constitute an industrial system element, then the industrial system correlation coefficient between the downstream industrial element and the upwardly associated industrial element in the industrial system element is used as the degree of the urban industrial system network.

[0060] It should be noted that the industrial system referenced in this embodiment is a unidirectional industrial system. Therefore, the two nodes presenting the industrial system elements in the urban industrial system network point from the upstream link to the downstream link. Therefore, the degree of the node is unique and directional.

[0061] Furthermore, this embodiment divides the urban industrial system network into communities through the Louvain algorithm to obtain several industrial communities, and uses the historical output of industrial elements in each community on all dates as the data basis for generating a digital twin, thereby obtaining several digital twins for urban planning, wherein the industrial relationships between the digital twins do not interfere with each other.

[0062] It should be noted that the Louvain algorithm described in this embodiment is a prior art technology and will not be described in detail in this embodiment.

[0063] It should be noted that the The model only shows negative correlation and the output of the constraint model is in In the range, As the input of this model, it can be replaced by other models with the same purpose in specific implementation. The model is used as an example for description without any specific limitation.

[0064] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for constructing a digital twin for urban planning, characterized in that: The method comprises the following steps. The system comprises the following modules: Get the daily output of each industrial element in the preset date sequence; According to the consistency of the output change of each industrial element with other industrial elements and the difference of time lag in the output under the consistency condition, the time lag consistency coefficient of each industrial element with other industrial elements is obtained; Based on the time-lag consistency coefficient, the downstream industrial elements of each industrial element are screened to obtain all the upstream industrial elements of each downstream industrial element; According to the covariation relationship of the output between all the upstream industrial elements of each downstream industrial element as the raw materials of the same link, the consistency of the output change of each upstream industrial element under each downstream industrial element is obtained; According to the consistency of the output change of each upstream industrial element under the downstream industrial element, combined with the time-lag consistency coefficient of the corresponding industrial element, several industrial system elements are obtained; An urban industrial system network is constructed based on all industrial system elements. After dividing the urban industrial system network into communities, several industrial communities are obtained. Based on each industrial community, several digital twins for urban planning are obtained.

2. A method for constructing a digital twin for urban planning according to claim 1, characterized in that: The specific steps of obtaining the time-lag consistency coefficient include: For the i-th industrial element, obtain all production date segments of the i-th industrial element and all matching production date segments in the m-th suspected downstream industrial element of the i-th industrial element; The time-lag consistency coefficient between the i-th industrial element and the m-th suspected downstream industrial element The calculation method is: in, represents the variance function; is an exponential function with a natural constant as base; is the median date of the k-th production date segment of the i-th industrial element, It is the median date of the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

3. A method for constructing a digital twin for urban planning according to claim 2, characterized in that: The specific steps of obtaining the production date segment include: Use the dynamic time warping algorithm to dynamically match the daily output of the i-th industrial element with the j-th industrial element except the i-th industrial element, and obtain the daily output matching day of the i-th industrial element on the j-th industrial element; If the date sequence value of the second day's output matching day of the i-th industrial element is greater than or equal to 2, the j-th industrial element is recorded as the suspected downstream industrial element of the i-th industrial element; Obtain all the minimum values ​​of the sequence composed of the daily output of the ith industrial element, and record the date sequence composed of all dates between two adjacent minimum values ​​and including the first minimum value of the two minimum values ​​as a production date segment of the ith industrial element.

4. A method for constructing a digital twin for urban planning according to claim 3, characterized in that: The specific steps of obtaining the matching production date segment include: The first date in the k-th production date segment of the i-th industrial element, the production matching date in the date sequence of the m-th suspected downstream industrial element of the i-th industrial element, is recorded as the first day of the k-th matching production date segment in the m-th suspected downstream industrial element of the i-th industrial element; The last date in the k-th production date segment, the production matching date in the date sequence of the m-th suspected downstream industrial element of the i-th industrial element, is recorded as the last day of the k-th matching production date segment in the m-th suspected downstream industrial element of the i-th industrial element; The sequence consisting of all dates between the first day and the last day of the kth matching production date segment in the mth suspected downstream industrial element is recorded as the kth matching production date segment in the mth suspected downstream industrial element of the ith industrial element.

5. A method for constructing a digital twin for urban planning according to claim 1, characterized in that: The specific steps of obtaining the upstream industrial elements include: A time lag consistency threshold is preset, and when the time lag consistency coefficient between the ith industrial element and the mth suspected downstream industrial element is greater than or equal to the time lag consistency threshold, the mth suspected downstream industrial element of the ith industrial element is recorded as the downstream industrial element of the ith industrial element; For the nth downstream industrial element of the ith industrial element, the ith industrial element is recorded as the nth downstream industrial element. An upstream industrial element.

6. A method for constructing a digital twin for urban planning according to claim 1, characterized in that: The specific steps for obtaining the consistency of the output change include: The time-lag consistency coefficient between the ith industrial element and the nth downstream industrial element of the ith industrial element is recorded as the time-lag consistency coefficient between the nth downstream industrial element and the nth downstream industrial element. The time-lag consistency coefficient of the upstream industrial elements; In the nth downstream industrial element In the date sequence of each covariant contributing industrial element of the upstream industrial element, obtain the date sequence of the first The sequence segment with the same production date segment date as the kth production date segment of the nth downstream industrial element is recorded as the The k-th production date segment under each covariant contributing industrial element of the upstream industrial elements; The nth downstream industrial element Consistency of production changes of upstream industrial elements The calculation method is: in, To find the function of Pearson correlation coefficient; The nth downstream industrial element The number of covariant contributing industrial elements of the upstream industrial elements, The nth downstream industrial element The number of production date segments of each upstream industrial element; The nth downstream industrial element The sequence of the production of the kth production date segment of the upstream industrial elements, The nth downstream industrial element The sequence of the covariant contributions of the upstream industrial elements and the production of the industrial elements in the k-th production date segment.

7. A method for constructing a digital twin for urban planning according to claim 6, characterized in that: The specific steps of obtaining the covariant contribution industrial elements include: Among all the upstream industrial elements of the nth downstream industrial element, except the The other upstream industrial elements except the first upstream industrial element are recorded as the nth downstream industrial element. The covariant contribution of upstream industrial elements to industrial elements.

8. A method for constructing a digital twin for urban planning according to claim 1, characterized in that: The specific steps of obtaining the industrial system elements include: Obtain the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element, perform linear normalization on the industrial system correlation coefficients of all upstream industrial elements under the nth downstream industrial element to obtain the normalized value of the industrial system correlation coefficient of each upstream industrial element under the nth downstream industrial element, use the Otsu threshold algorithm to perform maximum difference segmentation on the normalized value of the industrial system correlation coefficient, and use all upstream industrial elements in the largest part of the two after the segmentation as several upwardly associated industrial elements of the nth downstream industrial element; The industrial system correlation coefficient of the upstream industrial element corresponding to each upwardly associated industrial element of the nth downstream industrial element is used as the industrial system correlation coefficient of each upwardly associated industrial element of the nth downstream industrial element; each downstream industrial element and each of its upwardly associated industrial elements constitute an industrial system element.

9. A method for constructing a digital twin for urban planning according to claim 8, characterized in that: The specific steps of obtaining the industrial system correlation coefficient include: The nth downstream industrial element Industrial system correlation coefficient of upstream industrial elements The calculation method is: in, The relationship between the nth downstream industrial element and the nth downstream industrial element The time-lag consistency coefficient of the upstream industrial elements, For the nth downstream industrial element The consistency of output changes of upstream industrial elements.

10. A method for constructing a digital twin for urban planning according to claim 1, characterized in that: The urban industrial system network is constructed according to all industrial system elements, and a number of industrial communities are obtained after the urban industrial system network is divided into communities. A number of digital twins for urban planning are obtained based on each industrial community, including: Constructing an urban industrial system network based on a relational network structure, wherein the urban industrial system network is a topological structure, wherein each node in the urban industrial system network represents the first industrial element, and if two nodes can constitute an industrial system element, the industrial system correlation coefficient between the downstream industrial element and the upwardly associated industrial element in the industrial system element is used as the degree of the urban industrial system network; The urban industrial system network is divided into communities through the Leuven algorithm to obtain several industrial communities. The output of industrial elements in each community on all dates is used as the data basis for generating a digital twin, thus obtaining several digital twins for urban planning.

Images