Driving risk assessment method for automatic driving vehicle in expressway interleaving area

Abstract

The invention belongs to the technical field of traffic control systems, and particularly relates to an expressway interleaving area automatic driving vehicle driving risk assessment method, which comprises the steps of 1, obstacle space-time risk field modeling, 2, road boundary and road line risk field modeling, 3, expressway upper and lower ramp risk field modeling, 4, total risk field calculation, 5, expressway interleaving area automatic driving vehicle driving risk assessment, 6, expressway interleaving area automatic driving vehicle driving risk assessment, 7, expressway interleaving area automatic driving vehicle driving risk assessment and 7, expressway interleaving area automatic driving vehicle driving risk assessment. The method has the advantages that the physical quantity of the space-time distance is innovatively introduced, the dimension of a two-dimensional space risk field is raised to a three-dimensional time-space risk field, the real-world expressway vehicle track is extracted, and the real-world expressway vehicle track is extracted. And the influence of the prediction trajectory and the time effect on the risk is effectively quantified. A special geometric configuration field of the upper ramp and the lower ramp of the interlacing area is also constructed, and accurate limitation is applied to operation of a specific vehicle. And on the basis of a risk dynamic balance theory, YOLO and other machine vision technologies, model parameters are calibrated by using real aerial photography data.

Description

Technical Field

[0001] The present invention belongs to the technical field of traffic control systems, and specifically relates to a driving risk assessment method for autonomous driving vehicles in an expressway weaving area based on a spatiotemporal risk field. Background Art

[0002] Expressways, as the core arteries of urban transportation networks, bear a heavy traffic load, but traffic safety issues are particularly acute in weaving zones. Traffic flow in these zones exhibits a typical three-way interaction: on-ramps merge, off-ramps diverge, and vehicles on the main line travel through, creating multiple weaving flows. Frequent lane changes induce turbulence, significantly increasing accident rates and congestion indices, making them bottlenecks with high accident rates and major congestion. In recent years, the rapid development of automated driving (CAV) technology has provided a new technological path to addressing driving risks in weaving zones with its precise perception, decision-making, and control capabilities. However, in the current traffic environment where CAVs and human-driven vehicles (HDVs) coexist for a long time, the external complexity of driving risks has increased significantly. On the one hand, the inherent uncertainty of traditional traffic flow persists; on the other hand, the significant differences in driving behavior between CAVs and HDVs further exacerbate operating risks in weaving zones. This dual risk effect has, to a certain extent, limited the potential advantages of CAV technology. Therefore, conducting a systematic assessment and quantitative calibration study on the driving risks of autonomous vehicles in expressway weaving areas has important theoretical value and practical significance for ensuring the accurate operation of CAV lane-changing decisions and trajectory planning.

