Power system operators need new, efficient operational tools to use the flexibility of distributed resources and deal with the challenges of highly uncertain and variable power systems. Transmission system operators can consider the available flexibility in distribution systems (DSs) without breaching the DS constraints through flexibility areas. However, there is an absence of open-source packages for flexibility area estimation. This paper introduces TensorConvolutionPlus, a user-friendly Python-based package for flexibility area estimation. The main features of TensorConvolutionPlus include estimating flexibility areas using the TensorConvolution+ algorithm, the power flow-based algorithm, an exhaustive PF-based algorithm, and an optimal power flow-based algorithm. Additional features include adapting flexibility area estimations from different operating conditions and including flexibility service providers offering discrete setpoints of flexibility. The TensorConvolutionPlus package facilitates a broader adaptation of flexibility estimation algorithms by system operators and power system researchers.

Power system operators require advanced applications in the control centers to tackle increasingly variable power transfers effectively. One urgently needed application concerns estimating the feasible available aggregated flexibility from a power system network, which can be effectively deployed to mitigate issues in interconnected networks. This paper proposes the TensorConvolution+ algorithm to address the above application. Unlike related literature approaches, TensorConvolution+ estimates the density of feasible flexibility combinations to reach a new operating point within the p-q flexibility area. This density can improve the decision-making of system operators for efficient and safe flexibility deployment. The proposed algorithm applies to radial and meshed networks, is adaptable to new operational conditions, and can consider scenarios with disconnected flexibility areas. Using convolutions and tensors, the algorithm efficiently aggregates the combinations of flexibility providers' adjustable power output that can occur for each flexibility area set point. Simulations on the meshed Oberrhein and radial CIGRE test networks illustrate the effectiveness of TensorConvolution+ for flexibility estimation with high numerical confidence and a minor computing effort. Additional simulations highlight how system operators can interpret the estimated density of feasible flexibility combinations for decision-making purposes, the algorithm's capability to estimate disconnected flexibility areas, and adapt to new operating conditions.

Power electronic interfaced devices progressively enable the increasing provision of flexible operational actions in distribution networks. The feasible flexibility these devices can effectively provide requires estimation and quantification so the network operators can plan operations close to real- time. Existing approaches estimating the distribution network flexibility require the full observability of the system, meaning topological and state knowledge. However, the assumption of full observability is unrealistic and represents a barrier to system operators' adaptation. This paper proposes a definition of the distribution network flexibility problem that considers the limited observability in real-time operation. A critical review and assessment of the most prominent approaches are done, based on the proposed definition. This assessment showcases the limitations and benefits of existing approaches for estimating flexibility with low observability. A case study on the CIGRE MV distribution system highlights the drawbacks brought by low observability.