Method for Calculating Traffic Carbon Emissions of Light Passenger Vehicles Driven by Urban Land Use Functions

Abstract

In an embodiment of the present invention, a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions is provided, belonging to the field of computing technology, specifically including: Step 1, divide traffic zones and calculate the impact variables of each traffic zone according to the distribution and pattern of different urban land use functions; Step 2, calculate the number of passenger vehicle trips in each traffic zone according to the impact variables of each traffic zone; Step 3, calculate the traffic travel distribution between traffic zones according to the gravity model and the number of passenger vehicle trips; Step 4, calculate the urban traffic travel carbon emissions according to the traffic travel distribution; Step 5, perform multi-perspective analysis on the urban traffic travel carbon emissions according to preset rules to obtain the carbon emission regional-level allocation results. Through the solution of the present invention, the calculation efficiency, accuracy and adaptability are improved.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of computing technologies, and particularly to a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions. Background Art

[0002] At present, calculating traffic carbon emissions by estimating travel demand has become one of the important technical means for traffic carbon emission accounting, and the mixed-use development model provides an effective tool for connecting urban land use functions and travel demand. However, the existing methods still have deficiencies in depicting the travel interaction patterns between regions and localizing the calibration of model parameters. There is still a lack of a fine-scale traffic carbon emission calculation method that can accurately depict the travel interaction patterns between regions and calibrate parameters based on actual traffic behavior data. The main disadvantages of the existing technologies are as follows: The existing traffic carbon emission calculation method based on the mixed-use development model only considers the travel generation volume and average travel distance of a single region, ignoring the actual interaction patterns and travel flow distribution characteristics between regions, resulting in a significant deviation between the calculated result of the carbon emission spatial distribution and the actual situation; Most of the existing traffic carbon emission calculation methods analyze the carbon emission attribution from a single perspective, lacking a multi-perspective collaborative analysis framework of "departure place, arrival place, and moving process", making it difficult to comprehensively reflect the spatial attribution characteristics of traffic carbon emissions and unable to meet the diverse needs of different management entities for carbon emission responsibility allocation.

[0003] It can be seen that there is an urgent need for a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions with high calculation accuracy and adaptability. Summary of the Invention

[0004] In view of this, the embodiments of the present invention provide a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions, which at least partially solves the problem of poor calculation accuracy and adaptability in the existing technology.

[0005] The embodiments of the present invention provide a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions, including:

[0006] Step 1: Divide traffic zones according to the distribution and pattern of different urban land use functions, and calculate the influence variables of each traffic zone.

[0007] Step 2: Calculate the number of passenger vehicle trips of each traffic zone according to the influence variables of each traffic zone.

[0008] Step 3: Calculate the traffic travel distribution between each traffic zone according to the gravity model and the number of passenger vehicle trips.

[0009] Step 4: Calculate the urban traffic travel carbon emissions according to the traffic travel distribution.

[0010] Step 5: Analyze the carbon emissions of urban transportation from multiple perspectives according to preset rules to obtain the regional-level carbon emission allocation results.

[0011] According to a specific implementation manner of an embodiment of the present invention, step 1 specifically includes:

[0012] Step 1.1: After combining the road network, divide the urban area into multiple traffic zones according to the distribution and pattern of different land use functions in the city.

[0013] Step 1.2: Calculate the influence variables of each traffic zone, where the influence variables include land use, traffic accessibility, and socioeconomic factors.

[0014] According to a specific implementation manner of an embodiment of the present invention, step 2 specifically includes:

[0015] Step 2.1: Based on the influence variables, use the trip generation manual in ITE to calculate the total trip volume of each traffic zone.

[0016] Step 2.2: Decompose the total trip volume into internal trips and external trips and calculate the total external trip volume accordingly.

[0017] Step 2.3: Decompose the total external trip volume into walking trips, bicycle trips, public transportation trips, and vehicle trips, and calculate the probabilities of walking, cycling, and public transportation in external trips through a preset formula.

