Causal Operator Discovery in Partial Differential Equations via Counterfactual Physics-Informed Neural Networks