[0003] Common driving risk assessment methods currently fall into three main categories: time-based (such as time to collision (TTC), headway (THW), and time to approach (PET), distance-based (such as minimum safe distance (MSD)), and acceleration-based (such as deceleration rate (DRAC)). A Chinese patent application, "A Method for Driving Risk Assessment in Underground Tunnels Based on Dynamic Bayesian Networks," patent application number CN202411687491.1, proposes using a Bayesian network to calculate safe distances, thereby assessing rear-end collision risk. Another Chinese patent application, "A Dynamic Driving Risk Assessment Method and System Based on DNC and NGBoost-PA-WOE," patent application number CN202411663328.1, proposes using machine learning methods to learn from traffic accident-related data (such as time, distance, acceleration, etc.) to predict accident risk. These methods are computationally simple and efficient, and are widely used in risk identification and early warning systems. However, there are still significant deficiencies in practical applications: (1) The three types of methods generally assume that the vehicle state remains unchanged in a short period of time, which makes it difficult to reflect the uncertainty of vehicle movement and its surrounding environment, and easily leads to risk misjudgment; (2) Most methods focus on the identification of longitudinal rear-end collision risks, and it is difficult to effectively assess the more critical lateral collision risk during lane changing; (3) In scenarios with complex risk sources and diverse interactions, traditional methods tend to underestimate risks, which limits their application effect in highly complex traffic environments.

[0004] Risk field theory, an improved artificial potential field approach, originates from field theory in physics. Compared to traditional risk assessment methods, it offers advantages such as strong dynamic modeling capabilities, the ability to quantify risk distribution, and adaptability to the characteristics of diverse traffic participants. It can dynamically reflect the interactive risks between traffic participants across time and space, making it particularly suitable for assessing and predicting risk dynamics in complex, mixed traffic scenarios. For example, the article "Extension and Application of Driving Risk Field Models in a Collaborative Vehicle-Road Environment" proposes a risk assessment model for pure CAV scenarios that incorporates parameters such as TTC, vehicle geometry, and yaw angle. The Chinese patent application "Automated Driving Takeover Risk Assessment Model Considering Driving Risk Field and Its Construction Method" (patent application number: CN202211137773.5) proposes a method for constructing an automated driving takeover risk assessment model that considers driving risk fields, constructing a distance-based static risk field. The Chinese patent application "Vehicle Driving Risk Identification Method Based on Digital Twins and Predicted Risk Fields" (patent application number: CN202411018302.1) uses the predicted risk field to superimpose a series of potential static fields within the prediction time domain to obtain a risk field for the corresponding single-modal predicted trajectory, and uses this to identify driving risks. Although increasing research has improved risk fields, these constructed risk fields remain essentially two-dimensional static fields, accounting only for the risk impact of current risk sources and ignoring the impact of the future motion trends of dynamic obstacles on the current risk. Furthermore, the special fields in interweaving scenarios have rarely been studied. Due to the increased frequency of lane changes and the dominance of turbulence, the risk field in weaving areas is particularly unique, and the impact of the environment on risk is far greater than that of ordinary roads. However, there is still no relevant research to systematically explain this change. Summary of the Invention

[0005] In view of the above problems, the purpose of the present invention is to provide a driving risk assessment method for autonomous driving vehicles in the weaving area of expressways. By introducing the physical quantity of space-time distance, the two-dimensional spatial risk field is upgraded to a three-dimensional time-space risk field, effectively quantifying the impact of predicted trajectory and time effect on risk, and constructing a special geometric configuration field of the on- and off-ramps in the weaving area, which imposes accurate restrictions on the operation of specific vehicles.

[0006] The present invention provides a method for assessing driving risks of an autonomous driving vehicle in an expressway weaving zone, comprising the following steps: Step 1: Modeling the space-time risk field of obstacles, including: constructing a space-time coordinate system, calculating space-time distance, calculating the charge considering its own and motion properties, calculating the dielectric constant considering anisotropy and environmental complexity, and calculating the space-time risk of obstacles; Step 2: Modeling the risk field of road boundaries and road lines, including: calculating the field strength of road boundaries, calculating the field strength of road solid lines, calculating the field strength of road dashed lines, and calculating the risk field of road boundaries and road lines; Step 3: Modeling the risk field of on- and off-ramps on the expressway; Step 4: Calculate the total risk field; Step 5: Real-world freeway vehicle trajectory extraction, including: vehicle trajectory extraction based on YOLOv8, target tracking module based on Deep-SORT, video stabilization processing, model training and testing, and coordinate system conversion; Step 6: Parameter calibration based on the risk dynamic balance theory, including: determining the objective function, determining the upper and lower risk limits of the risk dynamic balance theory, and using a hybrid genetic algorithm and a variable neighborhood search algorithm to solve the problem and realize the driving risk assessment of autonomous vehicles in the expressway weaving area.

[0007] As a preferred embodiment of the present invention, step one also includes the following steps: A1. Construct a space-time coordinate system. A spatiotemporal three-dimensional coordinate system is established based on the Frenet coordinate system, and the direction of travel along the road is set as Axis, perpendicular to the road driving direction Axis, time is set to axis; A2. Calculate spatiotemporal distance, including: obtaining predicted obstacle trajectories, calculating spatial distances that take into account vehicle geometry, calculating equivalent time-based distances that take into account risk anisotropy, and obtaining spatiotemporal distances. A3. Calculate the amount of charge taking into account its own and motion properties. Two representative parameters are selected from the comprehensive attributes of vehicle functions and driving behavior as key influencing factors of charge: braking performance and reaction time; Charge The calculation formula is defined as follows; ; Where, is the mass of the obstacle, is the speed of the vehicle, is the maximum deceleration of the vehicle, is the vehicle reaction time, which is the sum of the driver's reaction time and the vehicle system delay, and is the unknown coefficient, is an exponential function; A4. Calculate the dielectric constant taking into account anisotropy and environmental complexity. The dielectric constant in the risk field is defined as the product of the anisotropic attenuation coefficient and the environmental attenuation coefficient. The calculation formula is as follows: ; in, is the anisotropic attenuation coefficient, is the environmental complexity, which is calculated as follows: ; ; in, The obstacle's direction of travel and The angle between is the position vector of the obstacle pointing to the observation point, and are the speed and acceleration of the vehicle, is the ratio of the square difference of the speed in the weaving section to that of the normal road, is the ratio of the lane change ratio in the weaving section to that in the normal road, 、 and is the coefficient to be determined; A5. Calculate the spatial and temporal risk of obstacles. Obstacle space-time risk The calculation formula is as follows: ; in, is the space-time distance; is the charge; is the dielectric constant.

[0008] As a preferred embodiment of the present invention, step A2 further includes the following steps: A2.1. Obtain the predicted trajectory of the obstacle. Obstacle trajectory scatter point collection Expressed as: ,in, express The space-time trajectory of the obstacle prediction at each moment, , , They are The position of the obstacle along the direction of travel of the road and the position perpendicular to the direction of travel of the road; A2.2. Calculate the spatial distance considering the vehicle's geometric outline. In the vehicle body coordinate system, the center of mass of the vehicle is taken as the origin and the direction of the vehicle is Axis, perpendicular to the direction of vehicle travel Axis, where, in the space-time coordinate system, point The coordinates are , the coordinates of the obstacle mass center are , in the vehicle body coordinate system, the four vertices of the rectangle The coordinates are: ,In addition, the vehicle's yaw angle , the speed of the vehicle , two types of coordinate systems are transformed using translation and rotation, as follows: ; Solving spatial distance based on vehicle body coordinate system The formula is as follows; ; ; Where, is the coordinate of the observation point in the vehicle body coordinate system; is the coordinate of the observation point; is the transformed coordinate; is the vehicle width; is the vehicle length, The four vertices of the rectangle With the transformed observation point spatial distance; A2.3. Calculate the equivalent time-base distance considering risk anisotropy. Introducing equivalent time base distance , equivalent time base distance Indicates: When the point of the interactive object Directly in front of the obstacle course, the point of the interactive object When not directly in front of the vehicle's path, assume that the point of the interactive object Points with interactive objects With equal collision risk, the point from the vehicle to the interactive object The time distance is equal to the point from the vehicle to the interactive object Equivalent time distance, equivalent time base distance Used to deal with the problem of unequal risks in different directions, equivalent time base distance It is based on the time dimension weighting of time distance and space-time distance. The specific expression is as follows: ; ; Where, is the coordinate of the observation point in the vehicle body coordinate system; is the vehicle width; is the vehicle length; The four vertices of the rectangle The corresponding equivalent time base distance, is the conversion parameter; A2.4. Combining temporal distance and time to obtain spatial-temporal distance, exist The space-time distance of moments Solve according to the following formula: ; In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, for The point closest to the observation point on the obstacle trajectory at that moment The equivalent time base distance from the observation point, The point closest to the observation point on the obstacle trajectory at that moment Recorded as , the time required for a vehicle to reach a designated location is called time-based distance, for Time trajectory scatter of axis position, is the weight of the time effect; Select the minimum value of the spatiotemporal distance of each scattering point as the point The space-time distance to the obstacle is calculated as follows: ; In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, .

[0009] As a preferred embodiment of the present invention, step 2 further includes the following steps: B1. Calculate the field strength at the road boundary. Road boundary field strength The calculation formula is as follows: ; in, is the road width; is the lateral position of the road boundary; is the lateral position of the observation point, is the coefficient to be determined; B2. Calculate the field strength of the solid road line. Field strength of solid road line The calculation formula is as follows: ; in, is the road width; is the lateral position of the solid line of the road; is the lateral position of the observation point, is the coefficient to be determined; B3. Calculate the field strength of the dotted line of the road. Road dotted line field strength The calculation formula is as follows: ; in, is the road width; is the lateral position of the road's imaginary line; is the lateral position of the observation point, is the coefficient to be determined; B4. Calculate the road boundary and road line risk field, Road boundary and road line risk areas The calculation formula is as follows: ; in, is the road boundary field strength; is the field strength of the solid road line; is the field strength of the dotted line of the road.

[0010] As a preferred embodiment of the present invention, step three also includes the following steps: Field strength in the interweaving area The calculation method is as follows: ; Where, The starting point of the off-ramp Axis coordinates, Starting position of the forced lane change area Axis coordinates, Add a spot to the vehicle location, The mandatory lane change area. and is a parameter to be determined, It is a collection of vehicles that are entering the ramp and intend to merge into the main line. It is a collection of vehicles on the main line that intend to get off the ramp. Gather for the remaining vehicles.

[0011] As a preferred embodiment of the present invention, step 4 further includes the following steps: The total risk field is considered as the superposition of multiple traffic element fields. The total risk field is obtained by adding the space-time risk field of obstacles, road boundaries and road line fields, and the risk field of on- and off-ramp on the expressway. The total risk of the field is calculated as follows: ; Where, Yes Total risk field, The obstacle point The applied field strength, Road boundaries and road line points The applied field strength, It is the intersection point Applied field strength.

[0012] As a preferred embodiment of the present invention, step five also includes the following steps: E1. Extract vehicle trajectory based on YOLOv8. E1.1 Vehicle detection module based on YOLOv8, The vehicle detection module based on YOLOv8 is used to directly predict the location, size and category of the target in a single network; E1.2, Deep-SORT-based target tracking module, The Deep-SORT-based target tracking module is used to describe the target's appearance information and distinguish different vehicles; E1.3, video stabilization processing, Use the After Effects system to pre-process the video data. By tracking multiple fixed buildings in the video, the rotation and displacement information of the entire video is extracted, and then the inverse motion is applied to eliminate the transformation information. E1.4, Model training and testing, Construct a dedicated dataset for vertical view and use the image annotation tool LabelImg for data annotation; E1.5, Coordinate system conversion, By obtaining the position information of road boundaries and road lines, converting them into coordinate information in the Frenet coordinate system, we can obtain the real-world expressway vehicle trajectory.

[0013] As a preferred embodiment of the present invention, step six also includes the following steps: F1. Determine the objective function, Based on the risk dynamic balance theory, the mean square error of the distance between the sample risk after decision making and the risk dynamic balance curve is used as the objective function of the calibration model. The specific formula is as follows: ; ; Where, is the objective function of the calibration model, that is, The distance between the risk after a sample decision and the acceptable risk space; For the The distance between the decision-making of each sample and the risk dynamic equilibrium curve; For the The risk value of a sample at a specific moment after the decision; and They are the upper limit risk and lower limit risk of the risk dynamic balance theory; is the number of samples; F2. Determine the upper limit risk of the risk dynamic balance theory and lower risk ; Calculate cap risk using a proprietary vertical viewpoint dataset and lower risk The vertical view perspective exclusive dataset includes high-risk dataset, low-risk dataset and decision dataset, which contains the trajectories of the observed vehicle and multiple surrounding vehicles. The high-risk dataset includes Includes collision time assessment, observing vehicles encountering high-risk scenarios, using 3 seconds as the TTC threshold, and high-risk datasets to determine the upper limit risk ,On the contrary, the low-risk dataset contains scenarios where the observed vehicle always remains safe according to the ,judgment of the TTC threshold, which is used to determine the optimal risk level. , upper limit risk and optimal risk level The calculation method is as follows, ; ; in To observe the level of risk experienced by the vehicle after making a decision; based on and , determine the lower limit risk , the calculation method is as follows, ; F3, using hybrid genetic algorithm and variable neighborhood search algorithm to solve, A hybrid genetic algorithm and a variable neighborhood search algorithm are used to solve the risk field parameters, and according to steps one to three, the driving risk assessment of autonomous vehicles in the expressway weaving area is realized.

[0014] The beneficial effects of this invention are as follows: By innovatively introducing the physical quantity of space-time distance, the two-dimensional spatial risk field is upgraded to a three-dimensional time-space risk field, effectively quantifying the impact of predicted trajectory and time effects on risk. A special geometric configuration field for on- and off-ramps in weaving areas is also constructed, which accurately restricts the operation of specific vehicles. Furthermore, based on the theory of dynamic risk balance and machine vision technologies such as YOLO, the model parameters are calibrated using real aerial data. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] By referring to the following description in conjunction with the accompanying drawings, and with a more complete understanding of the present invention, other objects and results of the present invention will become more clear and easy to understand. In the accompanying drawings: Figure 1 A schematic diagram of the vehicle geometric profile space distance of the present invention; Figure 2 A schematic diagram of a vehicle weaving field on an expressway on-ramp or off-ramp according to the present invention; Figure 3 This is a trajectory recognition effect diagram of the present invention; Figure 4 This is a flow chart for solving the hybrid genetic algorithm and variable neighborhood search algorithm of the present invention. DETAILED DESCRIPTION