[0018] Step 2.4: Subtract the probabilities of walking, cycling, and public transportation from the total probability to obtain the trip probability of vehicles in external trips, and calculate the number of vehicle trips in each traffic zone according to the total external trip volume and the trip probability of vehicles.

[0019] According to a specific implementation manner of an embodiment of the present invention, step 3 specifically includes:

[0020] Step 3.1: Calculate the spatial interaction intensity between any two traffic zones of all traffic zones according to the number of vehicle trips, where the expression of the spatial interaction intensity is

[0021] ;

[0022] where represents the number of trips from traffic zone to traffic zone , represents the number of trips from traffic zone to traffic zone , and respectively represent traffic zones and traffic zones the number of passenger vehicle trips;

[0023] Step 3.2, assume that the Euclidean distance between traffic zone and traffic zone is , and the spatial interaction intensity is . Perform logarithmic transformation on the spatial interaction intensity and the Euclidean distance , and perform linear fitting on the transformed data, and calculate the distance decay exponent accordingly, where the expression of the linear fitting is

[0024] ;

[0025] Among them, represents the Euclidean distance between regions and , represents the spatial interaction intensity, is the constant term, represents the distance decay exponent;

[0026] Step 3.3, based on the distance decay exponent, use the gravity model formula to allocate the number of passenger vehicle trips of each traffic zone to different traffic zones to obtain the traffic travel distribution between traffic zones

[0027] .

[0028] According to a specific implementation manner of an embodiment of the present invention, the specific steps of step 4 include:

[0029] Based on the traffic travel distribution between traffic zones, use a vehicle emission calculation model to calculate the urban traffic travel carbon emissions according to the driving distance and the number of trips of passenger vehicles

[0030] ;

[0031] Among them, represents the carbon dioxide emissions of passenger vehicles, represents the travel distribution from traffic zone to traffic zone , represents the proportion of gasoline vehicles, represents the energy intensity of fuel, represents the carbon dioxide emission factor.

[0032] According to a specific implementation manner of an embodiment of the present invention, the preset rules include attributing the carbon emissions of urban traffic trips to the area where the starting point is located from the perspective of the departure place, attributing the carbon emissions of urban traffic trips to the area where the arrival point is located from the perspective of the arrival place, and, from the perspective of the moving process, dividing the carbon emissions of urban traffic trips proportionally according to the areas passed by the actual path of the vehicle.

[0033] The light-duty passenger vehicle traffic carbon emission measurement scheme based on urban land use function drive in the embodiment of the present invention includes: Step 1, dividing traffic zones and calculating the influence variables of each traffic zone according to the distribution and pattern of different urban land uses; Step 2, calculating the number of passenger vehicle trips in each traffic zone according to the influence variables of each traffic zone; Step 3, calculating the traffic travel distribution between traffic zones according to the gravity model and the number of passenger vehicle trips; Step 4, calculating the carbon emissions of urban traffic trips according to the traffic travel distribution; Step 5, performing multi-perspective analysis on the carbon emissions of urban traffic trips according to the preset rules to obtain the carbon emission regional-level distribution results.

[0034] The beneficial effects of the embodiment of the present invention are as follows: Through the solution of the present invention, the traffic travel flow distribution between regions driven by urban land use function is accurately measured. On the basis of obtaining the regional travel generation volume through the mixed-use development model, the model parameters are locally calibrated by combining the gravity model and actual traffic behavior data to accurately depict the traffic travel interaction mode between regions, thereby improving the accuracy and reliability of the calculation of the spatial distribution of traffic carbon emissions; establishing a multi-perspective collaborative traffic carbon emission attribution analysis method to spatially allocate carbon emissions from three perspectives: the starting point, the ending point, and the moving process of the trip, enhancing the practical value of the carbon emission measurement results in urban low-carbon planning and carbon emission reduction decision-making. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required to be used in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0036] Figure 1 It is a flow schematic diagram of a light-duty passenger vehicle traffic carbon emission measurement method based on urban land use function drive provided by an embodiment of the present invention;

