I'm not sure about published references to this specific problem, but I'm pretty sure it can be solved in polynomial time via a reduction to minimum weight perfect matching, as follows.

Replace each vertex of degree $d$ by a complete bipartite graph $K_{d-2,d}$ where each of the original edges incident to the vertex becomes incident to one of the vertices on the $d$-vertex side of the bipartition. Set all edge weights in this complete bipartite graph to zero. Next, add a complete graph (with zero edge weights) connecting all of the vertices in the whole graph that are on the $(d-2)$-vertex sides of their bipartitions.

For each vertex $v$ in the original graph, a perfect matching in this modified graph has to include at least two of the original graph's edges incident to $v$, because at least two of the vertices on the $d$-vertex side of the complete bipartite graph that replaces $v$ are not matched within that complete bipartite graph. Because all the other edges have cost zero, the cost of the perfect matching is the same as the cost of the solution it leads to in the original graph.

On the other hand, whenever one has a subgraph of the original graph that includes at least two edges at each vertex, it can be completed to a perfect matching at no extra cost by matching the remaining vertices on the $d$-vertex sides of their bipartition to vertices on the $(d-2)$-vertex sides, and then using the complete graph edges to complete the matching among any remaining unused vertices on the $(d-2)$-vertex sides.

Therefore, the cost of the minimum weight perfect matching in this graph is the same as the cost of the optimal solution to your problem.

Added later: something like this appears to be in "An efficient reduction technique for degree-constrained subgraph and bidirected network flow problems", Hal Gabow, STOC 1983, doi:10.1145/800061.808776