[0016] See Figure 1-4 The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0017] An embodiment of the present invention provides a method for assessing driving risk of an autonomous driving vehicle in an expressway weaving area, comprising the following steps: Step 1: Modeling the space-time risk field of obstacles. A1. Construct a space-time coordinate system. A spatiotemporal three-dimensional coordinate system is established based on the Frenet coordinate system, and the direction of travel along the road is set as Axis, perpendicular to the road driving direction Axis, time is set to axis; A2. Calculate the time and space distance, Since the obstacle trajectory does not follow a specific function, the time-space distance cannot be calculated directly. exist Taking the time-space distance at the moment as an example, the time-space distance calculation method is used, including: obtaining the predicted trajectory of the obstacle, the spatial distance calculation method considering the geometric contour of the vehicle, the calculation method of the equivalent time-based distance considering the risk anisotropy, and obtaining the time-space distance; A2.1. Obtain the predicted trajectory of the obstacle. Obstacle trajectory scatter point collection Expressed as: ,in, express The space-time trajectory of the obstacle prediction at each moment, , , They are The position of the obstacle along the direction of road travel and the position perpendicular to the direction of road travel, that is, the position in the time and space coordinate system; A2.2. Calculate the spatial distance considering the vehicle's geometric outline. Treating obstacles as point masses with mass but no shape does not fully consider the impact of vehicle geometry, reducing the accuracy of risk quantification. To address this issue, a spatial distance solution method based on coordinate transformation is designed. In this method, the obstacle shape is regarded as a rectangle with fixed length and width values. This is because most obstacles during driving are interfering vehicles. because The axis does not change, so the establishment ignores The time-space coordinate system of the axis and the body coordinate system. In the body coordinate system, the center of mass of the vehicle is the origin and the direction of the vehicle is Axis, perpendicular to the direction of vehicle travel Axis, where, in the space-time coordinate system, point The coordinates are , the coordinates of the obstacle mass center are , in the vehicle body coordinate system, the four vertices of the rectangle The coordinates are: ,In addition, the vehicle's yaw angle , the speed of the vehicle , using the two steps of translation and rotation, we can get two types of coordinate systems, as follows: ; The formula for solving the spatial distance based on the vehicle body coordinate system is shown below; ; ; Where, is the coordinate of the observation point in the vehicle body coordinate system; is the coordinate of the observation point; is the transformed coordinate; is the vehicle width; is the vehicle length, The four vertices of the rectangle With the transformed observation point spatial distance; A2.3. Calculate the equivalent time-base distance considering risk anisotropy. Considering that the Euclidean distance in A2.2 has two major flaws: first, it assumes that obstacles pose equal risks to objects at different locations but equal distances, which obviously violates common sense; second, the Euclidean distance and time are not the same physical properties, and there is a disadvantage of being physically meaningless when weighting the time-space distance; To improve the above problems, the equivalent time-based distance is introduced. The concept of equivalent time base distance Indicates: When the point of the interactive object Directly in front of the obstacle course, the point of the interactive object When not directly in front of the vehicle's path, assume that the point of the interactive object Points with interactive objects With equal collision risk, the point from the vehicle to the interactive object The time distance is equal to the point from the vehicle to the interactive object Equivalent time distance, equivalent time base distance It is used to deal with the problem of unequal risks in different directions, and in terms of physical properties, the equivalent time base distance It is based on the time dimension weighting of time distance and space-time distance. The specific expression is as follows: ; ; Where, is the coordinate of the observation point in the vehicle body coordinate system; is the vehicle width; is the vehicle length; The four vertices of the rectangle The corresponding equivalent time base distance, is the conversion parameter; A2.4. Combining temporal distance and time to obtain spatial-temporal distance, exist The space-time distance of moments (i.e. the distance from the observation point to the vehicle trajectory scatter point) is solved according to the following formula, ; In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, for The point closest to the observation point on the obstacle trajectory at that moment The equivalent time base distance from the observation point, The point closest to the observation point on the obstacle trajectory at that moment Recorded as , the time required for a vehicle to reach a designated location is called time-based distance, for Time trajectory scatter of axis position, is the weight of the time effect; After obtaining the time-space distances of all scattering points, the minimum time-space distance of each scattering point is selected as the point The space-time distance to the obstacle is calculated as follows: ; In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, ; A3. Calculate the amount of charge taking into account its own and motion properties. In the electrostatic field, charge is a fundamental physical quantity that characterizes the electrical properties of an object (positive or negative) and its relationship with potential energy. Similarly, in the risk field, charge represents the comprehensive attributes of vehicle functions and driving behavior. Two representative parameters are selected from the comprehensive attributes of vehicle functions and driving behavior as key influencing factors of charge: braking performance and reaction time. The reasons for selecting these two parameters are: (1) Braking performance reflects the ability of the vehicle to avoid risks when faced with risks; (2) Reaction time, including driver response time and vehicle system delay, determines the driver and vehicle's perception and response efficiency to potential risks; (3) Charge is measured by combining the vehicle's own properties, motion properties, mass, and acceleration. Charge The calculation formula is defined as follows; ; Where, is the mass of the obstacle, is the speed of the vehicle, is the maximum deceleration of the vehicle, is the vehicle reaction time, which is the sum of the driver's reaction time and the vehicle system delay, and is the unknown coefficient, is an exponential function; A4. Calculate the dielectric constant taking into account anisotropy and environmental complexity. In electrostatic fields, the dielectric constant is a physical quantity that describes the ability of a material to enhance the charge storage capacity of a capacitor. It is defined as the product of the absolute dielectric constant in a vacuum and the relative dielectric constant in a given environment. In the context of risk fields, an analogy is established between the dielectric constant in risk fields and the dielectric constant in electrostatic fields. The dielectric constant in the risk field is defined as the product of the anisotropic attenuation coefficient and the environmental attenuation coefficient. The calculation formula is as follows: ; in, is the anisotropic attenuation coefficient, is the environmental complexity, which is calculated as follows: ; ; in, The obstacle's direction of travel and The angle between is the position vector of the obstacle pointing to the observation point, and are the speed and acceleration of the vehicle, is the ratio of the square difference of the speed in the weaving section to that of the normal road, is the ratio of the lane change ratio in the weaving section to that in the normal road, 、 and is the coefficient to be determined; A5. Calculate the spatial and temporal risk of obstacles. Obstacle space-time risk The calculation formula is as follows: ; in, is the space-time distance; is the charge; is the dielectric constant; Step 2: Modeling the road boundary and road line risk field. Road boundaries and lane markings impose restrictions on vehicles traveling normally. To avoid penalties associated with crossing lane markings, drivers constantly adjust the vehicle's lateral position throughout the driving process. Therefore, the risk field of road boundaries and lane markings varies with relative distance. The size of the risk field is only related to the line shape and lateral position of the road boundaries and lane markings. B1. Calculate the field strength at the road boundary. Road boundary field strength The calculation formula is as follows: ; in, is the road width; is the lateral position of the road boundary; is the lateral position of the observation point, is the coefficient to be determined; B2. Calculate the field strength of the solid road line. Field strength of solid road line The calculation formula is as follows: ; in, is the road width; is the lateral position of the solid line of the road; is the lateral position of the observation point, is the coefficient to be determined; B3. Calculate the field strength of the dotted line of the road. Road dotted line field strength The calculation formula is as follows: ; in, is the road width; is the lateral position of the road's imaginary line; is the lateral position of the observation point, is the coefficient to be determined; B4. Calculate the road boundary and road line risk field, Road boundary and road line risk areas The calculation formula is as follows: ; in, is the road boundary field strength; is the field strength of the solid road line; is the field strength of the dotted line of the road; Step 3: Modeling the risk field of on- and off-ramps on the expressway. In the weaving area of an expressway, vehicles on the ramp merging onto the main line must move from the acceleration lane to the inner lane within specified spatial constraints, while vehicles on the ramp off must move from the main line lane to the outer lane within specified spatial constraints. For a typical weaving area geometry, the forced lane change behavior on the ramps is considered an adaptive decision influenced by the risk field of the longitudinal acceleration and deceleration lanes. Field strength in the interweaving area The calculation method is as follows: ; Where, The starting point of the off-ramp Axis coordinates, Starting position of the forced lane change area Axis coordinates, Add a spot to the vehicle location, The mandatory lane change area. and is a parameter to be determined, It is a collection of vehicles that are entering the ramp and intend to merge into the main line. It is a collection of vehicles on the main line that intend to get off the ramp. Gather for the remaining vehicles; Step 4: Calculate the total risk field. The total risk field is considered as the superposition of multiple traffic element fields. The total risk field is obtained by adding the space-time risk field of obstacles, road boundaries and road line fields, and the risk field of on- and off-ramp on the expressway. The total risk of the field is calculated as follows: ; Where, Yes Total risk field, The obstacle point The applied field strength, Road boundaries and road line points The applied field strength, It is the intersection point applied field strength; Step 5: Real-world expressway vehicle trajectory extraction, E1. Extract vehicle trajectory based on YOLOv8. E1.1 Vehicle detection module based on YOLOv8 (vehicle detection algorithm), The YOLO family of algorithms can directly predict the location, size, and category of an object within a single network. Their high detection speed makes real-time detection easy. The YOLOv8 algorithm further improves detection accuracy and speed based on the YOLO family of algorithms, enabling more effective detection of vehicles in videos. The vehicle detection module based on YOLOv8 is used to directly predict the location, size and category of the target in a single network; E1.2, Deep-SORT-based target tracking module (target tracking algorithm), After vehicle detection, the results need to be combined with a multi-target tracking (MOT) algorithm to continuously track the vehicle's trajectory. Deep-SORT is an improvement on the SORT (Simple Online and Realtime Tracking) algorithm. It introduces deep learning features to describe the target's appearance, improving target matching accuracy and tracking stability. When vehicles are obscured or their appearance changes, Deep-SORT can use appearance features to better distinguish between different vehicles, reducing mismatches and target loss, and providing greater robustness in complex scenarios. The Deep-SORT-based target tracking module is used to describe the target's appearance information and distinguish different vehicles; E1.3, video stabilization processing, When drones shoot at high altitudes, they are affected by wind direction and air currents, causing the captured video to rotate and shift, inevitably resulting in slight jitter. This change in the video image can alter the vehicle's pixel coordinates, causing them to deviate from their actual position. To stabilize the video and extract the vehicle's trajectory, we use the After Effects system to pre-process the video data. By tracking multiple fixed buildings in the video, we extract the overall rotation and displacement information of the video. We then apply an inverse motion to eliminate this transformation information, thereby achieving video stabilization. E1.4, Model training and testing, Many existing object detection models rely on training datasets acquired from oblique perspectives, reflecting common shooting scenarios. Focusing on vehicle detection from a vertical bird's-eye view, while this perspective provides more comprehensive and unique vehicle information, existing models trained on oblique-perspective datasets are unable to adapt to these image characteristics, resulting in poor detection performance from this perspective. To address this issue, a dedicated vertical bird's-eye view dataset was constructed and annotated using the image annotation tool LabelImg. The selected altitude was approximately 350 meters, the minimum required to fully cover the intersecting area and help capture global features. The results of aerial vehicle recognition and tracking were also analyzed. E1.5, Coordinate system conversion, Due to limitations in the camera's battery life and human control, the shooting height and position of each video set cannot be completely consistent. As a result, the same position in different video sets may have discrepancies. Therefore, by obtaining the position information of road boundaries and road lines, and converting it into the Frenet coordinate system, we can obtain the real-world vehicle trajectory on the expressway. Step 6: Parameter calibration based on risk dynamic balance theory. F1. Determine the objective function, Based on the risk dynamic balance theory, the mean square error of the distance between the sample risk after decision making and the risk dynamic balance curve is used as the objective function of the calibration model. The specific formula is as follows: ; ; Where, is the objective function of the calibration model, that is, The distance between the risk after a sample decision and the acceptable risk space; For the The distance between the decision-making of each sample and the risk dynamic equilibrium curve; For the The risk value of a sample at a specific moment after the decision; and They are the upper limit risk and lower limit risk of the risk dynamic balance theory; is the number of samples; F2. Determine the upper limit risk of the risk dynamic balance theory and lower risk ; Calculate cap risk using a proprietary vertical viewpoint dataset and lower risk The vertical view perspective exclusive dataset includes high-risk dataset, low-risk dataset and decision dataset, which contains the trajectories of the observed vehicle (the vehicle) and multiple surrounding vehicles. The high-risk dataset includes Includes time-to-collision (TTC) assessment, observing vehicles encountering high-risk scenarios, using 3 seconds as the TTC threshold, and high-risk datasets to determine the upper limit of risk ,On the contrary, the low-risk dataset contains scenarios where the observed vehicle always remains safe according to the ,judgment of the TTC threshold, which is used to determine the optimal risk level. , upper limit risk and optimal risk level The calculation method is as follows, ; ; in To observe the level of risk experienced by the vehicle after making a decision; based on and , determine the lower limit risk , the calculation method is as follows, ; F3, using hybrid genetic algorithm and variable neighborhood search algorithm to solve, A hybrid genetic algorithm and a variable neighborhood search algorithm are used to solve the risk field parameters ( 、 、 、 ), and according to steps one to three, realize the driving risk assessment of autonomous driving vehicles in the expressway weaving area.