[0037] Figure 2 It is a specific implementation flow schematic diagram of a light-duty passenger vehicle traffic carbon emission measurement method based on urban land use function drive provided by an embodiment of the present invention;

[0038] Figure 3 It is a diagram of the calculation result of the number of vehicle trips provided by an embodiment of the present invention;

[0039] Figure 4 A calculation result graph of inter - regional traffic travel distribution provided by an embodiment of the present invention;

[0040] Figure 5 A calculation result graph of traffic carbon emissions from the perspective of the departure place provided by an embodiment of the present invention;

[0041] Figure 6 A calculation result graph of traffic carbon emissions from the perspective of the arrival place provided by an embodiment of the present invention;

[0042] Figure 7 A calculation result graph of traffic carbon emissions from the perspective of the moving process provided by an embodiment of the present invention. Detailed implementation manners

[0043] The embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0044] The following uses specific specific examples to illustrate the implementation manners of the present invention. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. The present invention can also be implemented or applied through other different specific implementation manners. Various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that, without conflict, the following embodiments and the features in the embodiments can be combined with each other. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without making creative efforts belong to the scope of protection of the present invention.

[0045] It should be noted that the following describes various aspects of the embodiments within the scope of the appended claims. It should be obvious that the aspects described herein can be embodied in a wide variety of forms, and any specific structure and / or function described herein is illustrative only. Based on the present invention, those skilled in the art should understand that one aspect described herein can be implemented independently of any other aspect, and two or more of these aspects can be combined in various ways. For example, any number of aspects described herein can be used to implement the device and / or practice the method. In addition, this device can be implemented and this method can be practiced using other structures and / or functionality in addition to one or more of the aspects described herein.

[0046] It should also be noted that the diagrams provided in the following embodiments only illustrate the basic concept of the present invention in a schematic manner. The components shown in the diagrams only relate to those of the present invention, rather than being drawn according to the number, shape, and size of the components in actual implementation. The types, quantities, and proportions of the components in actual implementation can be arbitrarily changed, and the component layout type may also be more complex.

[0047] In addition, in the following description, specific details are provided to facilitate a thorough understanding of the examples. However, those skilled in the art will understand that the described aspects can be practiced without these specific details.

[0048] An embodiment of the present invention provides a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions, and the method can be applied to the process of calculating traffic carbon emissions in traffic management.

[0049] See Figure 1 , which is a schematic flowchart of a method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions provided by an embodiment of the present invention. As Figure 1 and Figure 2 shown, the method mainly includes the following steps:

[0050] Step 1: Divide traffic zones and calculate the influence variables of each traffic zone according to the distribution and pattern of different urban land use functions.

[0051] Further, the specific implementation of Step 1 includes:

[0052] Step 1.1: Combine the road network end and divide the urban area into multiple traffic zones according to the distribution and pattern of different urban land use functions.

[0053] Step 1.2: Calculate the influence variables of each traffic zone, where the influence variables include land use, traffic accessibility, and socioeconomic factors.

[0054] In specific implementation, traffic zones can be divided according to the distribution and pattern of different urban land use functions. Based on the existing traffic network, full consideration is given to the road network constraints to ensure that each zone can accurately reflect the internal traffic behavior characteristics. Each traffic zone consists of two or more urban land use functions, and the sizes of the traffic zones are relatively consistent. By restricting the area of the region to less than 3 square kilometers, it is ensured that within each traffic zone, walking or cycling is the main mode of travel, while between traffic zones, passenger cars are the main mode of travel, thus facilitating the accurate estimation of the traffic volume between regions. At the same time, calculate the influence variables of each traffic zone. The influence variables include three types of parameters: land use, traffic accessibility, and socioeconomic factors, and the values of these indicators are determined by the spatial planning policies under the overall framework.

[0055] Step 2: Calculate the number of passenger vehicle trips for each traffic zone based on the influencing variables of each traffic zone.

[0056] Based on the above embodiments, step 2 specifically includes:

[0057] Step 2.1: Based on the influencing variables, use the trip rates in the ITE Trip Generation Manual to calculate the total trip volume for each traffic zone.

