The integration of stochastic renewable energy sources into energy systems presents challenges due to flexibility insufficiency and uncertainty modeling inaccuracy. To address these challenges, this paper explores a high-dimensional covariance matrix approach based on stochastic differential equations (SDEs) for energy and reserve co-dispatch in integrated electricity and heating systems. By capturing the long-term correlation of wind power forecast errors, the SDE model simultaneously models point forecasts and high-dimensional positive-definite covariance matrices, resulting in a compact uncertainty set. To ensure computational efficiency, we employ duality theory to transform the stochastic co-dispatch problem with the high-dimensional uncertainty set into deterministic convex programming. Additionally, a virtual heat storage model is introduced for the district heating network, leveraging the thermal inertia of the network to enhance reserve capacity. Simulations validate the calibration and sharpness of the proposed high-dimensional uncertainty set. The results demonstrate that co-dispatch using the SDE-based uncertainty set reduces operational costs by 4.9% compared to a set that does not consider the long-term correlation of wind power, under a 95% coverage rate. Moreover, the storage capabilities of district heating pipelines and smart buildings provide cost-free reserves, saving the operational costs of the two test systems by 2.8% and 1.3%, respectively.