[0018] The above are merely specific embodiments of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A method for assessing driving risk of an autonomous vehicle in an expressway weaving zone, characterized in that: The following steps are involved: Step 1: Modeling the space-time risk field of obstacles, including: constructing a space-time coordinate system, calculating space-time distance, calculating the charge considering its own and motion properties, calculating the dielectric constant considering anisotropy and environmental complexity, and calculating the space-time risk of obstacles; Step 2: Modeling the risk field of road boundaries and road lines, including: calculating the field strength of road boundaries, calculating the field strength of road solid lines, calculating the field strength of road dashed lines, and calculating the risk field of road boundaries and road lines; Step 3: Modeling the risk field of on- and off-ramps on the expressway; Step 4: Calculate the total risk field; Step 5: Real-world freeway vehicle trajectory extraction, including: vehicle trajectory extraction based on YOLOv8, target tracking module based on Deep-SORT, video stabilization processing, model training and testing, and coordinate system conversion; Step 6: Parameter calibration based on the risk dynamic balance theory, including: determining the objective function, determining the upper and lower risk limits of the risk dynamic balance theory, and using a hybrid genetic algorithm and a variable neighborhood search algorithm to solve the problem and realize the driving risk assessment of autonomous vehicles in the expressway weaving area.

2. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step 1 also includes the following steps: A1. Construct a space-time coordinate system. A spatiotemporal three-dimensional coordinate system is established based on the Frenet coordinate system, and the direction of travel along the road is set as Axis, perpendicular to the road driving direction Axis, time is set to axis; A2. Calculate spatiotemporal distance, including: obtaining predicted obstacle trajectories, calculating spatial distances that take into account vehicle geometry, calculating equivalent time-based distances that take into account risk anisotropy, and obtaining spatiotemporal distances. A3. Calculate the amount of charge taking into account its own and motion properties. Two representative parameters are selected from the comprehensive attributes of vehicle functions and driving behavior as key influencing factors of charge: braking performance and reaction time; Charge The calculation formula is defined as follows; ； Where, is the mass of the obstacle, is the speed of the vehicle, is the maximum deceleration of the vehicle, is the vehicle reaction time, which is the sum of the driver's reaction time and the vehicle system delay, and is the unknown coefficient, is an exponential function; A4. Calculate the dielectric constant taking into account anisotropy and environmental complexity. The dielectric constant in the risk field is defined as the product of the anisotropic attenuation coefficient and the environmental attenuation coefficient. The calculation formula is as follows: ； in, is the anisotropic attenuation coefficient, is the environmental complexity, which is calculated as follows: ； ； in, The obstacle's direction of travel and The angle between is the position vector of the obstacle pointing to the observation point, and are the speed and acceleration of the vehicle, is the ratio of the square difference of the speed in the weaving section to that of the normal road, is the ratio of the lane change ratio in the weaving section to that in the normal road, 、 and is the coefficient to be determined; A5. Calculate the spatial and temporal risk of obstacles. Obstacle space-time risk The calculation formula is as follows: ； in, is the space-time distance; is the charge; is the dielectric constant.

3. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step A2 also includes the following steps: A2.

1. Obtain the predicted trajectory of the obstacle. Obstacle trajectory scatter point collection Expressed as: ,in, express The space-time trajectory of the obstacle prediction at each moment, , , They are The position of the obstacle along the direction of travel of the road and the position perpendicular to the direction of travel of the road; A2.

2. Calculate the spatial distance considering the vehicle's geometric outline. In the vehicle body coordinate system, the center of mass of the vehicle is taken as the origin and the direction of the vehicle is Axis, perpendicular to the direction of vehicle travel Axis, where, in the space-time coordinate system, point The coordinates are , the coordinates of the obstacle mass center are , in the vehicle body coordinate system, the four vertices of the rectangle The coordinates are: ,In addition, the vehicle's yaw angle , the speed of the vehicle , two types of coordinate systems are transformed using translation and rotation, as follows: ； Solving spatial distance based on vehicle body coordinate system The formula is as follows; ； ； Where, is the coordinate of the observation point in the vehicle body coordinate system; is the coordinate of the observation point; is the transformed coordinate; is the vehicle width; is the vehicle length, The four vertices of the rectangle With the transformed observation point spatial distance; A2.

3. Calculate the equivalent time-base distance considering risk anisotropy. Introducing equivalent time base distance , equivalent time base distance Indicates: When the point of the interactive object Directly in front of the obstacle course, the point of the interactive object When not directly in front of the vehicle's path, assume that the point of the interactive object Points with interactive objects With equal collision risk, the point from the vehicle to the interactive object The time distance is equal to the point from the vehicle to the interactive object Equivalent time distance, equivalent time base distance Used to deal with the problem of unequal risks in different directions, equivalent time base distance It is based on the time dimension weighting of time distance and space-time distance. The specific expression is as follows: ； ； Where, is the coordinate of the observation point in the vehicle body coordinate system; is the vehicle width; is the vehicle length; The four vertices of the rectangle The corresponding equivalent time base distance, is the conversion parameter; A2.