[0058] Step 2.2: Decompose the total trip volume into internal trips and external trips and calculate the total external trip volume accordingly.

[0059] Step 2.3: Decompose the total external trip volume into walking trips, bicycle trips, public transportation trips, and passenger vehicle trips, and calculate the probabilities of walking, cycling, and public transportation in external trips through a preset formula.

[0060] Step 2.4: Subtract the probabilities of walking, cycling, and public transportation from the total probability to obtain the trip probability of passenger vehicles in external trips. Calculate the number of passenger vehicle trips for each traffic zone based on the total external trip volume and the trip probability of passenger vehicles.

[0061] In specific implementation, the specific process of calculating the number of passenger vehicle trips for each traffic zone can be as follows:

[0062] 2.1 Estimate the total vehicle trip volume:

[0063] Adopt the Institute of Transportation Engineers (ITE) model based on land use, and use the trip rates in the ITE Trip Generation Manual to estimate the total number of trips. The basic formula is as shown in formula (1) below:

[0064] ;

[0065] In the formula, is the total vehicle trip volume of zone , is the area of urban land use function in zone ; The function calculates the total trip volume of urban land use type starting from traffic zone . Its specific calculation formula refers to the introduction of the ITE model. This process estimates the trip volume of each land use type through the influencing degree of different land use function types on the trip volume.

[0066] 2.2 Decompose vehicle trips into internal and external vehicle trips

[0067] ①The vehicle trip types are divided into three categories according to their purposes: home-work trips, home-other trips, and non-home trips. The total number of trips is decomposed into the total number of trips for each travel purpose based on the NCHRP report. ②The estimated vehicle trips are decomposed into internal and external vehicle trips: the logarithmic odds ratio of internal trips is calculated through formula (2). ③By exponentiating the logarithmic odds ratio, the ratio of internal trips is obtained and converted into a probability, as shown in formula (3). ④The total external trips are obtained by subtracting the internal trips from the total trips estimated according to the ITE Trip Generation Manual, as shown in formula (4).

[0068] ;

[0069] ;

[0070] ;

[0071] In the formula, represents the internal trip odds ratio from traffic zone ; represents the internal trip probability from traffic zone ; represents the total external trips from traffic zone ; is the influencing variable; , are the correlation coefficients based on existing research. In existing research, the estimation of this coefficient is mainly obtained through regression of a large amount of real survey data. Since the survey data comes from cities in different regions and with different levels of development around the world, the reliability of this model is widely recognized and it is widely used in various cities around the world.

[0072] 2.3 Modal split of external trips

[0073] The estimated external trips are decomposed into walking trips, bicycle trips, public transport trips, and vehicle trips. The probabilities of walking, cycling, and public transport in external trips are calculated through the equations in the model. The formulas are as shown in (5) to (10):

[0074] ;

[0075] ;

[0076] ;

[0077] ;

[0078] ;

[0079] ;

[0080] In the formula, , , represent the occurrence ratios of walking trips, bicycle trips, and public transportation trips for external trips from traffic zone ; , , represent the probabilities of walking trips, bicycle trips, and public transportation trips for external trips from traffic zone ; is an influencing variable; , , , , , The calculation method is the same as above.

[0081] 2.4 Calculation of External Trips of Light Passenger Vehicles

[0082] The probability of trips by light passenger vehicles in external trips is obtained by subtracting the probabilities of walking, cycling, and public transportation from the total probability. The formula is shown in (11). The total external trip volume is multiplied by this probability to obtain the total external trip volume generated by light passenger vehicles. The formula is shown in (12).

[0083] ;

[0084] ;

[0085] In the formula, represents the occurrence ratio of light passenger vehicles in external trip distances; represents the number of trips by light passenger vehicles.

[0086] Step 3: Calculate the traffic trip distribution between each traffic zone according to the gravity model and the number of trips by passenger vehicles.

[0087] Based on the above embodiments, step 3 specifically includes:

[0088] Step 3.1: Calculate the spatial interaction intensity between every two traffic zones of all traffic zones according to the number of trips by passenger vehicles. Among them, the expression of the spatial interaction intensity is

[0089] ;

[0090] Among them, represents the number of trips from traffic zone to traffic zone , Indicates the number of trips from traffic zone to traffic zone . The number of trips and respectively represent the number of trips made by passenger vehicles in traffic zone and traffic zone ;

[0091] Step 3.2, assume that the Euclidean distance between traffic zone and traffic zone is , and the spatial interaction intensity is . Perform logarithmic transformation on the spatial interaction intensity and the Euclidean distance , and perform linear fitting on the transformed data, and calculate the distance decay exponent based on this. Among them, the expression of the linear fitting is

[0092] ;

[0093] Among them, represents the Euclidean distance between regions and , represents the spatial interaction intensity, is the constant term, represents the distance decay exponent;

[0094] Step 3.3, based on the distance decay exponent, use the gravity model formula to allocate the number of trips made by passenger vehicles in each traffic zone to different traffic zones, and obtain the traffic travel distribution between traffic zones

[0095] .

[0096] In specific implementation, the travel assignment method based on the gravity model is mainly used to calculate the travel number assignment between regions. The gravity model assumes that the travel demand is proportional to the attractions of the departure region and the arrival region, and inversely proportional to a certain power of the travel distance. By applying the gravity model, the travel demand assignment between regions can be estimated more accurately, thereby providing basic data for further calculation of traffic travel carbon emissions. This process can be decomposed into the following three steps:

[0097] 3.1 Calculation of spatial interaction between regions

[0098] Spatial interaction refers to the interaction intensity between two regions. Usually, the spatial interaction between regions can be estimated through actual taxi trajectory data. The specific method is to use taxi trajectory data to analyze the travel frequency between different regions, so as to quantify the interaction intensity between regions. Suppose there are two regions and , the interaction intensity The calculation formula is as shown in Equation (13):

[0099] ;

[0100] In the formula, represents the number of trips from area to area . represents the number of trips from area to area . and respectively represent the total number of trips in area and area . By calculating the interaction intensity between each area through this step, it provides basic data for the subsequent calculation of the gravity model.

[0101] 3.2 Calculation of distance decay exponent

[0102] The distance decay exponent is a measure of the rate of change of spatial interaction intensity with geographical distance. Specifically, the distance decay exponent can be calculated by linearly fitting the spatial interaction intensity and distance between different areas in taxi trajectory data. Assume that the Euclidean distance between two areas and is , and the spatial interaction intensity is , then the distance decay exponent can be calculated through the following steps: First, perform logarithmic transformation on and , and then perform linear fitting on the transformed data. The fitting equation is as shown in Formula (14):

[0103] ;

[0104] In the formula, is the Euclidean distance between area and ; is the spatial interaction intensity; is the constant term; is the distance decay exponent. The value obtained by fitting can be used for the calculation of the gravity model.

[0105] 3.3 Trip number allocation

[0106] Using the gravity model formula, allocate the number of trips to each area. The specific formula is as shown in Equation (15):

[0107] ;

[0108] In the formula, is the number of trips from zone to zone ; is the total number of trips in zone ; is the Euclidean distance from the center of zone to the center of zone ; is the distance decay exponent.

[0109] Step 4: Calculate the urban traffic travel carbon emissions according to the traffic travel distribution.

[0110] Furthermore, the specific steps of Step 4 include:

[0111] Based on the traffic travel distribution among traffic zones, use the vehicle emission calculation model to calculate the urban traffic travel carbon emissions according to the driving distance and number of trips of passenger vehicles.

[0112] ;

[0113] Among them, represents the carbon dioxide emissions of passenger vehicles, represents the travel distribution from traffic zone to traffic zone ; represents the proportion of gasoline vehicles, represents the energy intensity of fuel, represents the carbon dioxide emission factor.

[0114] In specific implementation, based on the estimation result of the spatial distribution of traffic travel volume according to urban land use functions, use the vehicle emission calculation model to calculate the carbon dioxide emissions through the driving distance and number of trips of light passenger vehicles. The formula for the total carbon dioxide emissions of light passenger vehicles is shown in Equation (16):