4. Combining temporal distance and time to obtain spatial-temporal distance, exist The space-time distance of moments Solve according to the following formula: ； In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, for The point closest to the observation point on the obstacle trajectory at that moment The equivalent time base distance from the observation point, The point closest to the observation point on the obstacle trajectory at that moment Recorded as , the time required for a vehicle to reach a designated location is called time-based distance, for Time trajectory scatter of axis position, is the weight of the time effect; Select the minimum value of the spatiotemporal distance of each scattering point as the point The space-time distance to the obstacle is calculated as follows: ； In the formula for The temporal and spatial distance between the observation point and the obstacle at any moment, .

4. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step 2 also includes the following steps: B1. Calculate the field strength at the road boundary. Road boundary field strength The calculation formula is as follows: ； in, is the road width; is the lateral position of the road boundary; is the lateral position of the observation point, is the coefficient to be determined; B2. Calculate the field strength of the solid road line. Field strength of solid road line The calculation formula is as follows: ； in, is the road width; is the lateral position of the solid line of the road; is the lateral position of the observation point, is the coefficient to be determined; B3. Calculate the field strength of the dotted line of the road. Road dotted line field strength The calculation formula is as follows: ； in, is the road width; is the lateral position of the road's imaginary line; is the lateral position of the observation point, is the coefficient to be determined; B4. Calculate the road boundary and road line risk field, Road boundary and road line risk areas The calculation formula is as follows: ； in, is the road boundary field strength; is the field strength of the solid road line; is the field strength of the dotted line of the road.

5. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step three also includes the following steps: Field strength in the interweaving area The calculation method is as follows: ； Where, The starting point of the off-ramp Axis coordinates, Starting position of the forced lane change area Axis coordinates, Add a spot to the vehicle location, The mandatory lane change area. and is a parameter to be determined, It is a collection of vehicles that are entering the ramp and intend to merge into the main line. It is a collection of vehicles on the main line that intend to get off the ramp. Gather for the remaining vehicles.

6. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step 4 also includes the following steps: The total risk field is considered as the superposition of multiple traffic element fields. The total risk field is obtained by adding the space-time risk field of obstacles, road boundaries and road line fields, and the risk field of on- and off-ramp on the expressway. The total risk of the field is calculated as follows: ； Where, Yes Total risk field, The obstacle point The applied field strength, Road boundaries and road line points The applied field strength, It is the intersection point Applied field strength.

7. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step five also includes the following steps: E1. Extract vehicle trajectory based on YOLOv8. E1.1 Vehicle detection module based on YOLOv8, The vehicle detection module based on YOLOv8 is used to directly predict the location, size and category of the target in a single network; E1.2, Deep-SORT-based target tracking module, The Deep-SORT-based target tracking module is used to describe the target's appearance information and distinguish different vehicles; E1.3, video stabilization processing, Use the After Effects system to pre-process the video data. By tracking multiple fixed buildings in the video, the rotation and displacement information of the entire video is extracted, and then the inverse motion is applied to eliminate the transformation information. E1.4, Model training and testing, Construct a dedicated dataset for vertical view and use the image annotation tool LabelImg for data annotation; E1.5, Coordinate system conversion, By obtaining the position information of road boundaries and road lines, converting them into coordinate information in the Frenet coordinate system, we can obtain the real-world expressway vehicle trajectory.