[0115] ;

[0116] In the formula, is the carbon dioxide emissions of light passenger vehicles; is the number of trips from zone to zone ; is the Euclidean distance from zone to zone ; is the proportion of gasoline vehicles; is the energy intensity of fuel; is the carbon dioxide emission factor.

[0117] Step 5: Analyze the carbon emissions of urban traffic trips from multiple perspectives according to preset rules to obtain the regional-level distribution results of carbon emissions.

[0118] Further, the preset rules include attributing the carbon emissions of urban traffic trips to the region where the starting point is located from the perspective of the starting point, attributing the carbon emissions of urban traffic trips to the region where the arrival point is located from the perspective of the arrival point, and, from the perspective of the moving process, dividing the carbon emissions of urban traffic trips proportionally according to the regions passed by the actual path of the vehicle.

[0119] In specific implementation, when analyzing the carbon emission attribution region, existing research mainly focuses on the road network. Its advantage is that it can more accurately reflect the actual driving path of the vehicle and its corresponding carbon emissions. However, since China's management and planning are mainly based on administrative units, such as blocks and streets, attributing carbon emissions to the road network is not convenient for management and planning. In addition, the regional summary of traffic carbon emissions can help identify high-emission regions, so as to concentrate resources for optimization and improvement and improve the efficiency of carbon emission management. Therefore, the method of the present disclosure is based on the traffic zone scale and considers three perspectives to analyze the carbon emissions of traffic trips: the starting point perspective, the arrival point perspective, and the moving process perspective. The starting point perspective attributes the generated carbon emissions of traffic trips to the region where the starting point is located; the arrival point perspective attributes the carbon emissions to the region where the arrival point is located; the moving process perspective divides the carbon emissions proportionally according to the regions passed by the actual path of the vehicle. The following will specifically explain and analyze the three perspectives:

[0120] ① Starting point perspective: Under this perspective, the carbon emissions of traffic trips are attributed to the region where the starting point is located. This perspective believes that regardless of the final arrival region of the vehicle, all emissions are considered to be contributed by the starting region. This perspective helps to understand the impact of traffic behavior starting from a certain region on carbon emissions and can effectively reflect the contribution of traffic behavior in each region to carbon emissions.

[0121] ② Arrival point perspective: Contrary to the starting point perspective, the arrival point perspective attributes the carbon emissions of traffic trips to the region where the arrival point is located. That is, regardless of where the vehicle starts, all emissions are attributed to the region where the vehicle finally arrives. This perspective helps to analyze the traffic carbon emission burden borne by a certain region when it is the destination. This attribution method is not only convenient for management and decision-making but also helps to formulate more targeted carbon emission reduction policies.

[0122] ③ Moving process perspective: This perspective is more refined and divides the carbon emissions proportionally according to the regions passed by the actual driving path of the vehicle on the road network. That is, the carbon emissions generated during the vehicle's driving in each region are allocated to the corresponding region. This perspective can more comprehensively reflect the carbon emission situation of each region during the traffic trip and provide more accurate carbon emission data.

[0123] The method for calculating carbon emissions from light passenger vehicle transportation based on urban land function drive provided in this embodiment, by combining a gravity model based on actual traffic behavior data such as taxi trajectory calibration, accurately depicts the interaction pattern of traffic travel between regions, achieves high-precision calculation of the spatial distribution of traffic carbon emissions, and improves the accuracy of carbon emission spatial attribution; it can achieve fine-scale transportation carbon emissions calculation based only on easily accessible urban land function distribution data and some travel behavior sample data without relying on a large amount of field survey data, thereby reducing calculation costs and improving calculation efficiency.

[0124] The method of the present invention will be further described below in conjunction with a specific embodiment, and the specific implementation process of the present invention will be described using a case study:

[0125] (1) In this example, Area A is selected as the research area. The area is about 78.66 square kilometers, and the permanent population is about 1.521 million by the end of 2023. As a political, economic, and cultural center, the area has rich traffic variables and diverse land use functions, making it suitable for the study of transportation carbon emissions.

[0126] (2) Division of research units

[0127] According to the land use distribution characteristics and traffic network structure of Area A, the study area is divided into 52 traffic communities, each with an area of 1.5-3 square kilometers and an average area of 2.2 square kilometers. The division process fully considers the separation effect of administrative boundaries, natural geographical boundaries and major traffic arteries to ensure that the land use function within each traffic community is relatively homogeneous and that light passenger vehicles are the main means of travel between traffic communities.

[0128] (3) Calculation of influencing variables

[0129] The following influencing variables are calculated for each traffic zone:

[0130] ① Land use mix: This factor is used to measure the diversity of land use types within the transportation zone. The calculation formula is:

[0131] ;

[0132] In the formula, Indicates the traffic zone The proportion of the coverage area of this type of land use type to the total area of the transportation zone.

[0133] ② Activity density: This factor is used to measure the density of crowd activities in a traffic zone. The calculation formula is:

[0134] ;

[0135] In the formula, Indicates traffic area The number of permanent residents in Indicates traffic area The number of employed people in Indicates traffic area area.

[0136] ③Intersection density: This factor is used to measure the density of road intersections in a traffic zone. The calculation formula is:

[0137] ;

[0138] In the formula, Indicates traffic area The number of permanent residents in Indicates traffic area The number of employed people in Indicates traffic area area.

[0139] ④ Ratio of employed population to permanent population: This factor is used to measure the ability of a transportation community to meet the travel needs of residents for work commuting. The value range is [0, 1]. When the value is 1, it means that the transportation community can meet the travel needs of all residents for work commuting.

[0140] ⑤ Accessibility index: calculate the average distance and accessibility of each transportation community to major city centers and public transportation stations.

[0141] (4) Calculation of vehicle trip times

[0142] Based on the above influencing variables and equations 1-12 mentioned above, the internal travel ratio, walking travel ratio, bicycle travel ratio, public transportation travel ratio and light passenger vehicle travel ratio of each traffic zone are calculated, and then the travel volume generated by each traffic zone is calculated. The travel volume calculation results of each traffic zone in 2018 are as follows: Figure 3 As shown in the figure, the overall travel volume in the urban area of Area A shows a trend of uniform distribution, with higher travel volumes in the northeast and southwest corners. The number of trips in traffic areas A and B in the northeast corner of the urban area of Area A and traffic area C in the south is significantly higher than that in other areas.

[0143] (5) Calculation of inter-regional travel flow distribution

[0144] Using the taxi trajectory data in September 2018, the origin and destination of each trip were extracted and matched to traffic zones to construct an actually observed OD matrix. Based on this matrix, the intensity of spatial interaction between regions was calculated, and the distance decay exponent was determined through linear fitting. Substituting the decay exponent into the gravity model and combining the total trips of light passenger vehicles in each traffic zone calculated in step (4), the distribution of light passenger vehicle trips between regions was calculated. The calculation results are as Figure 4 shown. The results show that the intensity of trips between regions exhibits obvious spatial heterogeneity, mainly concentrated among the central business district of Region A, Commercial Area B, and the surrounding areas of the bonded area in Region A, which is highly consistent with the actually observed travel patterns.

[0145] (6) Calculation and multi-perspective analysis of traffic carbon emissions

[0146] According to the distribution of trips between regions and the travel distance, the total driving distance of light passenger vehicles was calculated. Based on the current development status of new energy vehicles in Region A, it was set that the proportion of gasoline vehicles in light passenger vehicles was 80% and the proportion of new energy vehicles was 20%. Referring to the energy consumption parameters in relevant research literature, the energy intensity of gasoline vehicles was 0.08 kgce / km, and the CO2 emission factor was 2.08 kgCO2 / kgce; for new energy vehicles, the direct emissions were zero (indirect emissions were not considered). Based on Formula 18, the carbon dioxide emissions of light passenger vehicles generated in each traffic zone were calculated, and attribution analysis was carried out from three perspectives: the perspective of the departure place, the perspective of the arrival place, and the perspective of the moving process. The calculation results are as Figure 5 、 Figure 6 、 Figure 7 shown.