8. The method for risk assessment of autonomous driving vehicles in an expressway weaving zone according to claim 1 is characterized in that: Step six also includes the following steps: F1. Determine the objective function, Based on the risk dynamic balance theory, the mean square error of the distance between the sample risk after decision making and the risk dynamic balance curve is used as the objective function of the calibration model. The specific formula is as follows: ； ； Where, is the objective function of the calibration model, that is, The distance between the risk after a sample decision and the acceptable risk space; For the The distance between the decision-making of each sample and the risk dynamic equilibrium curve; For the The risk value of a sample at a specific moment after the decision; and They are the upper limit risk and lower limit risk of the risk dynamic balance theory; is the number of samples; F2. Determine the upper limit risk of the risk dynamic balance theory and lower risk ; Calculate cap risk using a proprietary vertical viewpoint dataset and lower risk The vertical view perspective exclusive dataset includes high-risk dataset, low-risk dataset and decision dataset, which contains the trajectories of the observed vehicle and multiple surrounding vehicles. The high-risk dataset includes Includes collision time assessment, observing vehicles encountering high-risk scenarios, using 3 seconds as the TTC threshold, and high-risk datasets to determine the upper limit risk ,On the contrary, the low-risk dataset contains scenarios where the observed vehicle always remains safe according to the ,judgment of the TTC threshold, which is used to determine the optimal risk level. , upper limit risk and optimal risk level The calculation method is as follows, ； ； in To observe the level of risk experienced by the vehicle after making a decision; based on and , determine the lower limit risk , the calculation method is as follows, ； F3, using hybrid genetic algorithm and variable neighborhood search algorithm to solve, A hybrid genetic algorithm and a variable neighborhood search algorithm are used to solve the risk field parameters, and according to steps one to three, the driving risk assessment of autonomous vehicles in the expressway weaving area is realized.

Images