[0147] It should be understood that each part of the present invention can be implemented by hardware, software, firmware, or a combination thereof.

[0148] The above is only the specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims

1. A method for calculating traffic carbon emissions of light passenger vehicles driven by urban land use functions, characterized in that, Including: Step 1: According to the distribution and pattern of different land use functions in the city, divide traffic zones and calculate the impact variables of each traffic zone; Step 2: Calculate the number of passenger vehicle trips of each traffic zone according to the impact variables of each traffic zone; Step 3: Calculate the traffic trip distribution between traffic zones according to the gravity model and the number of passenger vehicle trips; The specific content of Step 3 includes: Step 3.1: Calculate the spatial interaction intensity between every two traffic zones of all traffic zones according to the number of passenger vehicle trips. Among them, the expression of the spatial interaction intensity is Among them, represents the number of trips from traffic zone to traffic zone ; represents the number of trips from traffic zone to traffic zone ; and respectively represent the number of trips made by passenger vehicles in traffic zone and traffic zone . Step 3.2, assume traffic zones and traffic zone The Euclidean distance between them is , the spatial interaction intensity is , for the spatial interaction intensity and the Euclidean distance Perform logarithmic transformation, and perform linear fitting on the transformed data, and calculate the distance decay exponent accordingly , where the expression of the linear fitting is Among them, represents the Euclidean distance between and represents the spatial interaction intensity, is a constant term, represents the distance decay exponent; Step 3.3: Based on the distance decay index, use the gravity model formula to distribute the number of passenger vehicle trips of each traffic zone to different traffic zones to obtain the traffic trip distribution between traffic zones ； Step 4: Calculate the carbon emissions of urban traffic trips according to the traffic trip distribution; Step 5: Analyze the carbon emissions of urban traffic trips from multiple perspectives according to preset rules to obtain the carbon emission regional-level allocation results.

2. The method according to claim 1, characterized in that The specific content of Step 1 includes: Step 1.1: Combine with the road network end, and divide the urban area into multiple traffic zones according to the distribution and pattern of different land use functions in the city; Step 1.2: Calculate the impact variables of each traffic zone. Among them, the impact variables include land use, traffic accessibility and socioeconomic factors.

3. The method according to claim 2, wherein The specific content of Step 2 includes: Step 2.1: Based on the impact variables, use the trip rate in the ITE Trip Generation Manual to calculate the total trip volume of each traffic zone; Step 2.2: Decompose the total trip volume into internal trips and external trips and calculate the total external trip volume accordingly; Step 2.3: Decompose the total external trip volume into walking trips, bicycle trips, public transportation trips and passenger vehicle trips, and calculate the probabilities of walking, bicycle and public transportation in external trips through a preset formula; Step 2.4: Subtract the probabilities of walking, bicycle and public transportation from the total probability to obtain the trip probability of passenger vehicles in external trips, and calculate the number of passenger vehicle trips of each traffic zone according to the total external trip volume and the trip probability of passenger vehicles.

4. The method according to claim 3, wherein The specific content of Step 4 includes: Based on the traffic trip distribution between traffic zones, use the vehicle emission calculation model to calculate the carbon emissions of urban traffic trips according to the driving distance and the number of trips of passenger vehicles wherein, represents the carbon dioxide emissions of a passenger vehicle, represents a traffic zone to traffic zone for trip distribution, represents the proportion of gasoline vehicles, represents the energy intensity of fuel, represents the carbon dioxide emission factor.

5. The method according to claim 4, wherein The preset rules include attributing the carbon emissions of urban traffic trips to the region where the starting point is located from the perspective of the starting point, attributing the carbon emissions of urban traffic trips to the region where the arrival point is located from the perspective of the arrival point, and dividing the carbon emissions of urban traffic trips proportionally according to the regions actually passed by the vehicle's actual path from the perspective of the moving process.

Images